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ABSTRACT 


This  contract  supported  my^work  toward  a  monograph  on  semiconductor 
allovs.  n  The  activities  supported  included  a  one-year  sabbatical 
leave  to  Stanford  University,  the  collection  of  information  for  the 
book,  and  the  development  of  theoretical  interpretations.  All  these 
activities  have  been  accomplished,  and  the  monograph  is  expected  to 
be  completed  in  1986. -  The  chapters  to  be  included  are:  Crystal 
structure  and  bonding,  Alloy  statistics  and  phase  diagram.  Electronic 
structure,  Mechanical  properties  and  lattice  vibration,  Defect  states 
and  structures,  Alloy  scattering  effects,  and  Material  engineering  -- 
alloy  and  superlattice .  A  copy  of  the  monograph,  when  published, 
will  be  forwarded  to  ON R  to  acknowledge  the  support. 
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ACCOMPLISHMENTS 


This  contract  supported  my  work  toward  a  monograph  on 
semiconductor  alloys.  This  a  very  fast  growing  area  of  research. 
Although  we  have  collected  an  extensive  data  base  for  the  book, 
new  data  are  pouring  in  and  are  changing  every  day.  Theoretical 
work  has  also  been  increasing  substantially  in  past  two  years;  yet 
the  interpretations  for  many  alloy  topics  are  still  very 
rudimentary.  Among  the  most  striking  discoveries  are  the 
bi-modal  bond  length  distributions  in  pseudo-binary  alloys,  the 
long-range  ordering  of  epitaxially  grown  alloys  such  as  GaAlAs , 
and  new  results  for  the  band  offsets.  Against  this  background,  we 
had  to  modify  our  work  constantly  in  order  to  prevent  the  book 
from  becoming  obsolete.  As  a  result,  we  have  further  extended  the 
deadline  of  the  book  to  the  end  of  1986.  However,  all  the 
activities  supported  by  this  contract  have  been  accomplished, 
which  included  a  one-year  sabbatical  visit  to  Stanford  University, 
the  collection  of  the  data  base  for  the  book,  and  the  development 
of  theoretical  interpretations.  Instead  of  further  extending  the 
expiration  date  of  the  contract  to  coincide  with  the  deadline  of 
the  book,  I  should  conclude  the  contract  now.  A  copy  of  this 
monograph,  when  published,  will  be  forwarded  to  ONR  to  acknowledge 
the  support.  Below  is  an  outline  of  the  contents: 


Chapter  One:  Crystal  Structure  and  Bonding 


1.  Crystal  structures  of  important  semiconductors,  Tables: 
Lattice  parameters,  Atomic  term  values  and  covalent  radii. 


2.  Tetrahedral  bonding:  covalent  and  ionic  nature,  Harrison's 
model  and  modified  version,  first-principles  theory,  Table: 
experimental  and  theoretical  cohesive  energies. 


3.  Structure  of  pseudo-binary  alloys  revealed  from  EXAFS  and 
theoretical  models,  Table:  Alloy  bond  lengths. 


4.  Other  alloy  structures,  clustering,  and  long-range  order. 
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Chapter  Two:  Alloy  Statistics  and  Phase  Diagram 
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1.  Mixing  free  energy,  enthalpy  and  entropy,  miscibility  gap, 
spinodal  decomposition,  critical  temperature.  Table: 
Thermodynamics  data  for  binary  compounds. 
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2.  Statistical  models:  ideal  alloy,  zeroth  approximation, 
quasi -chemical  approximation  (QCA) ,  cluster  variational 
method  (CVM) .  Tables:  Comparison  of  theoretical  Tc . 
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3.  Models  for  mixing  enthalpies,  Table:  Estimated  and  measured 
values . 
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h .  Phase  diagram,  Vieland  method,  Binary  liquidus,  Ternary 
liquidus-solidus  curves. 


.  Generalized  QCA. 

6.  The  16-bond  clusters  and  non-random  distributions. 

Chapter  Three:  Electronic  Structure 

1.  Elements  of  pure -crystal  band  structure. 

3.  Alloys  as  means  for  tailoring  the  band  gaps. 

4.  Bowing  of  alloy  band  gaps,  Table:  experimental  and  theoretical 
bowing  parameters. 

5.  The  coherent  potential  approximation  and  the  realistic 
tight-binding  theory. 

6.  CPA  band  structures  for  SixGe1.x,  and  III-V  and  II -VI 
pseudo-binary  alloys. 

7.  Effects  of  non-random  distribution. 

8.  Magnetic  semiconductor  alloys. 

Chapter  Four:  Mechanical  Properties  and  Lattice  Vibration 

1.  Elasticity  of  pure  semiconductors,  Table:  Bulk  modulus  and 
stiff  coefficients. 

2 .  Elastic  models . 

3.  Hardness,  Table:  Measured  hardness  of  semiconductors. 

t.  Pure  semiconductor  vibrational  bands,  Table:  LO  and  TO  modes. 
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Abstract 

The  objective  of  this  paper  Is  to  identify  the  principal  microscopic 
phenomena  controlling  dislocation  denaltlea  in  bulk  grown  semiconductors. 

Then,  based  on  this  understanding,  a  strategy  for  selecting  materials  to  reduce 
dislocation  densities  is  offered.  The  relevant  quantities  are  calculated  from 
an  extension  of  Harrison's  bonding  theory,  which,  with  our  improved  accuracy 
relates  properties  of  the  solids  to  the  constituent  atoms’  valence  electron 
energy  states  and  wave  functions.  We  report  on  the  alloy  composition  varia¬ 
tion  of  bond  energies,  bond  lengths,  charge  redistribution  as»ng  constituents, 
vacancy  formation  energies,  dislocation  energies,  and  hardness.  Several 
III-V  and  II-VI  compound  semiconductors  are  treated  including,  GaAs,  GalnAs, 
HgCdTe,  and  ZnHgTe  . 


(a)  This  work  was  supported  In  part  by  AFOSR  Contract  F49620-81-K0012  and 
DARPA  Contract  MDA-903-83-C-0108 . 

<  b  >  Stanford  Aacherman  Professor  of  Engineering 
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Introduction 


According  to  a  currently  accaptad  nodal  (1)  of  the  mechsnlsm  leading  to 
dislocations  In  bulk-grown  material  at  a  given  temperature  and  tapersture 
gradient  behind  the  growth  front;  their  density  Is  governed  by: 

•  The  vacancy  density  that  depends  on  the  formation  energies, 

•  The  competition  between  vacancy  annealing  rates  and  vacancy 
Interaction  caused  clusters, 

•  The  condensation  rate  of  these  vacancy  clusters  Into  dislocation 
loops ,  and 

•  The  subsequent  growth  rate  of  these  loops. 

The  objective  of  this  paper  Is  to  Identify  some  of  the  principal  microscopic 
phenomena  controlling  these  features  In  semiconductors  compounds  and  their 
pseudobinary  alloys.  If  any  of  the  foregoing  steps  can  be  Inhibited,  then 
there  will  be  fewer  dislocations  In  bulk-grown  crystals. 

A  model  of  the  bonding  of  tetrahedrally  coordinated  semiconductor  com¬ 
pounds  due  to  Harrison  underlies  this  work  (2).  He  has  derived  expressions 
for  the  bond  energies  and  strain  coefficients  of  pure  semiconductor  compounds 
In  terms  of  the  constituent  atom's  valence  state  wave  functions  and  energies. 

There  are  four  contributions  to  the  bond  energies: 

a  A  covalent  energy  arising  from  the  Interaction  between  sp^  hybrids 
on  adjacent  sites,  which,  according  to  a  universal  rule  deduced 
by  Harrison,  Is  related  to  the  Inverse  square  of  the  bond  length,  d, 

•  An  Ionic  energy  which  Is  proportional  to  the  energy  difference 
between  the  spd  hybrid  energies  of  the  anion  and  cation, 

•  A  metallic  energy  arising  from  the  Interaction  between  filled 
bonding  and  unfilled  antibonding  states  on  adjacent  bonds,  and 

-4 

e  An  overlap  repulsion  energy  which  Is  taken  to  vary  as  d 

The  shear  coefficients,  which  play  an  important  role  in  dislocation  energies, 
sre  shown  to  vary  for  covalent  solids  as  the  covalent  energy  per  cell  unit 
volume  or  as  d-*.  The  Ionic  energy  arising  as  It  does  from  coulomb  inter¬ 
actions,  la  Insensitive  to  bond  angle  distortions  and  depends  only  on  bond 
lengths.  Thus,  in  the  Ionic  semiconductors  the  shear  coefficients  are  reduced 
from  those  of  equivalent  bond  length  covalent  compounds,  but  also  the  power 
law  dependence  on  bond  length  increases  to  d“^  for  the  extreme  Ionic  limit  (3). 
The  metallic  Interaction  is  bond  angle  sensitive  and  contributes  to  the  d”** 
power  lav. 

We  have  modified  Harrison's  theory  in  several  ways  (3)  (4).  Pirat, 
instead  of  using  atomic  term  values  as  the  input  energies  in  the  bonding 
calculation*  we  use  correlated  atomic  energies  that  take  account  of  all  the 
atomic  electronic  energy  changes  associated  with  a  state  change.  In  our  prior 
publications  these  energies  were  calculated  from  pseudopotentials  (3).  In 
this  paper  we  have  modified  the  procedure  and  taken  the  ground  state  energies 
from  free  atom  ionization  experiments,  and  onlv  calculate  the  excitation 
energies  from  the  pseudopotentials.  Second,  the  Harrison's  theory  has  been 
extended  to  alloys  in  which  each  type  of  the  Individual  bond's  energy  and 
length  changes  and  the  net  alloy  substitution  energies  are  calculated.  A 
theory  of  dislocation  energies  and  hardness  of  semiconductors  has  been  advanced. (4) 
In  this  paper  the  second  item  will  be  emphasized  since  item  one  is  less 
relevant  to  this  conference  and  these  will  be  published  elsewhere. 


Bond  Us|th  and  Imio  Chant**  In  an  Alloy 


Single  Impurity 

Start  by  considering  tha  simplest  defect,  one  cation  of  a  host  sealcon- 
ductor  compound  AC  la  replaced  by  an  looelactranlc  lapurlty  I,  an  akatchod 
In  Fig.  X.  Thl*  will  be  generalized  to  an  AB^.-C,  alloy  later  ana  the  cation 
subatltutad  caac  Aj.^B^C  follow*  by  ayaetry.  To  stud?  the  principal  effect 

of  cor*  lattice  dlatortlona,  we  allow 
the  poeltlona  of  the  first  ahell  of 
atoaa  labeled  by  bond  length  dj  to  ame 
but  fix  the  atoms  of  the  aecond  neigh¬ 
bor*  and  beyond  at  their  pure-cryatal 
poaltlana.  Eventually,  the  affects  of 
long  range  strain  fields,  faneratnd  by 
a  point  distortion,  nust  be  added  to 
this  core  distortion.  Here  we  allow 
only  three  different  bond  lengths 
d!  -  (l-6)d,  d2  -  (1  +  2/  36  +  62)1/2d, 
and  d3  •  d.  where  d  is  the  bond  length 
of  the  pure  host  AC  compound ,  and  6  is 
a  scale  parameter. 

There  are  four  bonds  with  energy 
E],  twelve  with  energy  E2,  and  36  with 
energy  Ey  These  energies  can  be 
written  in  dosed  fora  in  the  Baxrlaoo 
theory.  The  energies  Ej  and  E2  differ 
fron  the  host  bond  energy  EQ  because 
of  bond  length  and  angle  distortions 
as  well  as  difference  in  the  ionic 
energy  of  I  and  C  atou  induced  charge 
transfer.  The  energy  E}  differs  fron 
E„  only  by  the  charge  transfer  coupled 
through  the  metallization  terns.  There 
are  two  interesting  energies  we  can  calculate.  The  first  is  the  energy  d» 
required  to  substitute  the  I  a ton  for  a  C  atoo,  l.e.  bring  a  free  I  atom  fron 

infinity  and  replace  a  C  atoa  that  Is  taken  fron  the  crystal  and  removed  to 

Infinity. 

d,  *  4  d  Ej  ♦  12  d  Ej  +  36  d  Ej  -  (tc  -  Cj)  (1) 

where  dE,  ■  I,  ■  I  ,  ]  ■  1,2,3,  and  «c  and  «j  are  the  free  aton  valence  elec¬ 

tron  binning  energies.  The  second  is  the  bond  energy  chenge  of  the  lapurlty 
AI  bend  relative  to  the  bond  energy  of  a  pure  A1  compound,  denoted 

1 4Ei  '  i  (ti '  «c)  (2) 

This  energy  tells  us  If  the  AI  bond  Is  stabilized  (dj  <  0)  or  destabilized 
(d|  >  0)  when  It  Is  In  an  AC  host.  If  we  define  dfc  •*  the  bond  energy  dif¬ 
ference  between  the  LA  and  CA  bonds  each  In  their  respective  pure  crystals 

t ^  «  BE(1A)  -  BE(CA)  (3) 


Schematic  representation  of  the 
bonds  around  an  atom  of  interest 
(designated  by  the  square).  The 
first,  second,  and  third  neigh¬ 
bor  bonds  are  designated  d, ,  d., 
sod  dj 


then  one  can  write 


4*  ■  4Ab  +  v 


(4) 


where  d  .  1*  extra  energy  difference  caused  by  strain  and  charge  transfer. 


Tha  enargyea  E.  and  band  lengths  d.,  ]  <  1,  2,3,  ara  calculated  by  mini¬ 
mizing  4,.  If  ws  lit  dj  s  d0  (1  -  ^,)  Ware  dj  la  tha  bond  length  in  the 
pure  AI lattice,  than  the  approximate  expression 
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_  (5) 

4-  (nail  tana  due  to  charge  traaafar) 


can  be  deduced.  Tha  bulk  nodulll  B  of  the  boat  and  B.  of  the  impurity  lattices, 
and  a  shear  coefficient  Cxx  •  of  tha  boat  appear  in  tha  expraaalon.  While 
this  expression  Is  approximate,  the  effecte  of  the  serious  types  of  strains 
can  be  visualised  using  it.  Our  detailed  calculations  are  done  using  the  full 
theory.  If  one  fixes  the  bond  lengths  *  dj  -  dj  ■  dQ,  then  6  -  and  the 
strain  energy  la  large  and  positive.  If  one  allows  only  bond  length  dis¬ 
tortions  (Cxi  ~^12  “  O)  then  i  is  reduced  and  the  bonds  tend  to  adjust  so 
d{  »  dj.  This  reduces  ASc  considerably  fron  the  undlstorted  lattice  case. 

It  dj  >  d„  then  both  bonds  are  stretched  somewhat,  dt  -  dj  >  0  and  dj  -  dD  >  0. 
However,  this  configuration  produces  large  bond  angle  distortions.  When  the 
shear  coefficients  are  turned  on  the  lattice  relaxation  is  modified  and  the 
bond  lengths  cannot  adjust  as  nuch,  so  di  differs  from  d.  by  a  larger  amount. 

The  net  result  la  that  i,c  Is  Increased  since  one  pays  the  price  of  strain 
energy  either  as  a  bond  length  or  bond  angle  distortion. 

If  we  calculate  6  from  the  full  theory,  then  we  predict  values  that  are 
too  large.  This  occurs  because  the  theory  predicts  Cxx  -  Cxx  which  agree  well 
with  experiment  but  it  predicts  bulk  nodulll  with  the  proper  treads  fron  one 
compound  to  another  but  which  are  about  a  factor  of  2.4  too  small.  If  we  use 
the  experimental  values  of  the  strain  coefficients  in  Eg.  (5)  then  good  agree¬ 
ment  is  found  with  the  experiments  (5)  on  and  ZnSex-sTex.  The 

results  are  quoted  in  Table  I. 

Table  I.  Bond  Lengths  in  A  for  Impurities  in  the  Indicated  Host 

Ga  in  InAa _ In  In  CaAa _ Se  In  ZnTe _ Te  In  ZnSe 


Experiment  ^ 

2.487 

2.587 

2.496 

2.595 

Eq. (5)  and 
experimental  B&C 

2.499 

2.547 

2.510 

2.570 

FU11  Theory 

2.538 

2.518 

2.570 

2.512 

Alloy 

In  an  A,  ^BjC  alloy  the  four  cations  around  a  given  C  anion  can  be 
arranged  in  five  different  configurations,  denoted  by  A(4-TJ),  B(T|),  7)  •  0,1, 
2,3,4.  An  A(2)  B(2)  configuration,  for  example,  is  one  in  which  the  C  atom 
has  two  A  atom  and  two  B  atom  neighbors.  It  is  possible  to  solve  the  full 
alloy  problem  for  large  clusters,  but  for  now  we  have  restricted  the  cluster 
to  near  neighbors  only,  and  in  Fig.  1  the  central  atom  is  now  an  anion  C  and 
dj  ■  dj  -  deff  are  taken  to  be  an  effective  medium  bond  length  that  is  deter¬ 
mined  self-conslstently.  The  different  bonds  of  type  1  no  longer  have  the 
same  length.  The  procedure  is  as  follows:  First  assume  a  value  for  deff, 
say  the  virtual  crystal  bond  length  (1-x)  d  (AC)  +  xd  (BC) .  Next  calculate 
Che  dx  values  for  Che  various  conf iguradons  try  minimizing  the  energy  of  the 
configuration.  Then,  configuration  average  the  various  d  values  to  find  a 
new  deff.  Then  iterate  this  procedure  until  it  converges. 
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Figure  2:  Bond  length  (a)  and  energy  (b)  of  Ga^In^Aa 
aa  functions  of  x.  The  dashed  curses  are  the  average  values 
for  the  designated  bonds,  and  the  solid  curved  are  the  alloy 

averages . 

Because  the  bond  angle  distortion  terms  are  imphyslcally  large  we  have  left 
them  out  of  the  present  calculation.  Bence  bond  length  shifts  will  be  some¬ 
what  too  small. 

Baaulta  for  the  bond  energy  and  bond  length  shifts  with  concentration  are 
given  for  Ca1_sIn,As  In  Fig.  2  and  Bj.^ln^F  in  Fig.  }  as  examples  of  two 
different  behavior  patterns.  The  predicted  trends  for  GalnAs  agree  with 
experiment  but  the  bond  length  changes  are  too  small  as  expected  with  the  bond 
angle  distortion  terms  absent.  Notice  that  the  longer  bond  In  this  case  de¬ 
creases  In  length  and  the  shorter  ooe  increases  as  expected.  However,  in  the 
8InJ>  case  the  charge  shift  terms  are  so  large  that  the  longer  InP  bond  has  a 
nlnlm»  as  a  fucntlon  of  composition  rather  than  a  mono  tonic  decrease.  Be¬ 
cause  of  the  large  bond  length  difference  between  BF  and  InP  there  la  probably 
a  misaability  gap  In  this  alloy  that  prevents  these  compositions  from  being 
prepared.  However,  one  nay  be  able  to  examine  the  variation  of  the  anomalous 
behavior  of  the  InP  bond  la  a  BP  boat  (x  «  1). 


The  Hgi_zCdLTe  system  Is  completely  anamolous.  The  bond  lengths  of  CdTe 
(2.80V  A)  and  HgTe  (2.797- 1)  are  nearly  the  sane  by  an  accident.  CdTe  bonds 
are  dominated  by  covalent  and  Ionic  Interactions,  while  HgTe  is  more  covalent 
and  the  aetalllc  terns  are  Important.  The  nix  of  Interactions  In  the  two 
cases  leads  to  the  s asm  bond  lengths.  When  an  alloy  Is  made  the  charge  shift 
terms  dominate,  and  they  cause  the  longer  CdTe  bond  to  become  still  longer  and 
the  HgTe  to  contract  by  aaounts  that  are  large  compared  to  the  pure  crystal 
difference.  Moreover,  the  already  weak  HgTe  hat  Its  bond  energy  reduced  still 
m>re  by  the  presence  of  Cd.  Since  the  neltlng  point  of  HgTe  Increases  as  Cd 
Is  added  and  the  local  strength  of  HgTe  bonds  adjacent  to  a  Cd  decrease, 
vacancy  densities  will  Increase.  All  these  and  other  observed  phenomena  are 
predicted  by  tbe  theory.  A  complete  catalog  of  results  requires  more  pages 
than  we  have  been  allocated  in  this  article,  but  we  have  tried  to  offer  a 
representative  group. 
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(a)  0>) 

Figure  3:  Bond  length  (*)  and  energy  (b>  of  Bj_xIoxl>  aa 
functions  of  a.  The  dashed  curves  are  the  average  values 
for  the  designated  bonds  and  the  solid  curves  are  the  alloy 

averages. 

Conclusions 

We  have  deaonstrated  In  agrecnent  with  experiment  that  dislocation  energies 
and  hardness  of  scale onduc tors  are  proportional  to  the  shear  coefficients  which 
vary  roughly  aa  d~*.  Thus  to  decrease  dislocations  In  a  given  sealconductor 
one  should  Inhibit  their  formation  by  Introducing  sene  means  of  shortening 
bonds.  This  should  be  acconpllshed  without  also  decreasing  the  ratio  of  the 
vacancy  format  Ion  energy  to  the  aelt  temperature.  While  detailed  calculations 
to  support  the  following  contentions  are  still  Incomplete,  a  set  of  criteria 
on  an  impurity  (denoted  I)  In  a  host  semiconductor  (denoted  ■)  chat  are  likely 
to  satisfy  these  conditions  are:  the  bond  length  of  the  impurity  Is  smaller 
than  that  of  the  host  dj  <  dH,  the  smaller  bond  energy  Is  stabilised  In  the 
alloy  min  |ivEbH,  AEb]}  <  °*  In  these  clrcuastances,  Che  average  hood  length 
will  shrine  and  In  tne  vicinity  of  each  Impurity  the  four  surrounding  bonds 
nearly  have  the  length  of  the  Impurity  but  are  stretched  slightly.  The  next 
neighboring  host  bonds  are  also  stretched.  The  net  affect  la  an  arrangement 
which  la  more  rigid  than  the  unperturbed  lattice  and  coossquently  the  local 
shear  strain  energy  Increases,  causing  dislocation  anarglas  to  increase. 

Exaaples  of  this  case  are  Bj.^Ca^As,  and  Zn,  gRgxTe.  loth  B  In  GaAs  (6)  and 
Zn  In  HgTe  (7)  have  proven  to  be  effective  In  reducing  dislocation  densities. 
Unfortunately,  B  has  a  low  solubility  tuCaAs  and  it  Is  not  clear  that  enough 
can  be  gotten  in  to  mane  it  easy  to  prepare  dislocation-free  material. 

A  second  case  where  an  improvement  occurs  Is  If  dj  >  dB,  and  again 
■in  |^EbH,  aEbll  >  0.  In  this  esse  esch  of  the  four  longer  Impurity  atom 
bonds  are  conpressed  by  the  surrounding  host  bonds  and,  more  Importantly,  the 
twelve  next  neighbor  host  bonds  are  also  expressed.  Once  again  for  small 
Impurity  concentrations  the  net  dislocation  energy  should  be  Increased.  How¬ 
ever  ,  In  this  case  the  effect  Is  competing  against  a  nst  bond  lengthening 
trend  of  the  alloy  which  Is  tending  to  nake  It  less  rigid  against  a  shear,  so 
at  higher  concentrations  the  mechanism  should  cease  to  function.  An  example 
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of  this  case  la  Ge.-In^Aa,  where  In  baa  provan  to  be  effective  In  dislocation 
reduction  of  CaAa  (8). 

We  have  deaonatrated  that  our  Modification  of  the  Harrlaon  bonding  theory 
accruately  predicta  the  observed  change  of  bond  lengtha  In  semiconductor  alloya 
and  offera  guidance  to  means  for  dlalocatlon  reduction. 
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ABSTRACT 

Several  recent  theoretical  studies  of  the  local  structure  of  semiconductor  alloys  are 
summarized.  First,  dilute  limit  calculations  of  local  bond  lengths  and  mixing  enthalpies 
are  discussed.  These  calculations  include  effects  due  to  both  bond  length  and  bond- 
angle  distortions,  as  well  as  local  chemical  rearrangements.  Then,  a  new  statistical 
theory  of  concentrated  alloys  is  described.  Deviations  from  random  alloy  distributions 
(microcluslers)  are  predicted. 


INTRODUCTION 

This  paper  summarizes  our  recent  theoretical  studies  directed  toward  understand¬ 
ing  the  microscopic  structures  of  pseudo-binary  semiconductor  alloys  AjBj.jC.  We  first 
present  a  detailed  calculation  of  the  local  bond  length  relaxation  in  the  dilute  limit 
x  —  0,  i.e.  the  case  where  an  A  atom  is  substituted  for  a  B  atom  in  a  BC  compound. 
The  mixing  enthalpy  parameter  0  is  found  to  be  related  to  small  excess  substitution 
energies.  These  excess  energies  are  calculated  directly  through  a  minimization  pro¬ 
cedure.  Thus,  the  accuracy  of  the  predicted  0  is  not  limited  by  trying  to  find  small 
differences  between  large  numbers.  The  theory  is  then  generalized  to  concentrated 
alloys  using  statistics  based  on  combinations  of  tetrahedral  clusters  of  five  atoms.  Our 
results  predict  that  tnicroclustering  occurs  in  a  majority  of  alloys.  We  conclude  by  iden¬ 
tifying  systematic  correlations  between  the  theory  and  several  experiments. 

Before  discussing  the  calculation,  it  is  useful  to  provide  some  background  about 
the  structure  of  these  alloys.  It  was  customary  to  assume  that  these  alloys  have  two 
sublattices  in  which  the  C  atoms  occupy  one  sublallice.  and  A  and  B  atoms  are  ran¬ 
domly  distributed  on  the  other.  This  picture,  referred  to  as  the  virtual  crystal  approxi¬ 
mation  (VCA).  implies  that  the  nearest- neighbor  (nn)  bond  lengths  in  the  alloy  are  the 
concentration  weighed  average  values,  i.e.  dAr  =  dg<~  =  d  =  x  d$l  +  (1  x)  d^c'  where 
the  values  with  a  superscript  |0|  denote  the  pure-crystal  values.  On  the  other  hand, 
according  to  Pauling's  covalent  radii  approximation  (CRA),  the  local  bond  lengths 
retain  their  respective  pure-crystal  values,  i.e.  dA<~  =  d$l  and  dgt-  =  d^j-1 

If  we  define  <*0  =  (d  d $ '  /  d  and  =  Id  dAC)  /  d.  then  the  ratio  t>/in  in  VC \ 
is  /ero.  hut  in  CRA  it  is  I.  However.  Mihtol-.-n  and  Boyce’11  found  from  their  FXaFs 
experiment  on  (iaxIni_xA5  that  the  nn  bond  lengths  do  not  fit  either  VCA  or  CRA. 
Instead,  they  found  the  value  of  b! f„  to  he  close  to  3/4.  Since  then,  similar  experiments 
hate  been  done  for  a  number  of  zinc-blende  pseudo-binary  alloys.'11  and  the  3/1  rule 
appears  to  be  quite  general. 


Mu  til  Soc  $v*ti©  P»oe  Voi  M  'MS  Muitnn  •Ifwi'c"  Sooity 


138 


DILUTE  LIMIT 

The  dilute  limit  is  the  easiest  case  but  is  still  Dot  trivial.  Its  solution  provides 
both  end-point  results  (x  =  0  and  1),  as  well  as  insight  into  the  extension  to  the  concen¬ 
trated  alloy  case.  A  complete  description  of  this  case  is  being  reported  elsewhere;  here 
we  summarize  the  essential  results.  The  substitution  energy  A,  for  an  A  atom  replacing 
a  B  atom  in  a  BC  compound  is  calculated  and  minimized  to  find  the  relaxed 
configuration.  A,  can  be  written  as  A,  =  d  (AEb  +  AE,  +  AEck),  where  AEb  is  the 
binding  energy  difference  between  the  AC  and  BC  compounds,  AE,  is  the  strain  energy, 
and  AEch  is  a  chemical  energy  shift.  All  AE'  s  are  energies  per  bond.  Then, 
AE  =  AE,  +  AEcb  is  the  excess  energy  per  bond  for  the  impurity  substitution.  AE,  is 
calculated  by  dividing  the  crystal  into  two  regions.  Outside  R  (which  is  the  distance  of 
the  second-shell  atoms  to  the  impurity),  the  distorted  crystal  is  treated  as  an  elastic  con¬ 
tinuum  with  a  radial  displacement  field  which  is  inversely  proportional  to  the  square  of 
the  radial  distance,  so  AE,*0*11’  =  l/4RCuJ,  where  C  is  an  effective  shear  coefficient, 

C  =  *  (1.6  (C„  -  C,2)  +  4.8  C«,),  (1) 

and  u  is  the  magnitude  of  the  displacement  at  R.  Inside  R,  the  strain  energy  A,*“*’  is 
treated  with  a  valence  force  field  (VFF).***  Finally,  the  chemical  energy  shift  AEth  is  cal¬ 
culated  from  Harrison's  model  and  arises  from  changes  in  the  metalizalion  energies 
caused  by  different  bond  lengths  Ad  ss  dBC  -  dAc  and  covalent  energies 
AV3  ==  Vj(AC)  -  Vj(BC).  Note  that  -  d$)  /  d$  and 

4  =  (d$  -  dAC)  /  d$  in  this  dilute  limit,  so  the  excess  energy  AE  can  be  expanded  up 
to  second  order  in  4,  u.  and  AV3  For  a  given  pair  A  and  B,  AE  is  an  explicit  function  of 
4  and  u.  Minimization  of  AE  with  respect  to  4  and  U  leads  to  the  equilibrium  local 
bond  length  dAC  and  energy  AE.  Then,  AE  is  used  to  estimate  the  mixing  enthalpy 
parameter  fl  in  the  mixing  eDthalpv  AH  =  x(l-x)fl  by 

n  =  2  (AE  (A  in  BC)  +  AE  (B  in  AC)).  (2) 


A  systematic  comparison  with  other  models  based  on  strain  energy  alone  shows 
that  an  increase  of  the  range  of  the  fixed  boundary  R  increases  the  relaxation  of  dAc,  i.e. 
it  causes  6/Sq  to  increase.  The  inclusion  of  the  bond  angle  restoring  force,  on  the  other 
hand,  reduces  the  relaxation.  It  turns  out  that  a  delicate  cancellation  of  these  two 
effects  causes  a  simple  spring  model  pointed  out  by  Shih  et  al.  (SSHS)*3’  to  yield  accu¬ 
rate  results.  In  this  model  4  =  40  /  1 1  +  1/3o/Q|),  w  here  a  and  Q|  are  the  bond¬ 
stretching  force  constants  for  the  host  (BC)  and  the  impurity  (AC)  crystals.  With 
a  =:=  <>|.  this  model  predicts  4  =  3/34^  for  a  zinc-blende  alloy.  Although  our  full  pertur¬ 
bation  theory  (FPT)  and  the  \TF  model  of  Martin  and  Zunger  (MZ)*4’  predict  dAC  with 
an  average  absolute  deviation  comparable  to  the  experimental  uncertainty  of  0.01A.  the 
simple  spring  model  is  e\en  better. 

We  note  that  while  in  our  theory.  MZ  and  -SSHS.  the  II  values  are  directly  calcu¬ 
lated  without  any  adjustable  parameters,  our  theory  and  JsSHS  agree  with  the  experi¬ 
ment  as  well  or  even  slightly  better  than  the  one-parameter  theories.1561  Althoug'  our 
theory  predicts  a  negative  fi  value  for  all  three  (Ga.  Al)  alloys,  the  magnitude 
(fl  =  -0.17  kcal/mole)  is  too  small  to  account  for  the  ordering  of  Ga,.\l,  ,As  grown  at 
800  to  700 'C  found  recently.*7’  The  calculated  0  values  also  provide  guidance  in 
separating  miscible  from  immiscible  alloys.  In  a  random  alloy,  the  criterion  for  alloy 
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mixing  for  all  x  is  T  >  Tc,  where  the  critical  temperature  T  is  given  by  Tt  =  n/2Rg 
with  Rj  being  the  universal  gas  constant.  Figure  1  is  a  plot  of  Tc/T2  against 
I  I  /  I  I  >  w*»«re  Tj  is  the  lower  of  the  two  constituent's  melting  temperatures,  and 
6m  =  1.63xm  with  Xm  being  the  ratio  of  rms  bond  length  amplitude  fluctuation  to  the 
bond  length  at  Tj.  The  simple  spring  model  gives  Tc/T2  =  (A)/^m)2  88  indicated  by  the 
solid  curves. 

There  is  an  empirical  rule  stating  that  a  miscibility  gap  will  occur  if  [  ^  | 
between  two  alloy  components  exceeds  7.5%.  However,  if  Tc/T2  is  plotted  against 
|  50  |  ,  the  simple  spring  model  would  not  exhibit  a  smooth  simple  quadratic  curve,  and 
our  theoretical  points  would  be  much  more  scattered.  This  suggests  that 
|  |  /  |  |  >  1  is  a  better  criterion  than  |  6$  |  >  0.075.  Figure  2  also  clearly  shows 
the  chemical  effects,  namely  negative  and  positive  chemical  energies  AEch  for  cation  and 
anion  substitutions  respectively.  The  full  theory  and  the  experiments  correlate  within 
the  experimental  uncertainties.  The  simple  SSHS  model  clearly  is  an  excellent  universal 
representation.  However,  Tc/T2  varies  faster  than  quadratically  for  larger  |  i^/6a  | 
values,  as  born  out  from  both  the  experimental  data  and  the  full  theory. 


CONCENTRATED  ALLOYS 

Turn  now  to  the  concentrated  alloy  case.  First,  an  improved  statistical  model  is 
required.  We  have  extended  regular  solution  theory  based  on  pair  energies  to  one  for 
five-atom  clusters.  For  an  A,Bj.,C  alloy,  the  building  blocks  are  clusters  of 
A(m)D(-l-m)C,  where  m  ranges  from  0  to  4.  For  a  given  alloy  concentration  x  and  for  a 
given  set  of  energies  associated  with  these  clusters,  we  have  derived  expressions  for 
the  cluster  population  distribution  xm  S  nm  /  N,  where  N  is  the  total  number  of  unit 
cells  and  nm  is  the  averaged  number  of  cells  with  A(m)B(4-m)C  clusters.  The  partition 
function  Z  is  obtained  using  a  steepest  descents  argument  which  then  yields  the  mixing 
Helmholz  free  energy  AF.  The  result  reduces  to  Guggenheim's  tetrahedron  case17'  if 
pair  potentials  (for  the  second-neighbors)  are  assumed.  Another  major  difference  is  that 
we  only  need  to  solve  a  single  quartic  equation,  while  Guggenheim  needed  to  solve  fotir 
simultaneous  quartic  equations. 

The  key  to  the  problem,  however,  lies  in  the  calculation  of  the  energies  <m.  If  one 
assumes  that  the  size  of  the  tetrahedra  for  all  m-clusters  at  a  given  alloy  concentration 
takes  on  the  corresponding  VCA  values  but  allows  the  central  C  atom  to  relax,  then  the 
energies  as  functions  of  x  behave  like  those  shown  in  Figure  2(a).  There  are  at  least  two 
major  flaws  in  this  result.  First,  the  energies  a r»  too  large  and  would  correspond  to  f! 
values  many  times  the  experimental  values.  Seennd,  at  x  =  0.75,  0.5.  and  0.2.  these 
energies  imply  compound  formation  for  AsB,C ,.  A2B2C4.  and  AiBjC4,  respectively, 
which  is  opposite  to  the  known  tendency  for  spinoda!  decomposition  of  GaJn^jAs  at 
low  T.  However,  if  the  local  cell  volume  of  eaeh  cluster  is  allowed  to  be  in  mechanical 
equilibrium  with  a  continuous  medium  with  an  effective  shear  coefficient 
('  =  \t\?-  •*  !1  where  the  C  value  for  'he  pure  material  is  given  by  Eq.  (1). 

'hen.  til-  corresponding  energies  tm  as  a  function  of  \  are  given  in  Figure  2(b).  which 
now  yields  a  reasonable  value  of  mixiric  enthalpy  and  correctly  predicts  the  tendency 
’  ward  spinoda)  decomposition  at  low  temperature.  With  this  set  of  energies,  one  can 
•heii  calculate  the  cluster  distribution  and  compare  them  with  the  corresponding 
value'  for  a  random  alloy,  i.e.  xi01  =  (4m  tx^l  1  -x  l4* m.  Figure  3  shows  the  deviation 
from  randomness  Axm  =  xm  -  x^0  as  a  function  of  x  for  four  arbitrarily  chosen  growth 
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FIGURE  1  PLOT  OF  Tc/T2  AS  A  FUNCTION  i60  <6m  FOR  Tc  OBTAINED  FROM  III  PRESENT  THEORY 
ANO  Ibl  THE  EXPERIMENTAL  0  IN  111.  THE  DOTS  ARE  FOR  ANION  ALLOYS  AND  CROSSES 
FOR  CATION  SUBSTITUTION  THE  SOLID  LINES  IN  BOTH  PARTS  CORRESPOND  TO  THE  0 
BASED  ON  THE  SIMPLE  THEORY  OF  REF  3  THE  DASHED  LINES  AT  Tc<  Tj  ■  0  SEPARATE 
THE  MISCIBLE  FROM  IMMISCIBLE  GROUPS 
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The  dislocation  energies  and  hardness  of  semiconductors  are  calculated  by  an  extension  of 
Harrison’s  method.  It  is  demonstrated  in  agreement  with  experiment  that  dislocation  energies  per 
unit  length  are  proportional  tod  d  ~9,  where  d  is  the  bond  length  and  hardness  is  proportional 
tod  ~  d  ~ 1 The  hardness  is  related  to  the  interaction  energies  among  dislocations.  It  is  argued 
that  dislocation  densities  of  semiconductors  will  be  reduced  if  they  are  alloyed  with  a  second 
constituent  that  has  a  shorter  bond  length.  Experimental  evidence  supporting  this  strategy  is 
noted. 

Dislocations  in  semiconductors  are  detrimental  to  de¬ 
vice  function;  they  serve  as  channels  for  impurity  migration 
and  trapping,  which  cause  nonuniform  doping  and  degrades 
p-n  junctions. 1  They  also  decrease  the  material’s  resistance 
to  plastic  deformation.  The  aim  of  this  letter  is  to  provide 
insights  into  the  underlying  physical  mechanisms  control¬ 
ling  dislocations  and  semiconductor  hardness,  and  then  to 
suggest  strategies  for  decreasing  dislocation  densities.  It  is 
well  established  that  the  hardness  of  tetrahedrally  coordi¬ 
nated  semiconductor  materials — groups  IV,  III-V,  and  II- 
VI  compounds — exhibits  a  sharp  variation  with  their  near¬ 
neighbor  distance d,  approximately  proportional  tod  9  for 
one  group  of  seven  compounds.2  Thus,  semiconductors  with 
small  lattice  constants  tend  to  be  harder  materials.  These 
same  materials  have  larger  stiffness  coefficients3  and  have 
fewer  dislocations  in  as-grown  crystals. 

The  shear  coefficients  (combinations  of  C44  and 
CM  -  C,2  in  the  Schoenflies  notation)7  depend  on  crystal 
orientation  and  (in  Harrison's  notation")  are  proportional  to 
y'2/d'(V]  +  y\)V2,  where  U2  ad  2  is  the  covalent  and  F, 
is  the  ionic  energy.  The  metallic  interaction  modifies  the 
functional  dependence  of  the  shear  coefficient  on  V2  and 
but  introduces  no  explicit  dependence  on  the  hopping  inte¬ 
grals,  denoted  U,  by  Harrison. 1 "  In  a  pure  covalent  material, 
the  bond  energy  is  proportional  to  V2  (or  d  2),  and  the  bond 
volume  is  ad',  hence,  in  this  case,  the  shear  coefficient  var¬ 
ies  as d  \  In  the  limit,  F,>  V2,  C, ,  —  Cu  ad  ".  For  most 
polar  semiconductors,  d  9  is  a  good  approximation. 

Hardness  is  determined  by  applying  a  known  force  F  to 
a  probe  of  a  prescribed  shape  driving  it  into  the  surface  of  the 
sample. 7  The  area  A  of  the  resulting  indentation  is  measured, 
and  the  hardness  is  the  force  per  unit  indented  area.  Many 
dislocations  must  be  formed  to  allow  the  probe  to  indent  the 
semiconductor  If  the  indenter  is  a  rectangular  pyramid, 
then  the  hardness  is  H  =  F /A  =  Fh /Ah  =  (T/Ah,  where  eT 
is  the  work  required  to  cause  the  indenter  to  penetrate  to  a 
depth  h.  A  side  view  of  the  indentation  in  a  cut  through  its 
center  is  illustrated  schematically  in  Fig.  1 .  The  top  of  the 
indentation  has  side  length  W,  thus,  A  =  W7.  The  Burger's 
vector  has  magnitude  b,  proportional  to  the  bond  length  d. 

The  number  of  dislocations  Nh  required  to  accommodate  an 
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indentation  to  depth  h  is  Nh  =  h  /b  —  j  If  cos  d  /b,  where  o 
is  the  angle  between  the  normal  to  the  tip  of  the  indenter  and 
aside.  Figure  1  also  shows  a  model  of  one  possible  configura¬ 
tion  of  the  dislocations.  The  edges  of  the  extra  atom  planes 
that  are  driven  from  the  indented  volume  into  the  bulk  of  the 
semiconductor  are  shown  as  lines  terminated  by  dots.  The 
dotted  ends  of  these  lines  are  the  positions  of  the  disloca¬ 
tions,  which  are  perpendicular  to  the  plane  of  the  figure.  The 
planes  driven  to  the  sides  each  have  a  finite  extent  and  a 
trapezoidal  shape.  The  planes  driven  down  under  the  indent¬ 
er  have  a  square  shape. 

Much  of  the  work  done  on  the  indenter  goes  into  the 
energy  to  form  the  indicated  dislocation  configuration,  al¬ 
though  some  certainly  goes  into  heat.  There  are  two  major 
contributions  to  this  formation  energy.  The  first  is  the  enegy 
needed  to  generate  each  dislocation  as  an  isolated  entity,  and 
the  second  is  the  interaction  energy  among  these  disloca¬ 
tions.  Because  the  interaction  term  dominates  H,  approxi¬ 
mations  made  to  simplify  the  first  term  are  relatively  unim¬ 
portant.  The  extra  planes  driven  to  the  sides  of  the 
indentation  have  a  finite  extent;  accordingly  (in  this  idealized 
picture),  there  are  both  edge  dislocations  at  their  base  and 
screw  dislocations  associated  with  their  termination.  The 
square  planes  driven  below  the  indentation  have  edge  dislo¬ 
cations  around  the  sides  and  screw  dislocations  at  the 
comers  to  make  the  turns.  Moreover,  there  are  interactions 


FIG  I  Schematic  representation  of  an  idealized  minimum-energy  disloca¬ 
tion  configuration  produced  by  a  square  cross-section  indenter  in  a  hard¬ 
ness  measurement  IF is  the  side  length  and  h  is  the  depth  of  the  indentation 
In  this  ideal  case,  half  the  material  from  the  indented  region  is  displaced 
along  the  glide  planes  (indicated  by  the  dashed  arrows)  lo  the  sides  and  half 
is  displaced  below  the  indenter 


among  the  dislocations,  which  can  produce  a  minimum-en¬ 
ergy  configuration.  For  the  arrangement  depicted  in  Fig  1, 
it  always  costs  energy  to  position  a  second  dislocation  on  a 
parallel  glide  plane  to  one  already  present.  However,  the 
magnitude  of  this  extra  energy  can  be  minimized  and,  for 
proper  configurations  of  the  dislocations,  there  are  attractive 
forces  along  the  glide  planes  that  will  tend  to  position  the 
dislocations  into  the  minimum-energy  configuration.  The 
minimum  configuration  arises  when  half  of  the  atoms  from 
the  indented  volume  go  respectively  to  the  side  and  below  the 
indenter  Then,  in  both  regions,  the  maximum  angle  made 
between  successive  close-spaced  dislocation  lines  and  their 
glide  planes  is  d  =  n/2  —  d  /2,  as  shown  in  the  figure.  This  is 
the  minimum  realistic  energy  configuration.  If  the  disloca¬ 
tions  are  separated  more  than  shown  in  Fig.  1,  then  there  is 
more  volume  of  strained  material  and  the  interaction  energy 
would  be  larger  still. 

An  approximate  expression  for  the  energy  required  to 
indent  the  material  is7 


*r  =  2-4 
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In  the  first  term,  using  an  isotropic  medium  approximation 
and  neglecting  core  terms,  the  energy  per  unit  length  to  form 
an  edge  dislocation  is’ 
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the  shear  coefficient  is  G,  the  Burger’s  vector  b  -  d  / 3  for  an 
indentation  along  a  { 100)  axis,  the  range  of  the  elastic  defor¬ 
mation  of  a  dislocation  R  is  taken  equal  to  W  (for  want  of  a 
better  approximation),  r„~ d  is  the  dislocation  core  radius, 
the  Poisson  ratio  is  v  ~0  2  for  most  semiconductors,  and  L 
is  the  length  of  the  nli  dislocation  In  the  second  term,  £',(  is 
the  interaction  energy  per  unit  length  between  dislocations  / 
and  j.  Assuming  they  have  parallel  glide  planes  and  their 
Burger’s  vectors  have  the  same  sign,  Eu  is  given  by 
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where  r  is  the  separation,  and  <f4)  is  the  angle  that  a  line 
perpendicular  to  and  joining  i  he  dislocations  makes  w  ith  the 
glide  plane,  as  shown  in  f  ig  I  Because  the  various  disloca- 
tions  in  a  region  have  different  lengths,  the  net  interaction 
energy  iv  approximated  by  multiplying  the  energy  per  unit 
length  by  the  length  oi  the  shorter  one  I  he  upper  limits  on 
the  sums  \  arc  the  number  of  dislocations  in  one  region  (side 
or  bo! tom )  associated  with  one  edge  For  the  minimum-ener- 
gv  i  onliguration.  \  Y„ ,  2  The  four  that  multiplies  the 
bracket  accounts  for  the  four  sides,  and  'he  two  for  the  two 
regions  fo:  each  side 

We  now  eix-ouMtci  (nil  first  surprise  As  we  can  see 
from  Pq  ,2i,ut<f  >•.  h  and  E  have  comparable  magnitudes 
Because  thet •  .nr  approximately  N  terms  in  the  interaction 
energv  sum.  only  V  terms  in  l  he  formation  energy  sum.  and 
V  is  a  l.i  r  go  number  A  >  1  '•  the  interaction  energy  completely 
dominates  the  hardness  In  fact.  A  is  tvpically  of  the  order 
I  < 1  Hem  e,  terms  ow  ing  ti  >  sc  tew  dislocations,  core  energies, 
heat  dissipation  the  disli>calions  propagate  to  their  places, 
and  other  effects  associated  wrh  the  first  term  arc  unimpor¬ 


tant.  However,  care  must  be  taken  with  the  interaction 
terms:  Eq.  (1)  neglects  a  number  of  secondary  interactions, 
some  positive  and  others  negative;  these  will  be  added  later 
The  principal  neglected  terms  are  the  interactions  between 
the  dislocations  in  the  different  regions  on  each  side  (posi¬ 
tive),  and  the  interactions  between  adjacent  side  and  bottom 
regions  (positive),  and  the  interaction  between  the  opposite 
sides  in  the  bottom  regions  (negative).  Comparison  of  the 
results  with  experiment  will  indicate  how  important  these 
neglected  terms  are  likely  to  be.  The  length/.,  is  W  [X  -  i)/N 
for  i  from  1  to  .V  This  is  the  largest  length  of  the  side-inserted 
planes  and  its  choice  partially  accounts  for  interacitons 
between  the  otherwise  neglected  screw  dislocations.  The  dis 
tance  r;l  -  r  r;  is  given  by  2V2  b  \i  -  j\  for  /  and  j  ranging 
from  1  to  X  for  the  minimum-energy  configuration  and  a 
tetrahedrally  bonded  semiconductor.  Finally,  in  the  indical 
ed  configuration,  n/2  -  d  /2  Inserting  these  ex¬ 

pressions  into  Eq  1 1 1  and  retaining  only  terms  of  order  .V  ’ 
yields 
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One  can  also  get  a  number  for  the  hardness  of  a  disloca¬ 
tion  in  which  all  the  material  is  pushed  along  the  same  glide 
plane,  e  g  .  to  the  side,  to  the  bottom,  or  normal  to  the  face  of 
the  indenter  la  possibility  not  depicted  in  Fig  ll  In  this  case, 
the  factor  of  2  in  front  of  Eq  ill  is  removed,  X  and 
i i  n/2  -  d  Then  a  higher  nonequilibrium  hardness  in  t he 
context  of  this  model  (denoted  //,!,  is  obtained 
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The  proper  answer  for  most  materials,  and  depending  on 
crystal  orientation,  probably  lies  somewhere  between  H 
and  //,  For  an  indenter  with  d  n/4.  we  have  // . 


0  09696' /l  1  v|  and  //,///„ 


2  39  Harrison*  has 


show  n  that  one  contribution  to  the  shear  coefficient  tactually 
C, ,  C,.|  is  6  2.38  frma'/md \  where  m  is  the  free 

electron  mass,  is  the  covalence,  rr,  I  . -'(I';  ,  F  j  |'  , 
and  d  is  the  bond  length  We  will  approximate  6  bv  this 
expression  l  sing  this  6  and  v  0  2.  and  changing  the  di 
mensnms  to  those  in  terms  of  w  hich  experimental  hardness 
numbers  are  customarily  quoted  gives  2  38  *  104 

in'  ;d'  i  kgm/mm'.  w  here  d  is  in  angstroms  Calculated  v  a  I 
ues  T  //„„„  and  //,  are  ploited  against  experimental  results 
in  I  ig  2  for  a  number  of  semiconductors 

Figure  2  has  ihe  theoretical  //mm  and  H ,  values  eon 
nected  by  arrows  from  //m,„  to  //,  for  each  compound,  plot 
ted  as  a  function  of  Ihe  corresponding  experimental  val 
ues  ”  If  the  theory  were  perfect  and  ihe  experimental  values 
were  accurate,  the  points  would  fall  on  the  indicated  unity 
slope  line.  Several  conclusions  can  be  drawn  Firstly,  the 
order  of  magnitude  of  the  predicted  and  measured  values  are 
the  same,  a  result  obtained  with  no  adjustable  parameters  in 
ihe  theory  Secondly,  ue  trends  from  one  compound  to  an 
other  are  properly  given  by  the  theory  Although  the  //„„„ 
values  are  generally  too  small,  they  fit  the  soft  materials  bet¬ 
ter,  and  the  //,  values  fit  the  harder  materials  better.  Thirdly, 
from  Eqs.  (1|  and  (3),  H  is  given  in  a  rough  but  revealing 
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Mi»  2  I  hcoretical  vs  experimental  hardness  of  several  semiconductors  in 
Refs  2  and  9  I  he  two  theoretical  values  for  each  semiconductor  are  con- 
rut  ted  h>  an  arrow  from  t« »  H  Perfect  agreement  would  correspond 
|H»irtts  being  on  the  units  slope  solid  lin** 

approximation  by  4.V  A  /  /Ah,  where  the  averge  disloca¬ 
tion  length  /  H  / 2  anil  F  is  the  average  dislocation  pair 
interaction  energy  per  unit  length  Notice  that  .V  lor  W I  can¬ 
cels  from  this  expression,  thus,  //  is  independent  of  B7  lor  F I 
and  therefore.  //  is  truly  a  measure  of  the  properties  of  the 
material  This  result  would  not  be  found  if  the  dislocation 
energies  I  >  .V  I  were  to  dominate  //  rather  than  the  pair  inter¬ 
action  energy  I  ■  .V  i  Finally,  the  Berger's  vector  cancels 
front  the  leading  term  and  appears  only  in  the  argument  of 
the  logarithm  in  F.q  1 4 1  I  hus.  the  answers  are  also  insensi¬ 
tive  lo  Ms  choice 

Dislocations  are  often  found  in  materials  as  they  are 
grown  Their  density  is  determined  by  the  thermal  and  me¬ 
chanical  stresses  to  which  they  are  subjected  in  the  growth 
process  A  dislocation  constitutes  a  metastable  excitation 
relative  to  the  perfect  crystal  ground  state  At  the  elevated 
growth  temperatures  and  temperature  gradient  behind  the 
growth  front,  the  number  of  dislocations  present  is  con¬ 
trolled  by  the  relative  rate  at  which  vacancies  anneal  or  con¬ 
dense  into  dislocations  "  The  dislocation  formation  rate  will 
be  slower  in  a  material  grown  at  the  same  temperature  if  E,  is 
higher  If  an  alloy  is  formed  from  a  materia!  of  interest  and  a 
second  constituent  with  a  shorter  bond  length,  one  expects 
the  average  bond  energy  land  thus  both  the  melt  temperature 
and  vacancy  formation  energy!  to  increase  proportional  to  a 
low  inverse  power  of  the  average  bond  length  "  Hence,  the 
equilibrium  vacancy  density  |ust  below  the  melting  point 
tends  to  be  the  same  in  lowest  order  for  all  materials,  inde¬ 
pendent  of  the  bond  lengths  of  the  constituents.  Flowever, 
the  shear  coefficient  and  dislocation  energy  per  unit  length 
will  increase  with  much  higher  inverse  powers  of  the  bond 
length  Consequently,  dislocation  densities  should  be  re¬ 
duced  in  such  alloys  relative  to  those  found  in  the  longer 
bond  length  pure  constituent  This  expectation  is  confirmed 


in  the  recent  work  reported  on  Zn,  ,Cd,Te  bulk  material.4 
The  best  CdTe  that  has  ever  been  grown  has  dislocation  den¬ 
sities  in  excess  of  5  x  10'  cm  \  The  addition  of  only  4%  Zn 
reduced  the  dislocation  count  to  less  than  5  X  104  cm  2.  The 
ZnTe  bond  length  is  2.643  A,  while  the  CdTe  bond  length  is 
2.805  A,  a  6%  difference.  This  6%  difference  in  bond  length 
translates  into  a  2%  difference  in  the  dislocation  energy  per 
unit  length  for  1  -  x  =  0.04.  Dislocation  energies  per  unit 
length  are  typically  10  eV  per  lattice  spacing;  accordingly,  a 
2%  increase  can  be  expected  to  slow  their  formation  rate 
considerably. 

The  argument  just  presented  naturally  leads  to  a  strate¬ 
gy  for  decreasing  dislocations  in  other  semiconductors.  If  an 
alloy  is  made  of  the  material  of  interest  with  another  com¬ 
pound  with  a  shorter  bond  length,  then  the  dislocation  den¬ 
sity  should  be  reduced.  For  example,  this  suggests  that  the 
addition  of  a  small  amount  of  GaP(r/  =  2.359  A)  may  signifi¬ 
cantly  reduce  the  dislocation  density  of  bulk  grown  GaAs 
Id  =  2.448  A).  It  has  been  demonstrated  that  the  addition  of 
approximately  \%  GaN  id  =  1  946  A)'  *  or  of  a  10‘*  cm  ’ 
BAs|rf  =  2  069  A  (concentration"  to  GaAs  can  yield  a  large 
volume  of  dislocation-free  material.  An  InAs  additive  with 
its  longer  bond  length^  -  2  623  A  (serves  the  same  function 
indirectly,  by  causing  GaAs  bonds  in  its  neighborhood  to  be 
compressed.  This  indirect  mechanism  should  be  less  effec¬ 
tive  than  substituting  short  bond  length  additives 

We  have  demonstrated  that  the  dislocation  energies 
and  hardness  of  tetrahedrally  bonded  semiconductors  are 
rapid  functions  of  the  reciprocal  of  the  bond  length.  This 
rapid  d  dependence  of  dislocation  energies  provides  a  ratio¬ 
nale  for  the  dramatic  decrease  of  the  dislocation  density  in 
bulk  grown  Zn,)mCd(1,6Te  material  relative  to  that  found  in 
CdTe,  and  suggest  means  for  accomplishing  the  same  ends  in 
other  materials 
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The  sensitivity  of  defect  energy  levels  in  semiconductors  to  the  host  band  structures  and  impurity 
potentials  has  been  studied  for  approximately  30  impurities  in  CdTe  using  four  different  band- 
structure  models  The  discrepancies  in  the  defect  levels  between  two  different  sets  of  band  struc¬ 
tures  and  impurity  potentials  are  found  to  range  from  less  than  0.1  eV  to  the  whole  band  gap  (1.6 
eV).  The  band-structure  effects  are  analyzed  here  in  terms  of  detailed  partial  densities  of  states. 

Examples  of  contradictory  predictions  from  different  band  structures  are  illustrated,  and  ways  to 
improve  the  theory  are  suggested 


I  INTRODUCTION 

In  several  of  our  recent  papers,1”7  we  have  applied  a 
method  to  calculate  the  band  structure  of  semiconductors 
that  is  both  efficient  and  accurate.  Because  the  procedure 
involves  casting  the  basis  functions  into  orthonormal  local 
orbitals6  (OLO),  our  method  has  the  advantages  common 
to  empirical  tight-binding  (ETB)  calculations,8 “ 10  except 
that  the  Hamiltonian  matrix  elements  to  all  ranges  are  re¬ 
tained.  The  inclusion  of  these  higher  coefficients  makes  it 
possible  to  produce  excellent  band  structures  including 
conduction  bands  and  effective  masses.  The  method  also 
yields  wave  functions  for  optical  property  calculations.7 
Moreover,  its  OLO  description  also  permits  its  extension, 
through  the  coherent-potential  approximation,  to  al¬ 
loys.2  " ' 

The  recent  attention  focused  on  defects  in  semiconduc¬ 
tors  motivated  us  to  apply  our  method  to  this  problem. 
The  theories  of  defects  have  ranged  from  very  sophisticat¬ 
ed  self-consistent  density-functional  theory11”13  (SCDFI 
to  simple  ETB  calculations.  It  is  generally  recognized 
that  SCDF  is  as  accurate  in  defects  for  the  ground-state 
properties  as  it  is  for  pure  semiconductors,  but  less  certain 
in  assigning  excited  energy  levels.  ETB,  because  it  can 
produce  results  for  many  systems  in  one  study,  claims  to 
predict  the  trends  of  deep  levels’"  even  if  the  accuracy  for 
a  given  impurity  may  be  poor.  However,  this  contention 
remains  to  be  verified. 

To  assess  this  concern,  we  ask  the  following  question 
"How  sensitive  arc  defect  levels  to  host  band  structures 
and  impurity  potentials’"  To  this  end,  we  have  adopted 
the  simple  yet  nontrivial  defect  model,  that  of  site- 
diagonal  substitutional  defects  often  used  in  ETB  studies. 
CdTe  was  selected  in  this  study  because  its  hand  structure 
has  been  examined  in  great  detail  by  us,  and  there  are 
three  published  band-structure  models*  "  that  we  could 
easily  generate  for  comparison.  There  is  also  a  consider¬ 
able  body  of  expenmental  data  on  deep  states  in  this  svs- 
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II.  CALCULATION AL  PROCEDURE 

In  the  simple  site-diagonal  substitutional  defect  model, 
the  impurity  energy  levels  E  are  determined  by  the  equa¬ 
tion 

l-eag,(£)  =  0  ,  (1) 

where  a  designates  the  symmetry  of  a  local  state,  e.g.,  T6, 
T 7,  and  r8  on  an  atomic  site  in  the  zinc-blende  structure, 
and  ga  is  the  real  part  of  the  diagonal  matrix  element  of 
the  host-crystal  Green  function.  ga  can  be  calculated 
from  the  partial  density  of  states  (PDOS)  by 

gn(£)  =  J pa(e)/IE -c.)de  .  (2) 

The  PDOS  is  given  by 

pa(€  1  =  2  I  a“(*>  i  .  (3) 

n,k 

where  en(k)  are  band  energies  and  a%(k)  are  the  probabili¬ 
ty  amplitudes  of  the  band  state  in  the  Bloch  basis  con¬ 
structed  from  the  OLO  labeled  by  a.  The  Brillouin-zone 
integration  in  Eq.  (3)  is  calculated  using  an  accurate  ray 
scheme. 18 

Because  a  principal  concern  of  this  paper  is  the  sensi¬ 
tivity  of  impurity  levels  to  the  host  band  structures,  we 
should  emphasize  the  difference  between  our  method  and 
ET  B.  Our  method  consists  of  four  steps. 

(1)  We  start  with  four  Gaussian  orbitals  per  atom  and 
empirical  pseudopotentials,10  and  compute  the  Hamiltoni¬ 
an  matrix  Htk)  and  overlap  matrix  Sik)  as  was  done  by 
Kane2"  and  Chadi.21 

(2)  The  Gaussian  orbitals  are  transformed  into  OLO,6 
so  H  (k)  is  transformed  into  //„(  k  )  and  S  into  the  identity- 
matrix.  The  band  structures  calculated  from  H„(k  t  are 
accurate  to  5 Cr  as  compared  to  more  sophisticated 
methods  using  the  same  potential.1 

(3)  A  spin-orbit  Hamiltonian  in  the  OLO  basis4  is  in¬ 
corporated  to  deal  with  this  interaction. 

(4)  To  compensate  for  the  effects  of  truncated  basis  and 
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nonlocal  potentials,  a  perturbation  Hamiltonian  Hx  is 
added.  Ht  has  the  same  form  as  a  truncated  ETB  Hamil¬ 
tonian.  The  parameters  in  Hx  are  adjusted  to  fine  tune 
the  important  band  energies  and  effective  masses.1-4 

Although  both  ETB  and  our  methods  are  empirical, 
there  are  two  major  differences. 

(1)  While  most  ETB  retains  the  H  matrix  elements  only 
to  the  first-  or  second-neighbor  shell,  ours  extends  to  all 
ranges,  so  that  the  high  Fourier  components  needed  to 
produce  the  sharp  band  curvatures  are  properly  given. 

(2)  Our  method  can  directly  generate  wave  functions  for 
calculation  of  other  properties. 

Thus,  while  our  method  yields  more  accurate  band  struc¬ 
tures,  it  retains  much  of  the  advantage  of  ETB,  namely 
the  computational  speed  and  a  simple  direct-space 
description  of  the  Hamiltonian. 


III.  BAND  STRUCTURES  AND  PARTIAL 
DENSITIES  OF  STATES 

Figure  1  depicts  the  four  band  structures  to  be  con¬ 
sidered  for  CdTe.  Our  result  is  in  panel  (a);  panels  (b) 
(Ref.  8)  and  (c)  (Ref.  9)  are  two  ETB  band  structures  with 
the  Hamiltonian  matrix  elements  truncated  at  second 
neighbors.  (Because  different  parameters  were  selected, 
these  two  band  structures  are  not  identical.)  Panel  (d) 
(Ref.  10*  results  from  the  use  of  five  basis  orbitals  per 
atom;  the  extra  one  is  an  excited  s  state.  All  these  band 
structures  are  adjusted  to  have  the  proper  fundamental 
band  gap  of  1.6  eV.  The  principal  differences  one  sees  on 
first  inspection  are  in  the  band  curvatures,  especially  the 
conduction  bands.  The  effective  mass  at  the  bottom  of 
the  conduction  band  in  panel  (a)  is  0.1  times  the  free- 
electron  mass,  in  agreement  with  experiment,17  while  in 
other  panels  it  is  more  than  twice  as  large. 

Figure  2  shows  the  densities  of  states  (DOS)  for  each  of 
the  band  structures  in  Fig.  1.  While  the  valence  bands  at 
least  exhibit  general  common  features,  the  conduction 
bands  are  almost  unrecognizable  as  representing  the  same 
compound.  In  panels  (c)  and  (d),  for  example,  there  is  a 
second  band  gap  above  the  fundamental  gap.  Also  note 
that  there  are  two  extra  narrow  peaks  associated  with  the 
two  extra  excited  s  orbitals  (one  for  Cd  and  the  other  for 
Te)  included  in  the  calculation. 

To  analyze  the  band  effects  on  defect  levels  [see  Eqs.  (1) 
and  (2)],  the  DOS  is  further  decomposed  into  partial  den¬ 
sities  of  states  for  rb( s >,  T 7(pl/2),  and  T 8(pJ''2)  states  on 
the  Cd  and  Te  sites,  as  shown  in  Figs.  3—6.  The  T 8 
PDOS  are  not  shown  because  they  are  nearly  the  same  as 
r7  with  only  a  slight  upward  energy  shift.  These  PDOS 
show  how  the  "atomic"  levels  evolve  into  band  states. 
These  curves  contain  useful  information  about  many 
properties,  e  g.,  the  relation  between  the  crystal  bonding 
and  atomic  energies,  and  how  potential  disorder  in  alloys 
affects  different  parts  of  the  bands,2  “ 5  in  addition  to  de¬ 
fect  levels  studied  here. 

The  T6(Cd>  PDOS  shown  in  Fig.  3  split  between  the 
conduction  and  valence  bands.  It  is  generally  assumed 
that  the  cation  s  states  in  III- V  and  II- VI  compounds 


r  x  w  l  r  k.u  x  r  x  w  l  r  k.u  x 

k  (2ir/a) 


FIG.  1.  Four  band  structures  of  CdTe  used  for  comparative 
studies:  (a)  present  work,  (b)  Ref.  8,  (c)  Ref.  9,  and  (d)  Ref.  10. 


FIG.  2.  Densities  of  states  calculated  from  the  four  band 
structures  in  Fig  I. 
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evolve  into  the  conduction  bands,  while  the  anion  p  states 
make  up  most  of  the  major  valence  bands  just  below  the 
gap.  Thus  it  is  perhaps  a  surprise  to  see  a  prominent  peak 
derived  from  the  cation  j  states  at  the  bottom  of  the  ma¬ 
jor  valence-band  structure.  However,  this  is  a  general 
feature  for  all  sp’-based  compound  semiconductors. 
These  are  the  states  responsible  for  the  first  observed 
breakdown  of  the  virtual-crystal  approximation  for  a 
semiconductor  alloy:  Hgi^Cd^Te  (which  is  caused  by 
the  large  s-energy  shift  between  the  Cd  and  Hg  sites). 4,5,22 

A  more  detailed  examination  draws  attention  to  some 
important  differences  among  the  four  panels  in  Fig.  3: 
the  valence-band  peak  in  panel  (c)  is  about  2  eV  higher 
than  the  rest,  and  it  is  also  high  compared  to  experi¬ 
ment.22  Our  conduction-band  PDOS  in  panel  (a)  is 
broader  than  the  others.  The  ratio  of  the  integrated 
PDOS  in  the  conduction  bands  to  that  in  the  valence 
bands  in  our  model  is  larger  than  those  in  other  panels. 
Also  our  PDOS  just  below  the  valence-band  edge  is  obvi¬ 
ously  smaller  than  that  found  in  other  models. 

Figure  4  shows  that  the  Cd  p  states  are  concentrated  in 
the  conduction-band  states.  This  is  particularly  true  in 
panel  (a),  where  their  contribution  to  the  valence-band 
states  shrinks  almost  to  nothing.  In  other  panels,  there 
are  still  sizable  (~20%)  valence-band  states.  In  contrast, 
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all  four  panels  in  Fig.  5  show  that  the  Te  s  states  are  con¬ 
fined  to  the  deep  valence-band  states,  as  generally  recog¬ 
nized.  Finally,  Fig.  6  shows  that  the  Te  p  states  dominate 
the  upper  valence-band  states.  Panel  (a)  has  much  less 
conduction-band  content  than  the  other  three  panels.  As 
we  will  see,  these  differences  can  result  in  quantitatively 
or  even  qualitatively  different  predictions  about  the  deep 
levels. 

IV.  IMPURITY-LEVEL  DETERMINATION 

A  convenient  way  to  study  the  impurity  energy  levels 
using  Eq.  (1)  is  to  rewrite  it  as  ua  =  l/ga(£)  and  plot  E  as 
a  function  of  v.  Once  this  E-v  curve  is  deduced  for  each 
a,  the  deep  levels  Ea  for  a  given  impurity  can  be  read  off 
the  curve  by  drawing  a  vertical  line  at  the  appropriate 
value  of  va  for  the  impurity.  We  set  the  zero  of  energy  at 
the  top  of  the  valence  bands.  Because  the  gap  is  1.6  eV, 
we  will  focus  on  levels  in  the  energy  range  from  —0.5  to 
2.0  eV. 

Calculations  have  been  performed  for  all  neutral  impur¬ 
ities  listed  in  Table  I.  Because  we  do  not  believe  that 
there  exists  a  uniformly  accepted  table  for  v  we  have 
adopted  a  table  that  we  used  for  structural  studies.23,24 
Table  I  lists  the  term  values,  which  we  obtained  from  to- 


TABLE  I.  s-  and  p-state  correlated  term  values  in  units  of  —  eV.  The  top  entry  is  the  s-state,  the 
second  the  p,/2- state,  and  the  third  the  p)/2- state  energy.  (All  energies  are  negafiue.) 
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V  leV) 

FIG.  7.  The  E-v  curves  for  the  1%  states  on  a  Cd  site. 

tal  energy  differences  between  atomic  configurations  cal¬ 
culated  using  the  norm-conserved  pseudopotentials25  and 
self-consistent  charge-density-functional  theory,  with  the 
first  ionization  energies  adjusted  to  be  the  experimental 
values.26  These  term  values  are  found  to  yield  consistent¬ 
ly  better  structural  properties23  in  Harrison’s  theory27  28 
than  those  based  on  Mann’s  values29  adopted  by  Har¬ 
rison.2*  The  impurity-potential  parameters  will  then  be 
taken  as  the  difference  of  the  term  values  between  the  im¬ 
purity  atom  and  Cd  (or  Tel.  To  study  the  sensitivity  of 
Ea  to  va,  we  shift  va  by  ±0.5  eV  and  compute  the  corre¬ 
sponding  changes  in  the  energy  levels. 

Figures  7—10  display  the  E-v  curves  for  several  a. 
Each  figure  has  four  curves,  corresponding  to  the  four 
panels  of  PDOS  in  each  of  Figs.  3—6.  The  functional 
behavior  of  these  curves  can  be  understood  qualitatively 
using  Eq.  (2)  and  Figs.  3—6.  If  E  lies  in  the  gap,  the  con¬ 
tribution  from  conduction  bands  is  negative,  but  positive 
from  the  valence  bands.  The  closer  the  PDOS  to  the  E  in 
question,  the  larger  will  be  its  influence.  Applying  this 
argument  to  the  r6(Cd)  representation,  we  see  that  the 
curves  in  Fig.  7  are  negative  in  the  gap  region  because  the 
PDOS  in  Fig.  3  near  the  bottom  of  the  conduction  bands 
are  much  larger  than  those  near  the  valence-band  top. 
Thus,  on  the  Cd  site,  only  impurities  with  an  s  energy 
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FIG.  8  The  E-r  curves  for  the  f  -  slates  on  a  Cd  site 
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FIG.  9.  The  E-v  curves  for  the  T*  states  on  a  Te  site. 

below  the  Cd  s  level  (  —  8.99  eV)  will  produce  a  T6  level  in 
the  gap.  However,  we  note  that  in  Fig.  7,  ga(£)  =  0  for 
models  (b)  and  (c)  just  below  the  valence-band  edge  be¬ 
cause  of  cancellation  between  the  conduction-  and 
valence-band  contributions.  At  this  E  value,  the  E-v 
curve  switches  from  u  =  —  oo  to  v=  <x>  (not  shown);  an 
ideal  vacancy  level  (corresponding  to  oa  =  oo )  is  located  at 
this  E.  A  similar  consideration,  but  with  the  conduction 
and  valence  bands  interchanged,  leads  to  an  understanding 
of  the  curves  in  Fig.  10.  Using  the  same  principle,  we  can 
easily  understand  why  all  curves  in  Fig.  9  for  the  T6(Te) 
representation  are  positive,  but  the  reasons  for  the  large 
displacements  between  these  curves  are  not  easy  to 
deduce.  In  Fig.  8,  the  curve  labeled  a  is  distinctively  dif¬ 
ferent  from  other  curves,  because  the  PDOS  in  panel  (a)  in 
Fig.  4  is  completely  dominated  by  the  conduction  band; 
however,  for  the  other  panels  the  PDOS  just  below  the 
valence-band  edge  are  as  large  as  those  just  above  the 
conduction-band  edge.  This  produces  a  very  sharp  nega¬ 
tive  E-v  curve  for  (a),  but  split  behavior  for  (b),  (c),  and 
(d). 

These  E-v  curves  provide  a  clear  picture  of  how  dif¬ 
ferent  host  band  structures  may  affect  the  deep  levels. 
Numerical  values  for  the  impurity  levels  can  be  obtained 
from  these  figures  by  drawing  vertical  lines  at  the  ap- 
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FIG  10  The  H  i  curves  for  the  T-  Males  on  a  Tc  sue. 
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propnate  impurity  potentials  (i.e.,  differences  between  the 
term  values  listed  in  Table  I),  as  has  been  shown  for 
several  representative  impurities.  To  provide  a  more 
quantitative  comparison,  Table  II  lists  some  calculated 
impurity  levels  Ea  and  the  corresponding  changes  A Ea 
due  to  the  1-eV  change  in  va. 

V  RESULTS  AND  CONCLUSION 

To  summarize  we  recall  that  band  models  (b)  and  (c) 
are  the  same  second- neighbor  ETB  with  two  different  sets 
of  parameters,  and  model  (d)  is  a  first-neighbor  ETB  with 
one  extra  s  orbital  per  atom.  Our  model  [model  (a)]  has 
the  form  of  ETB  but  is  derived  in  a  very  different  manner 
and  includes  all  the  long-range  interactions.  Therefore, 
we  expect  that  the  results  from  models  (b)  and  (c)  will  be 
close,  model  (d)  will  have  larger  discrepancies  from  (b) 
and  (c)  than  that  between  (b)  and  (c),  and  model  (a)  will 
differ  even  more.  This  is  evident  from  Figs.  7—10  and 
Table  II.  We  found  the  energies  for  the  T6(Cd),  T7(Te), 
and  rg(Te)  states  produced  by  models  (b)  and  (c)  agree 
within  0.1  eV.  For  the  other  states,  i.e.,  r6(Te),  r7(Cd), 
and  Tg(Cd),  the  energies  from  (b)  and  (c)  are  qualitatively 
similar,  but  the  difference  can  be  as  large  as  0.4  eV.  The 
largest  discrepancy  between  models  (d)  and  (b)  [or  (c)]  is 
more  than  0.5  eV,  and  that  between  (a)  and  other  models 
is  more  than  1  eV.  The  largest  difference  comes  from  the 


p  levels  on  a  Cd  site.  For  example,  the  filled  p  level  of  C 
on  a  Cd  site  in  model  (a)  is  a  resonance  state  just  below 
the  valence-band  edge  but  is  a  donor  state  in  the  other 
models.  Similarly,  model  (a)  puts  the  neutral  Te  antisite 
defect  p  levels  at  about  7  and  7  of  the  gap  [£(T7  )  =  0.48 
eV  and  Ell'll =0.95  eV),  while  other  models  assign  them 
as  resonance  states  inside  the  conduction  bands.  We  also 
note  that  the  discrepancies  between  different  models  are 
not  uniform,  but  vary  with  va.  Consider  r6(Cd)  for  ex¬ 
ample.  All  four  models  yield  the  same  ordering  and  about 
the  same  energies  for  the  group-III  impurities  Al,  In,  and 
Ga.  However,  as  u  becomes  more  negative,  the  splitting 
between  the  curves  increases,  so  the  discrepancies  become 
larger  [~  1  eV  difference  between  models  (a)  and  (d)  for  I 
impurity].  Similarly,  for  the  T7(Te)  states,  all  four  models 
put  the  Sn  impurity  energies  close  to  the  valence-band 
edge,  but  the  agreement  deteriorates  as  va  increases. 

Regarding  the  sensitivity  of  energy  levels  to  impurity 
potentials.  Table  II  shows  that  a  1-eV  shift  in  va  produces 
a  change  in  Ea  ranging  from  less  than  0.1  to  0.65  eV. 
Very  little  is  known  about  the  size  or  trends  in  errors  in¬ 
troduced  in  va  from  the  use  of  atomic  term  values.  How¬ 
ever,  we  know  that  the  discrepancy  of  va  between  two  dif¬ 
ferent  tables  of  atomic  term  values  can  be  larger  than  2 
eV.  This  discrepancy  translates  into  an  uncertainty  of  less 
than  0. 1  to  more  than  1  eV  in  the  impurity  energy  levels. 


TABLE  II.  Defect  energy  levels  E  and  changes  A£  due  to  a  1-eV  change  in  the  impurity-potential  parameter.  All  energies  are  in 
units  of  eV  Fn  stands  for  ideal  vacancy. 


Defect 

Model 

E 

(a) 

AE 

Model  (b) 
E 

AE 

Model  (c) 
E 

A  £ 

Model  (d) 

£ 

A  £ 

r6  on 

Cd  site 

Ga 

1.29 

0  39 

1.42 

0.24 

1.33 

0.23 

1.57 

0.18 

C 

0.21 

0.09 

0.38 

0.09 

0.36 

0.13 

0.74 

0.08 

Si 

0.67 

0.30 

1.02 

0.10 

0.93 

0.19 

1.27 

0.15 

P 

-0  19 

0.11 

0.39 

0.09 

0.38 

0.08 

0.75 

0.08 

O 

<  -0  s 

-0.02 

0.02 

0.04 

0.01 

0.32 

0.02 

Te 

-  0. 1  3 

0.13 

0.44 

0.10 

0.42 

0.08 

0.79 

0.09 

Cl 

<  0  5 

0.06 

0.03 

0.10 

0.02 

0  41 

0.04 

yn 

<  0.5 

<  -0.5 

-0.30 

-0.20 

T  ;  on 

Cd  site 

c 

-  0.02 

0.37 

1.32 

0.22 

1.59 

0.20 

1.39 

0.19 

Si 

1  57 

0.65 

>2.0 

>2.0 

>  2.0 

P 

0.16 

OriS 

1.48 

0.26 

1.73 

0.23 

1.52 

0.21 

O 

,  (IT 

0.89 

0.14 

1.22 

0  1 3 

1.03 

0.12 

Te 

0.48 

0.55 

1.60 

0.29 

1.88 

0.23 

1.66 

0.24 

Cl 

.  0  5 

096 

0  17 

1.29 

0.14 

1.09 

0.14 

V,. 

n  * 

0.00 

0.21 

006 

T  „  on 

Te  site 

I.i 

<i  14 

0.2‘t 

1.28 

0.22 

1.15 

0.35 

0.76 

0  25 

Cu 

0.5 

0.54 

0.42 

0.12 

0  52 

0.03 

0.32 

I  on 

Te  site 

Ag 

1  Ml 

0  t2 

1  26 

0.22 

1.21 

0.23 

0.69 

0.20 

Cd 

1  66 

0  34 

1  1  1 

0.26 

1.05 

0.26 

0.85 

0.22 

Ga 

0  48 

O  4M 

0  61 

0.33 

0.55 

0.32 

0  40 

0.30 

Si 

0  07 

0  40 

Oil 

0.36 

on 

0.38 

0.38 

0.72 

Sn 

0  28 

0  4' 

O  15 

0.31 

0  13 

0.24 

0.02 

0.28 
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which  is  comparable  lo  that  due  to  different  host  band 
structures. 

Putting  this  large  uncertainty  in  the  deep  levels  against 
a  band  gap  of  1.6  eV,  we  are  left  with  great  doubts  about 
the  predictability  of  this  oversimplified  theory.  Unfor¬ 
tunately,  the  experimental  means  available  for  identifying 
microdefects  in  semiconductors  are  still  very  limited,  and 
the  ab  initio  band  theory  is  still  not  capable  of  accurately 
predicting  the  energy  levels.  Thus,  there  is  a  great  temp¬ 
tation  to  use  simple  theories  like  the  one  carried  out  here 
to  help  with  the  identifications.  To  illustrate  this  point, 
consider  the  following  examples:  Table  II  shows  that  Li 
on  a  Te  site  has  an  s  level  of  0.14  eV  in  model  (a),  so  one 
may  be  tempted  to  relate  it  to  the  acceptor  state  identified 
experimentally.14  However,  this  is  not  the  hydrogenic  ac¬ 
ceptor  state  on  a  Cd  site,  as  one  might  anticipate.  One 
might  also  want  to  assign  the  y  and  y  gaP  states  for  the 
Te  antisite  p  levels  on  the  Cd  site  found  from  model  (a)  as 
those  seen  in  experiments. I5  lh  Because  of  the  large  uncer¬ 
tainty  in  the  calculation,  these  results  should  be  regarded 
as  suspicious  surprises  rather  than  theoretical  confirma¬ 
tions. 

The  results  presented  here  should  not  discourage  con¬ 
tinued  research  on  the  ETB  approach,  but  improvement  is 
clearly  needed.  Work  ranging  from  universal23,27,30  to 
specific24  31,32  structural  studies  to  our  band  calculations 
and  alloy  studies'  7  indicates  that  the  ETB  type  of  theory 
is  practical  for  both  bonding  properties  and  electronic 
structures.  The  reason  that  ETB  works  well  for  some 
properties,  e  g.,  photoemission  spectra  and  bonding  prop¬ 
erties,  but  not  for  impurity  levels,  is  that  the  former  de¬ 
pend  only  on  the  gross  total  density  of  states,  while  the 
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latter  have  been  shown  to  be  sensitive  to  the  details  of  the 
partial  densities  of  states. 

To  establish  the  credibility  of  ETB  in  defect  studies, 
one  needs  to  look  at  the  problem  more  seriously.  The 
most  difficult  and  yet  important  task  is  to  develop  a  better 
way  for  determining  the  Hamiltonian  matrix  elements. 
Haas  el  al*  and  Harrison27,28  have  suggested  using  the 
atomic  term  values  as  the  diagonal  matrix  elements.  Our 
work1-4  has  suggested  using  a  universal  long-range  in¬ 
teraction  to  improve  the  accuracy  of  the  conduction 
bands.  Several  studies1,27,28,33  have  also  pointed  out  scal¬ 
ing  rules  of  the  matrix  elements.  A  combination  of  these 
ideas  may  lead  to  an  acceptable  model.  Secondly,  both 
the  bonding  and  deep-level  states  of  impurities  should  be 
studied  at  the  same  time  in  order  to  provide  correlated  in¬ 
formation  for  defect  identification.  Finally,  more  realistic 
models  should  be  examined.  Besides  the  substitutional 
site-diagonal  defects,  one  should  consider  the  possibility 
of  interstitial,  paired,  and  even  more  complex  defects. 
One  also  needs  to  deal  with  long-range  impurity  poten¬ 
tials,  possible  charge  shifts,  and  lattice  distortions.  Pro¬ 
gress  in  all  these  areas  can  be  expected  if  the  calculation  is 
constantly  correlated  with  experiments  and  available 
ab  initio  theory. 
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A  formula  for  alloy-scattering-limited  electron  mobility  in  semiconductors  is  obtained  for 
indirect  gap  systems  with  multiple  band  minima.  All  the  input  parameters  needed  are  defined 
explicitly.  The  drift  mobility  of  Si,  Ge,  _  ,  which  has  a  dip  at  x  —0.13  and  a  broader  minimum  at 
x  —0.5  is  calculated  by  adding  alloy  scattering  to  other  scattering  mechanisms  and  correlates  well 
with  the  measured  Hall  mobility. 


The  electron  and  hole  mobilities  in  semiconductors  are 
determined  by  the  band  structure  and  various  scattering 
mechanisms,  predominately  impurity  and  phonon  scatter¬ 
ing.  For  alloys,  the  mobility  is  also  affected  by  disorder  aris¬ 
ing  from  aperiodic  atomic  potentials  and  atomic  positions. 
Many  years  ago,  Nordheim1  and  Brooks3  obtained  an 
expression  for  alloy-scattering-limited  electron  mobilities  in 
metals  and  semiconductors,  respectively.  Brooks’  well- 
known  formula  reads 

,j2irefi4N0 _ 1  ,|, 

fi/>  3*|1  -  x)m*'n  (AE)2jKT  ' 
where  Nu  is  the  number  of  atoms  per  unit  volume,  m*  is  a 
band-edge  effective  mass,  x  is  the  fractional  concentration  of 
one  of  the  species,  and  AE\%  an  energy  parameter  character¬ 
izing  the  alloy  potential  fluctuations.  Although  this  formula 
has  been  widely  and,  to  some  extent,  successfully  used  for 
direct  gap  materials, the  identification  of  the  alloy  disor¬ 
der  parameter  AE  remains  uncertain.  Various  suggestions 
have  previously  been  made  for  AE,  e.g.,  and  band-edge  dis¬ 
continuity*  or  band-gap  differences.7  Any  of  these  simple 
choices  is  bound  to  fail  when  one  applies  Eq.  (1)  to  more 
complicated  indirect  gap  systems  such  as  Si,  Ge,  .  ,  alloys, 
where  one  encounters  conduction-band  minima  transferring 
between  the  X  and  L  points  of  the  Brillouin  zone.  For  exam¬ 
ple,  if  AE  is  taken  to  be  the  difference  in  corresponding  band 
edges,  then  one  finds  that  AE  —0.1  eV  for  the  X  (A  )  valley 
and  -  1.2  eV  for  the  L  valley.  The  values  that  fit  the  experi¬ 
ment  are  about  half  this  value  for  L  and  —  0.5  eV  for  T.*  The 
purpose  of  this  letter  is  to  resolve  the  identity  of  AE  for  indi¬ 
rect  gap  materials. 

Moreover,  there  is  a  problem  with  the  m*  that  enters 
Eq.  ( 1 ).  For  direct  gap  alloys,  the  band-edge  effective  mass  at 
r  naturally  enters  Eq.  |l|.  For  the  indirect  gap  alloys,  the 
effective  mass  is  anisotropic  and  hence  an  appropriate  mass 
must  be  chosen  Previous  authors*  7  have  chosen  m*  to  be 
the  effective  conductivity  mass  m*.  We  shall  show  that  dif¬ 
ferent  masses  enter  for  different  cases. 

The  first  unambiguous  assignment  for  A  E  in  a  direct  gap 
alloy  was  given  by  Hass  el  al .'  To  estimate  the  limiting  elec¬ 
tron  mobility  in  Hg,  Cd,  ,  Te  based  on  a  tight-binding  (TB) 
band  description,  they  defined  AE  to  be  fAE,,  where  f  is 
the  r  fraction  in  the  density  of  states  and  AE,  is  the  difference 
between  the  s  atomic  term  values  of  the  Hg  and  Cd  atoms.  By 
extending  this  approach  to  alloys  with  indirect  gaps  and 


multiple  bands,  we  show  that  all  the  uncertainties  identified 
above  are  resolved.  Our  generalized  Brooks’  formula  will 
then  be  applied  to  Si,Ge,  _,  systems  to  explain  their  ob¬ 
served  mobility.7,8 

Because  Brooks’  formula  has  never  been  derived  expli¬ 
citly  in  the  literature,  we  rederive  it  first  and  then  generalize 
it.  Consider  the  case  of  a  single  band  with  an  isotropic  effec¬ 
tive  mass.  The  dc  electronic  conductivity  based  on  the  linear 
response  theory9  is  given  by  : 

where  the  energy-dependent  off)  in  the  weak  alloy  scattering 
limit  is 

(3) 

D(e)  is  the  density  of  states  (DOS)  per  unit  volume  for  both 
spins,  so  D  (e)  =  2NA  p(e),  with  NA  =  Nq/2  being  the  number 
of  unit  cells  per  volume  (for  the  diamond  structure,  half  the 
number  of  atoms  N0  per  unit  volume)  and  />(e)  being  the  DOS 
per  unit  cell  per  spin.  The  mean  square  velocity  i>3(e)  for 
carriers  with  energy  e  is  given  by 

o3(f)=  Y|t;(k)|3  g[*~fk)1.  (4) 

T  P(f) 

The  scattering  lifetime  for  carriers  with  energy  e,  rif),  is  re¬ 
lated  to  the  alloy  broadening  A  (f )  by  r{()  =  ft/2A  (e),  where 
the  energy  A  (f )  is  the  imaginary  part  of  the  self-energy  in  the 
averaged  alloy  Green's  function.  For  weak  scattering  A  (f  |  is 

A  (e)  =  nx[  1  -  x)  [AE  )3p(e),  (5) 

where  in  a  tight  binding  (TB)  description  AE  is  the  difference 
in  the  term  values  of  the  constituents.  Then  the  mobility  is 
p  t  =  a/ne  with  the  electron  density  given  by 

n  =  2jf[()p[e\de.  (6) 

For  a  nondegenerate  semiconductor,  f\e)  is  the  Boltzmann 
distribution  and  f[e) «  e  "  *’ ,/kT.  Furthermore,  for  a  para- 
bolicbandfjkl  =  tfk  2/2 m*,p(e)  =  |2m*),/3f ''v^r^.then 
all  the  above  equations  can  be  combined  to  arrive  at  Eq.  1 1 ). 

For  a  real  semiconductor  alloy  in  a  TB  description,  the 
alloy  scattering  can  be  characterized  by  two  parameters  J  E , 
and  AEP ,  the  differences  in  s  and  p  term  values  between  two 
substitutional  atoms  Then  an  effective  broadening  is  given 
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by 


TABLE  I  Calculation  parameter*. 


A  (f)  =  \A,p,  +  Appp\/p,  (7) 

where  p, ,  pp  are  partial  density  of  states  (PDOS)  and  A ,  and 
Ap  are  similar  to  Eq.  (5),  with  p  replaced  by  p,  and  pp ,  re¬ 
spectively.  For  Hg,  ,  Cd,  Te,  the  s  disorder  is  predomin- 
ent’ 10  and  one  can  neglect  AE?.  Defining  p,  =f%p  (and 
pp  —  fpp),  one  arrives  at 

d=;jr*(l  - x)(f,AE,)2p. 

Thus  as  was  pointed  out  by  Hass  et  al.,*  f,AE,  plays  the  role 
of  AE  in  this  special  case  where  AEP  can  be  neglected. 

For  an  alloy  with  a  single  indirect  gap  minimum,  one 
has  to  consider  both  s  and  p  contributions  to  the  alloy  broad¬ 
ening  and  the  masses  that  enter  p  and  u2.  Again,  Eqs.  (2H6) 
can  be  combined  to  yield 


Pa 


_ {ertSM*) _ 

[3*(1  - x]mtm:(m*)'n(kT)',2NAf)AE]  +f\AE2p)]' 

(8) 


where  m*  and  m*  are  respectively  the  longitudinal  and  the 
transverse  mass  at  the  band  edge,  and  Nu  is  the  number  of 
equivalent  minima,  e.g.,  6  for  Si.  The  conductivity  mass  m* 
comes  from  averaging  v2  in  Eq.  (2)  and  is  given  by  3(2/m* 
-f  1/mf)  1  Equation  (8)  clearly  identifies  the  masses  and 
the  energy  parameter  that  enter  Brooks’  formula. 

Next  we  consider  a  still  more  complicated  case  where 
the  contribution  to  the  mobility  comes  from  more  than  one 
band.  For  example,  in  Si,  Ge,  _x  the  X  and  L  minima  cross 
ncarx  =0.15."  There  are  now  two  contributions  to  the  net 
conductivity,  so  a  =  lat ,  where  i  is  X  or  L.  The  quantities 
ufle).  D,  (ej,  and  \,{f)  now  take  different  values  for  different 
bands.  The  structure  of  r,  (e|  requires  more  careful  consider¬ 
ation.  The  complication  comes  from  the  fact  that  the  effec¬ 
tive  broadening  A  is  still  given  by  Eq.  (7),  but  p,,pp,  and p 
contain  contributions  from  both  the  bands.  The  proper  ex¬ 
pressions  ar tp  =  i,p,N'„  and pa  =  XJ„ ,  N  'vp ,  where  /  =  X 
or  L,  a  =  s  or p,  and  N  *  =  6,  N  p  =  4.  The  equation  for  A  is 

A  (e)  =  irx[  \  —  x)  I {l,fa.N\,PM)AEa  W 


The  mobility  associated  with  the  /th  band  is  defined  as 

p,  =  a,/(n,e\,  then 


P.  =■ 


netfN., 


3x|l  —  x I  [mrV(2m*|';:l, 


(10) 


ep,\f)fcN’vp,\f\)e 

/,  =  - — - - - "*'  *.  (  11) 

I  (ZL.HWUE.,) 

,t  \  / 


Thus,  the  generalized  formula  no  longer  has  the  explicit 
xand  T  dependences  of  the  original  Brooks' form.  However, 
all  the  quantities  needed— the  masses,  the  scattering  param¬ 
eters  AEC, ,  the  band  gaps,  and  the  fractions/,,  — can  be  eval¬ 
uated  theoretically  without  resorting  to  experimentally  fit¬ 
ted  parameters  To  demonstrate,  we  shall  apply  Eq.  (10)  to 
Si,  Ge,  ,  The  band  quantities  are  obtained  from  our  CPA 
calculation  1 1  We  found  that  the  effective  masses  vary  weak¬ 


Parameter 

Si,  Ge,  ,  systems 

mfl*) 

0.97m„ 

(XI 

0  19m,, 

m»|Z.| 

l.64m„ 

mT\L) 

0.082m„ 

E'.W 

0  8941  +  0.042lx  4  0.169U-’ 

£>l 

0.7596  +  1  0860«  +  0.3306xJ 

A rW 

0.333  +  0.05*  (0<x<0.3| 

0.339  +  0.03jc  |0.3<xc  1.01 

fuM 

0  632  +  0  13* 

ly  with  the  concentration,  so  m f  and  m*  are  assumed  to  be 
constant  and  assigned  the  values  0.97  and  0.19  for  the  X 
minima  and  1.64  and  0.082  for  the  L  minima,  respectively. 
The  calculated  energy  gaps  for  the  X  (A  )  follows  the  func¬ 
tional  form  E  1  =  a  +  bx  +  cx1  and  for  L  is  given  by  E  1 
=  A  +  Bx  +  Cx2.  All  the  parameters  of  our  calculations  are 
listed  in  Table  I. 

To  correlate  the  calculation  with  the  measured  mobili¬ 
ties,  we  need  to  have  an  estimate  of  scattering  rates  l/r„  due 
to  impurities  and  phonons.  A  crude  approximation  is  to  as¬ 
sume  1  /r0  for  a  given  valley  to  be  the  same  as  the  appropriate 
constituent’s  values  and  add  to  it  the  alloy  scattering  rate  1/ 
rA .  Then  the  average  mobility  and  the  mobility  from  the  ith 
minimum  in  the  alloy  are 

*  =  I 

w-r'-itfr'+wr'.  U2) 

p*  is  given  by  Eq.  (10)  and  p°  are  the  measured  drift  mobili¬ 
ties  for  Si  or  Ge.lz  The  drift  mobility,  calculated  from  Eq. 
( 1 2),  is  plotted  as  a  function  of  alloy  concentration  x  in  Fig.  1 . 

For  x <0.05  and  x>0.20,  the  energy  difference  between 
the  X  and  L  edges  is  large  enough  so  there  is  a  negligible 
contribution  to  the  mobility  from  the  higher  minima.  In  the 
Si,  Ge,  _ ,  system,  the  s  scattering  is  predominent.  Because 
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FIG  I  Calculated  (solid  line!  electron  drift  mobility  and  ihe  experimental 
Hall  mobility  idashed  linel  from  Refs  6  and  8  are  plotted  as  a  function  of 
alloy  concentration 
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the  L  edges  have  a  larger  s  content,  alloy  scattering  reduces 
the  average  mobility  substantially  for  small  x.  Even  though 
the  s  content  is  almost  the  same  for  allx>0.20  at  the  X edge, 
the  mobility  still  decreases  to  x=;0.5  as  shown  in  Fig.  1  be¬ 
cause  of  thex(l  —  x)  term  in  Eq.  (10). 

An  interesting  feature  is  obtained  for  the  compositions 
0. 1 3 < jc  <0.18.  The  average  mobility  attains  a  local  mini¬ 
mum  near  x  =  0. 14  and  a  smaller  maximum  at  x  =  0.17. 
This  feature  occurs  because  of  the  X  to  L  crossover. ' 1  For 
x<0. 14,  the  major  contribution  to/z  comes  from  L  minima. 
Near  x  =  0. 14,  the  density  of  states  increases  because  the  X 
and  L  minima  merge.  So  the  alloy  scattering  increases  there 
and  the  average  mobility  decreases.  For x>0. 14,  the  X bands 
have  the  lower  minima.  As  the  5  content  is  small  at  the  X 
minima,  the  reduced  alloy  scattering  increases  the  average 
mobility.  For  larger  values  of  x,  the  x{  1  —  jc)  term  takes  over 
and  the  mobility  varies  as  shown.  The  values  of  measured 
Hall  mobility  in  Si,  Ge ,  ,  systems  are  also  plotted  in  Fig.  1 . 
The  interesting  feature  near  x  =  0. 14  is  clearly  seen.  Since 
the  experimental  drift  mobility  for  Si,Ge,  .  „  is  not 
available  and  the  generalization  of  Eq.  (10)  to  Hall  mobility 
fiH  is  less  clear,  we  present  the  calculated  fiD  and  experimen¬ 
tal  fiH  (Ref  7,8)  here.  While  we  do  not  expect  quantitative 
agreement,  because  can  range  from  1  to  2, 11,14  we  do 

expect  them  to  display  the  same  qualitative  x  dependence.  It 
is  rewarding  to  note  the  similarity  in  the  trend  in  Fig.  1. 
Previous  authors  explained  the  dip  in  the  mobility  curve  by 
including  intervalley  scattering  with  an  arbitrary  adjustable 
coupling  constant  Our  calculations  automatically  include 
that  portion  of  intervalley  scattering  that  results  from  alloy 
disorder  with  a  coupling  constant  set  by  the  atomic  proper¬ 


ties  of  the  constituents  However,  the  additon  of  intervalley 
scattering  mediated  by  phonons  and  impurities  is  expected 
to  increase  the  dip  near  crossover. 

In  summary,  an  expression  for  alloy-scattering-limited 
charge  earner  mobilities  is  derived  for  indirect  gap  alloys 
with  multiple  bands.  This  expression  reduces  to  Brooks'  for¬ 
mula  for  direct  gap  alloys.  The  quantities  m*  and  AE  can  be 
calculated  exactly.  Alloy  scattering  accounts  for  the  ob¬ 
served  mobility  features  in  the  Si,  Ge,  _  „  alloy,  including 
the  anomaly  near  the  L  to  X  (d  )  crossover. 
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Binding  Energy  and  Spectral  Width  of  Si  2 p  Core  Excitons  in  Si*Gej  _  x  Alloys 
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and 
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A  calculation  is  presented  to  explain  the  anomalous  experimental  behavior  of  the  Si  2 p  core- 
exciton  binding  energy  and  linewidth  in  SixGei-x  alloys.  The  observed  minimum  in  the  linewidth 
near  x  =  0.15  can  be  explained  as  the  result  of  a  competition  between  intrinsic  broadening  due  to 
screening  and  extrinsic  alloy  broadening  For  pure  Si,  the  binding  energy  is  estimated  to  be 
0  15  ±  0  05  eV  and  the  width  is  shown  to  be  smaller  than  that  observed  at  x  —  0.15. 


PACS  numbers  71  55  Fr,  71  35  +z,  78.70  Dm 

Until  1984,  the  Si  2p  core  exciton  was  believed  to 
have  an  anomalously  large  binding  energy. I‘,°  Later, 
Newman  and  Dow11  proposed  a  radically  different  pic¬ 
ture  in  which  the  Si  2 p  core  exciton  is  in  fact  a  reso¬ 
nance  with  a  negative  binding  energy.  They  further 
predicted  that  the  exciton  binding  energy  remains  neg¬ 
ative  throughout  most  of  the  SixGe|_x  alloy  composi¬ 
tion  range,  except  near  x  —  0.20  where  it  becomes 
positive  In  a  recent  experiment,13  Bunker  et  al.  found 
an  anomalous  sharpening  of  the  exciton  spectra  near 
x-0. 15,  the  data  were  interpreted  to  support  the 
Newman-Dow  point  of  view  Yet  the  most  recent  ex¬ 
periment1’  still  suggests  a  positive  value  for  the  bind¬ 
ing  energy  Eb  in  silicon 

In  this  l  etter,  we  present  a  calculation  that  offers  a 
plausible  resolution  to  the  above  problem  In  our 
theory,  the  calculated  Si  2p  core-exciton  binding  ener¬ 
gy  Eb  (  x  i  ,ind  the  linewidth  A  (  x )  in  SixGei-x  alloys 
are  sensitively  dependent  on  the  parameter  Eb  A 
comparison  of  the  calculated  A (x)  with  the  experi¬ 
ment13  suggests  a  positive  value  0  15  ±0  05  for  Eb 
The  anomalous  experimental  spectrum13  near  x  =015 
is  explained  as  a  result  of  a  competition  between  an  in¬ 
trinsic  broadening  A  /  due  to  screening  and  an  extrinsic 
alloy  broadening  A  4  In  the  present  theory,  there  is  no 
need  to  suppose  that  the  exciton  suddenly  changes  its 
character  from  an  extended  effective-mass-like  state  to 
a  deep  localized  state 

We  need  to  calculate  Eb  and  A  =  A/  +  A^  as  a  func¬ 
tion  of  alloy  concentration  x  The  calculations  are 
based  on  a  quantitative  coherent-potential-approx- 
imation  H  P  \ >  band  structure  Details  of  the  CPA 
calculations  will  he  presented  elsewhere  Below,  we 
briefly  discuss  a  Green's-function  method  for  calculat¬ 
ing  Eb  and  A  , 

f he  one-particle  effective  Green's  (unction  in  C  PA 
takes  the  lorm 

(,<  /.)-[/  H  It  /  )  1  (I  » 

where  //  is  (he  virtual-crystal  approximation  Hamil- 

120 


Ionian  and  1(E)  is  the  self-energy.  The  site-diagonal 
Green's  function  is  denoted  as 

£„(£)«  <*.|G< £)!*.).  (2) 

where  <t>„  is  a  localized  orbital  of  specified  symmetry. 
Here  we  only  need  to  consider  a  =  sfor  A]  symmetry. 
The  corresponding  function  in  pure  Si  is  denoted  as 
£,°(£).  Following  the  theoretical  treatment  of  deep 
substitutional-impurity  levels,®  the  core-exciton  level 
for  pure  Si  is  determined  by 

£,°(£)  =  (  y-E,s,)~\  (3) 

where  £,Sl  is  the  site  potential  seen  by  an  s  electron  in 
bulk  silicon,  and  V  is  a  central-potential  parameter 
For  a  chosen  value  of  Eq.  (3)  can  be  solved  for  £, 
and  vice  versa  Then  £°=£f°-£.  where  £f°  is  the 
conduction-band  edge  in  pure  silicon.  Because  of  the 
uncertainties  in  the  value  of  experimental  Eb  and 
theoretical  V,  we  treat  Eb  (or  V)  as  a  parameter.  The 
binding  energy  Eb  in  a  SixGei_x  alloy  can  be  calculat¬ 
ed  by  solving 

£,(£)»lL-£,-  I.ff)]"1.  (4) 

where 

E,  =  x£,Sl  +  ( 1  -  x) £,0t  (5) 

Then  Eb  is  given  by 

£»  =  £,-£  (6) 

The  calculated  values  of  the  conduction-band  edge 
and  the  exciton  level  measured  relative  to  the  top  of 
the  valence  hand  are  plotted  in  Fig  1  The  band  gap 
increases  with  x  with  a  slope  discontinuity  at  =0.15. 
The  dashed  lines  a.  b,  and  c  represent  exciton  levels 
obtained  with  Eb  =  0  1,  0  15,  and  0  30  eV,  respective¬ 
ly  The  binding  energy  Eb  is  also  an  increasing  func¬ 
tion  of  x,  with  a  slope  discontinuity  near  x  =  0  15 
The  CPA  introduces  a  slight  bowing  in  £f  and  Eb 
Stnnati 10  has  calculated  the  variation  of  A;  with  Eb 
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FICi  I  Variation  of  the  band  gap  (solid  line)  and  the  Si 
2p  core-exciton  level  with  a  in  Si^Gei.,  alloys.  The  energy 
is  measured  from  the  top  of  the  valence  band  The  dashed 
curves  a.  b.  and  <  represent  exciton  levels  calculated  with 
£*°  -  0  1,0  15,  and  0  3  eV.  respectively 

by  replacing  the  short-range  Coulomb  potential  with  a 
spherical  square  well  of  variable  depth  and  a  screened 
Coulomb  tail  Strmati's  results  can  be  used  to  estimate 
A,  corresponding  to  the  calculated  Eb  A/  decreases 
rapidly  with  Eb.  then  saturates  for  larger  Eb. 

The  contribution  to  the  natural  linewidth  from  the 
alloy  broadening  is  calculated  by  a  consideration  of  the 
electron  part  of  the  exciton  wave  function,  »//,.  The  t bs 
is  expanded  in  a  linear  combination  of  the  s  part  of  the 
conduction-band  wave  functions  <t>s„ 

|(k).  (7) 

i»,  s 

We  found  that  alloy  scattering  is  only  moderate  and  s 
scattering  is  dominant,  thus,  the  alloy  broadening 
A  t  (  £1  is  well  approximated  by 

A  4  (  E)  --  \  ( I  -  vin/fm /,<£)  (8) 

where  h,  is  the  difference  between  £,Sl  and  E^e. 
Hence,  the  alloy-broadening  contribution  to  A  is  relat¬ 
ed  to  the  alloy  broadening  of  the  band  slates, 
A  „  <  k .  F,  l 

a,  -  —  £  <  •/< ,  <  k  >  ;a  4  </•;)[./*,(  k  >  > 

^  k 

■■■{  no.fk.f 

k  n 

=  J"p,(  /.'A  (  /.></£ 

-  v  (  I  v)h,:rr  f„;(E)<JE  (9) 

The  integral  in  Fq  <9)  is  evaluated  numerically 
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FIG  2  Variation  of  A  (solid  lines)  and  A,  (dashed  lines) 
with  xfor  three  Eb  values. 

The  calculated  A,  which  is  the  sum  of  A,  and  A^,  is 
plotted  against  v  in  Fig.  2  for  three  values  of  Eb  In  all 
three  panels,  the  dashed  curve  represents  A /  and  the 
solid  line  represents  A.  It  is  seen  from  Fig.  1  that  the  V 
exciton  level  follows  the  X  edge  of  the  conduction 
band.  Fience  the  binding  energy  Eb,  relative  to  the 
conduction  band  edge,  remains  almost  constant  (for  a 
given  Eb)  until  the  minimum  switches  from  the  V 
edge  to  the  L  edge  Because  of  the  change  in  the  slope 
of  Eg,  Eb  decreases  rapidly  when  L  becomes  the 
minimum.  Correspondingly,  A/  vanes  slowly  until  the 
X  to  L  crossover  and  then  increases  rapidly.  This 
feature  is  clearly  seen  in  Fig.  2. 

For  Eb  =015,  the  A/  and  A^  are  comparable  near 
v  =  0.50,  and  A/  dominates  for  all  small  x  and  large  \ 
These  two  competing  mechanisms  give  a  relative 
minimum  near  ,v  =  0.15,  a  broader  maximum  near 
v  =  0  50.  and  a  smaller  minimum  for  pure  silicon  As 
Eb  is  decreased,  the  relative  minimum  is  shifted  to 
larger  v.  e  g  ,  the  minimum  shifts  to  v  =  0  20  for 
Eb  =  0  10  eV  For  Eb=  0  15  eV.  the  position  of  the 
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relative  niinmium  is  in  agreemeni  with  the  experi¬ 
ment  (Bx  measuring  the  relative  width  at  a  =  0  15 
in  that  x  =  I,  one  can  make  a  better  estimate  of  Fb  • 
To  correlate  the  theory  with  experiment,  the  calculated 
1/A*’  is  compared  with  the  measured12  (AmI_i(c/m/ 
dh.  )mjl  in  Fig  3,  where  Am  is  the  edge  step  and 
<  dti  'Jf') m4V  arc  the  maximum  values  of  the  dertvatise 
of  absorption  spectra  with  respect  to  photon  energy 
Because  the  experimental  values  arc  given  in  arbitrary 
units,  the  values  arc  normalized  to  agree  at  a  =0  5 
The  observed  anomalous  behavior  near  v  =  0  15  and 
the  qualitative  x  dependence  in  that  region  is  clearly 
replicated  by  the  theory  However,  the  calculation 
predicts  a  larger  maximum  at  x  =  1  It  would  be  in¬ 
teresting  to  have  experiments  that  cover  the  entire 
range  of  x  to  further  test  this  prediction 

For  larger  values  of  A*1,  the  calculated  F.b  is  also 
large  and  hence  A/  decreases  slowly  with  x  Because 
the  broadening  is  determined  mainly  by  A<,  the 
linewidth  is  expected  to  be  small  for  x  =0  and  x  -  1 
o  > .  this  occurs  for  Fb  -  fl  30  eV  For  negative  values 
of  Fh  remains  negative  for  all  values  ol  x  Accord¬ 
ingly.  the  linewidth  is  broad  lor  all  x.  and  there  would 
be  no  such  anomaly  as  in  Fig  3 
The  calculations  presented  in  this  letter  are  slightly 
different  from  alchemy  approximations.*  We  treat  the 
central-cell  potential  t  as  a  parameter  and  narrow  its 
range  from  other  considerations  We  examine  values 
of  l  /,S|  ol  -8  49,  7  09.  and  6  56  eV.  corre¬ 

sponding  to  /•.’(!'  values  of  0  30.  0  15.  and  0  10  eV, 
respect  ely  II  the  strict  alchemv  approximation  were 
taken,  the  value  ol  I  /,S|  would  he  /.f  /.,s  =  -  4  59 

eV  in  the  tight-binding  approximation,  and  a  negative 
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£®  (  —  -0  10  eV)  would  be  obtained  When  long- 
range  interactions  are  included,  however,  the  above 
resonance  state  becomes  a  shallow  donor  level,  which 
is  the  experimental  situation  for  a  P  impurity  in  Si. 
Our  results  suggest  that  V  for  core  excitons  in  Si  is 
deeper  than  those  implied  by  alchemy  approximations. 
However,  if  we  use  the  alchemy  approximation  as  a 
means  of  scaling,  the  value  of  V  for  Ge  3 p  core  exci¬ 
tons  should  be  deeper  than  that  for  Si  2/icore  excitons 
Hence,  the  curve  corresponding  to  £4°=0.30  in  Fig  2 
should  be  a  reasonable  estimate  for  Ge  3/rcore-exciton 
binding  energy  in  alloys.  Therefore,  we  do  not  expect 
to  see  an  anomalous  behavior  of  A  in  alloys  for  this 
case. 

In  summary,  the  present  calculations  of  the  St  2 p 
core-exciton  binding  energy  and  linewidth  suggest  that 
the  exciton  level  is  about  0  15  ±0  05  eV  below  the 
conduction-band  edge  for  pure  St.  It  follows  the  X 
edge  for  x  >0  15  in  the  Si^Gci^,  alloys,  and  Eb  may 
eventually  reach  zero  in  the  dilute  limit  v  — ■  0.  Our 
value  for  Eb  represents  the  lower  end  of  the  previous 
measurements,1'6  but  is  in  good  agreement  with  a  re¬ 
cent  experimental15  value  of  0  1 20  ±  0.03  eV  By  con¬ 
sidering  the  intrinsic  linewidth  and  the  alloy  broaden¬ 
ing,  we  can  explain  the  observed  relative  minimum  in 
the  linewidth  near  a  —  0. 1 5.  without  requiring  a  sud-  % 
den  change  of  the  exciton  character  On  the  basis  of 
this  calculation,  we  expect  the  corresponding  width  in 
pure  Si  to  be  even  smaller  than  that  observed  near 
x  =  0  15  We  further  argue  that  the  binding  energy  of 
Cie  3/i  core  excitons  should  be  larger  than  that  of  Si  2 /> 
core  excitons  anil  there  should  be  no  anomaly  in  the 
Ge  3/i  linewidth  in  alloys. 
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Harrison’s  bonding  theory,  the  valence  force  field  (VFF),  and  an  elastic  continuum  are  combined 
in  a  study  of  the  substitution  energies  A,  and  local  (first-shell)  bond  lengths  dx  of  isoelectronic  im¬ 
purities  in  semiconductors.  Explicit  expressions  for  A,  and  d ,  are  derived,  which  enable  us  to  ab¬ 
sorb  measured  elastic  constants  into  the  calculation  and  to  study  the  chemical  effects  arising  from 
differences  in  the  covalent  radii  and  polarities.  Several  models  based  on  VFF  alone  are  also  derived 
for  comparison.  The  full  theory  and  at  least  five  VFF  models  are  found  to  produce  impurity  bond 
lengths  in  excellent  agreement  with  experiment.  The  substitution  energies  are  shown  to  provide 
good  estimates  of  the  mixing  enthalpies  ft  of  pseudobinary  alloys  and  to  predict  miscibility  gaps 
properly.  The  chemical  shifts  in  O  are  found  to  be  negative  for  most  cation  alloys  but  positive  for 
anion  substitutions. 


1.  INTRODUCTION 


n.  IMPURITY-SUBSTITUTION  ENERGY 


The  discovery  of  a  bimodal  distribution  of  the  nearest- 
neighbor  bond  lengths'  in  Ga^In,  _x  As  has  sparked  con¬ 
siderable  interest  in  the  bonding  nature  of  semiconductor 
alloys.2-6  This  finding  has  changed  the  conventional  pic¬ 
ture  of  the  alloy  crystal  bond  configuration,  which  has 
far-reaching  implications  about  the  electronic  structure, 
structural  stability,  and  thermodynamics  of  these  materi¬ 
als.  Because  of  the  complexity  of  both  the  structural  and 
the  potential  disorder  in  these  alloys,  ab  initio  band- 
structure  techniques  have  not  yet  evolved  to  a  stage  suit¬ 
able  for  direct  calculations.  Therefore,  we  have  extended 
Harrison’s  bonding  theory7  to  study  the  alloy  structural 
properties.’  6  In  this  paper,  we  apply  an  intermediate  ver¬ 
sion  of  the  theory  to  the  dilute-limit  case  of  an  isoelect¬ 
ronic  impurity. 

A  particularly  useful  application  of  the  theory  is  its 
perturbation-expansion  form,  in  which  measured  elastic 
constants  are  incorporated  to  obtain  accurate  results. 
This  form  is  also  useful  for  comparison  with  other  previ¬ 
ously  published  models3'8,9  that  are  based  on  the  valence- 
force-field  (VFF)  (Ref.  10)  model  alone.  Thus,  all  the 
factors  influencing  bond-length  relaxation,  e.g.,  strains, 
boundary  conditions,  and  chemical  effects,  can  be  studied. 
’» he  ability  to  incorporate  the  chemical  effects  is  one  ma¬ 
jor  difference  between  this  theory  and  other  VFF  models. 

The  remainder  of  the  paper  contains  the  following  sec¬ 
tions:  Sec.  II  describes  a  theory  for  calculating  impurity 
substitution  energies.  Section  III  casts  the  theory  into 
perturbation  form  and  combines  it  with  a  valence  force 
field  and  an  elastic  continuum.  Several  VFF  models  are 
derived  in  Sec.  IV.  The  modifications  of  numerical  re¬ 
sults  due  to  chemical  effects  on  local  bond  lengths  and  al¬ 
loy  mixing  enthalpies  are  summarized  and  discussed  in 
Sec.  V.  Conclusions  are  drawn  in  the  last  section,  Sec.  VI. 


Consider  the  problem  of  substituting  an  isoelectronic 
atom  A  for  a  B  atom  in  a  zinc-blende  compound  BC  (e.g., 
In  substitutes  for  Ga  in  GaAs,  as  shown  in  Fig.  1).  In 
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FIG.  1.  A  sketch  of  the  flattened  picture  of  a  single  impurity 
A  in  a  BC  compound.  The  positions  and  displacements  for 
those  atoms  labeled  are  used  in  Appendix  B. 
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general,  the  bond  lengths  du  d2,  di, .  . .  for  the  first-, 
second-,  and  third-shell  bonds  surroundjng  the  impurity 
are  different  from  the  equilibrium  values  of  either  the 
pure  BC  compound,  denoted  as  d,  or  the  “impurity”  com¬ 
pound  AC,  denoted  as  d /  =  d(  1  —  &o)-  ^  A  starts  being  a 
free  atom  and  B  also  ends  being  a  free  atom,  then  the  en¬ 
ergy  difference  between  the  final  and  initial  states  is  de¬ 
fined  as  the  substitution  energy  and  is  given  by 

A,  =(£def +  £»  )  — (Epurt +E4  )  ,  (1) 

where  eA  and  es  ant  free-atom  energies  for  A  and  B, 
respectively,  E& f  is  the  total  energy  of  the  semiconductor 
with  a  defect  as  sketched  in  Fig.  1,  and  E ^  is  that  of  the 
pure  host  BC  crystal.  Equation  (1)  can  be  written  as 

A,  =  I  £&{  —  Eft,)  —  (tA  £g  )  -hlEjii  £pure  1 

=  A /t  —  ( e  a  —  £g)  +  Aa,  >  (2) 

where  we  have  added  and  subtracted  a  term  £*,,  which  is 
the  total  energy  of  a  distorted  BC  compound  with  all  the 
atoms  held  at  the  positions  specified  in  Fig.  1,  except  that 
the  central  atom  is  a  6  atom.  Clearly,  A*  =£dtf  — ^dis *s 
a  replacement  energy,  and  the  distortion  energy 
Ad,, =  £  d,,  —  £  pure  >s  the  energy  required  to  deform  a  pure 
BC  crystal  from  its  equilibrium  lattice  to  that  specified  in 
Fig.  1.  Ak  contains  all  the  chemical  terms  that  arise  from 
different  bond  lengths  and  polarities  between  AC  and  BC 
bonds 

A„  can  be  treated  most  easily  by  Harrison’s  bonding 
theory.7  In  this  theory,  the  energy  per  bond  relative  to  the 
vacuum  state  is 

£*  =2 -f  +  6eJ,  +-6t:m  ,  (3) 

where  is  the  energy  of  the  bonding  state  constructed 
from  the  two  hybrid  orbitals  facing  each  other  along  the 
bond  direction 


^A£Ck)-(V\  +  V\)'n  , 


(4) 


with  e*  and  rj,  bang  (he  energies  of  the  anion  and  cation 
hybrid  orbitals,  respectively.  The  antibonding  energy  efl 
has  the  same  form  as  in  Eq.  (4),  except  with  a  plus  sign. 
I’;  i'  called  the  covalent  energy,  which  is  the  total  elec¬ 
tronic  Hamiltonian  matrix  element  between  the  two  hy¬ 
brids  in  question,  and  the  polar  energy  V}  is  the  differ¬ 
ence  V\-  -!  (e*  e*).  The  e^  and  E„i,  the  metallization 

energies,  are  the  shifts  of  the  bonding  level  caused  by  in¬ 
teractions  with  the  neighboring  antibonding  states,  where 
*  and  indicate  whether  the  common  adjacent  atom  is 
a  cation  or  anion  For  example,  this  term  for  an  AC  bond 
labeled  by  d ,  111  Fig  1  due  to  an  antibonding  state  labeled 
by  d ;  is  given  by 


1 2.  n 


i  t/*  ■  1 ) ! 7  I  f/.*(2) 


Cs'l) 


ff(2) 


(5) 


where  4  and  B  denote  AC  and  BC  bonds,  respectively, 
tp '.  with  t  j  and  e*  being  the  s-  and  p- term 
values  of  the  common  adjacent  atom  C.  lb{  1  >  is  the 
probability  amplitude  of  finding  an  electron  in  the  hybrid 
orbital  of  the  C  atom  in  the  bonding  state  of  an  AC  bond 


with  a  bond  length  d  t,  whereas  t//(  2)  is  the  correspond¬ 
ing  probability  amplitude  for  the  antibonding  state  of  a 
BC  bond  of  bond  length  d2.  Finally,  V0  is  a  repulsive 
pair  potential  required  to  prevent  the  crystal  from  collaps¬ 
ing  and  to  guarantee  a  correct  equilibrium  bond  length. 
The  local  perturbation,  Eq.  (S),  is  applicable  because  the 
square  root  of  the  numerator  is  much  smaller  than  the 
separation  between  the  antibonding  and  bonding  levels 
and  the  valence  band  is  completely  filled,  so  the  interac¬ 
tion  between  the  bonding  states  only  spreads  the  eb  levels 
into  bands  without  affecting  the  center  of  gravity  of  the 
occupied  states. 

The  replacement  energy  A*  of  Eq.  (2)  can  now  be  writ¬ 
ten  explicitly: 

A*  =  4[  2tg( l)+^(l)+6ei^(l,l)  +  6 e*£( 2, 1 ) 

-  2ef(  1 ) -  Fg(  1 )  - 6e*$<  1,1)-  6e*£( 2, 1 ) 

1,2)— 6e*£<  1,2)]  •  (6) 

The  distortion  energy  A*,  of  Eq.  (2)  now  involves  only  BC 
bonds  of  different  bond  lengths.  It  can  be  treated  with  ex¬ 
actly  the  same  procedure  for  any  given  set  of  bond-length 
distributions.  Thus,  a  straightforward  energy  minimiza¬ 
tion  procedure  can  be  carried  out.  The  accuracy  of  this 
procedure,  however,  depends  in  turn  on  the  accuracy  of 
scaling  rules  for  V2  and  V0  and  the  input  parameters.  At 
present,  Harrison’s  model7  with  V2  a  1/d7  and  F0  <*  1/d4 
and  his  universal  parameters  are  only  semiquantitative. 
We  are  improving  the  quantitative  nature  of  the  theory  so 
that  the  full  theory  will  yield  accurate  predictions  of  the 
structural  and  thermal  properties  of  semiconductor  de¬ 
fects  and  alloys. 


III.  PERTURBATION  EXPANSION, 
VALENCE  FORCE  FIELD, 

AND  ELASTIC  MEDIUM 

As  pointed  out  earlier,  a  perturbation  expansion  of  the 
theory  is  instructive.  This  is  feasible  because  the  differ¬ 
ences  EV2  =  V2(AC)-V2(BC)  and  AF,  =  V}(AC) 
—  F,(5C)  are  small  compared  to  each  individual  value  for 
many  of  the  isoelectronic  impurities  in  III-IV  and  Il-Vl 
compounds.  To  this  end,  Eq.  (6)  is  rewritten  as 

7  As  =  EbA(  1 )  -  £**(  1 )  +  6[  e*£<  2, 1 )  -  e4;£<  1 , 1 )] 

-6[e4;"  (2,1  )  —  £*"(  1,1)] 

+  6(  £*,*(  1,2)  —  e*’®  ( 2,2)] 

- 6{ e®’ 1,2)  —  £®',c( 2, 2 ) J  ,  (7) 

where  £^(1>  and  Eb(\)  are  energies  per  bond  in  Eq.  (3) 
for  AC  and  BC  compounds,  respectively,  with  the  relaxed 
bond  length  dt=d(l— 6).  The  difference  between  these 
energies  Eb  and  the  corresponding  values  at  their  respec¬ 
tive  equilibrium  bond  lengths  dt  and  d  are  just  the  strain 
energies  per  bond  in  uniform  deformation: 

£**(1)  =  £*<'(<#,  I +  2F1*, «/,(</, -d,)2  , 

(8) 

£*(!>-  E”id)  f  2v  3 /Mid,  d)2  , 
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where  B/  and  B  are  the  bulk  moduli  for  the  impurity  AC 
and  host  BC  crystals.  The  rest  of  the  terms  in  Eq.  (7)  are 
all  due  to  changes  in  tm  caused  by  the  differences  AK3 
and  AF2.  We  shall  use  Harrison’s  scaling  rules  to  deduce 
them.7  Expanding  Eq.  (7)  to  second  order  in  AK3  and 
Ad~d2-  d\,  we  write 

6Ui',c(2, 1 )  —  £m,c(  1  •  1  >1 

=/,  Ad-g/AF3+/t/(Ad)2-IF/AdAFj  +  l//(AFJ)2 , 

(9) 

where  //,  gi,  and  so  on,  are  appropriate  derivatives 
evaluated  for  the  impurity  crystal  AC.  When  similar  ex¬ 
pansions  are  made  for  the  rest  of  the  terms  in  Eq.  (7),  it 
becomes  [with  d,  —d ( 1—6),  dj=d(  I  — 5o>] 

7  A*  =  A£*  +  (/,  -f)Ad  —  (g;  — g)A  V}  +  ( h,  +  h )< Ad  )2 

-  (  Wt  +  W)Ad  A  V,  +  ( U,  +  UK  A  Kj  )2 

+  2\/3fl/d3(8-80)2-2V/3&/38J  ,  (10) 

where 

A  Eb=E?(d,)-E?(d)-iuA-EB)  (11) 

is  just  the  difference  in  the  binding  energy  per  bond  be¬ 
tween  the  “BC’  and  “AC’  crystals.  In  Eq.  (10)  the  coef¬ 
ficients  f,g  without  a  subscript  are  those  for  the  host  BC 
system.  It  is  convenient  to  define  an  excess  energy 
AE  —  A,/4  —  AEb,  which  is  the  extra  energy  per  bond  re¬ 
quired  for  the  impurity  substitution  over  and  above  the 
binding-energy  difference  between  the  BC  and  AC  crys¬ 
tals.  The  binding-energy  difference  accounts  for  much  of 
the  substitution  energy;  however,  the  correction  measured 
by  the  excess  energy  can  be  significant.  The  excess  energy 
results  from  strain  energies  and  chemically  driven  charge 
redistributions  around  the  defect.  Using  Eqs.  (2)  and  (10) 
and  defining  F  —fi—f  and  G  =g,  —  g,  we  can  write  A £ 
up  to  second  order  in  A  V}  and  A d  as 

A  £  =  2  VlB,d  3(  6  -  80)2  -  2 V^Bd  J82  +  F  Ad  -  G  A  F3 

+  H(\d)2+W  Ad  AVi  +  U(\V2)2  +  jA*,  ,  (12) 

where  H  =8,  +h,  W  =  w2  +w,  and  V  =  ut  +  u. 

To  treat  the  distortion  energy  A*,,  we  divide  the  crystal 
into  two  regions.  Inside  a  sphere  of  some  radius  R  mea¬ 
sured  from  the  impurity,  the  strain  energy  is  taken  to  be 
the  valence-force-field10  value: 

AdS'-  ~2a,[A(d,d,)]2 

M  / 

+  ~il'  2  /MAldfdy)]2,  (13) 

8d  t  ) 

i 

where  i  sums  over  all  the  bonds  inside  R  and  the  pairs  in 
the  0  terms  include  those  that  have  adjacent  atoms  inside 
R  and  on  the  boundary.  The  parameters  a  and  0  are 
force  constants  to  be  considered  later  A(d-dy)  =  d,dy 
C'-d'l"  measures  the  change  of  the  dot  product  be¬ 
tween  bond  vectors  due  to  distortions.  Outside  R  we  as¬ 
sume  an  elastic  continuum  with  radial  displacements  in¬ 


versely  proportional  to  the  square  of  the  distance  from  the 
center.  The  elastic  energy  in  this  medium  can  be  shown 
to  be  (see  Appendix  A) 

A^"  =  RCu2 ,  (14) 

where  the  effective  shear  coefficient  is 

C  =  4ir[0.4<C„-C12)+  1.2C*] 

and  u  is  the  displacement  at  R.  In  view  of  the  fact  that 
the  bonds  dx  and  d2  are  coupled  through  the  chemical 
terms  in  Eq.  (10),  the  smallest  logical  radius  R  is  the 
second-shell  atomic  distance,  namely  R  =2v/2d/v/3. 
Atoms  on  this  boundary  have  displacements  of  the  forms 
u=d(y,y,0)/t/5,....  Thus,  u=V2rd/s/3  and  the 
elastic  energy  in  the  continuum  is 

a (15) 

The  distortion  energy  represented  by  Eq.  (14)  contains 
six  different  contributions  (see  Appendix  B):  the  bond¬ 
stretching  energy  of  the  four  first-shell  bonds  6a6 2d2,  the 
0  terms  from  the  first-shell  bonds,  082d2,  the  a  terms 
from  the  second-shell  bonds,  2a(6+2y)2d2,  the  0  terms 
between  the  first-  and  second-shell  bonds,  20(8+y)2</2, 
the  0  terms  among  the  second-shell  bonds,  208 2d2,  and  fi¬ 
nally  the  0  terms  between  the  second-shell  bonds  and 
those  in  the  continuum, 

|0rf2[(36+ A.ir)2+(6+A.2r)2]  , 

where  X,=40v^/(19vT9)  and  A.2  =  2-40v'2/(HV'Tl). 

To  assemble  all  the  contributions  to  Eq.  (12),  we  need  to 
consider  the  assignments  of  the  elastic  constants  and  the 
force  constants  a  and  0  in  VFF.  While  the  experimental 
values"  of  Cn,  C12,  and  Cu  can  be  used  for  the  elastic 
constants,  a  and  0  have  to  be  deduced.  If  Martin’s  origi¬ 
nal  procedure10  (also  followed  by  Martins  and  Zunger3)  is 
used,  then  Eq.  (13)  alone  will  not  produce  the  correct  (ex¬ 
perimental)  bulk  moduli.  There  are  small  corrections  due 
to  Madelung  terms,  which  are  hard  to  treat  in  the  case  of 
non  uniform  distortions.  A  simpler  procedure  is  adopted 
here.  We  use  the  experimental  bulk  moduli  for  0,  and  B 
in  Eq.  (12)  and  experimental  elastic  constants  to  calculate 
C  of  Eq.  (14!  and  then  force  a  and  0  in  the  VFF  to  pro¬ 
duce  the  correct  bulk  moduli  B  and  shear  coefficients 
C||—  Ct2.  Such  an  approach  is  also  consistent  with 
Harrison’s  bonding  theory7  and  other  approaches  in  which 
the  Coulomb  forces  are  automatically  incorporated  in  the 
band  and  bond  energies,  and  do  not  need  to  be  redundant¬ 
ly  treated.  With  our  procedure,  the  bulk  modulus  is  sim¬ 
ply  related  to  the  force  constants  by  B  =(3a  +  0)/(4v/3d). 
Table  I  lists  our  a  and  0  values.  We  want  to  point  out  in 
advance,  however,  that  the  numerical  results  deduced 
from  our  sets  and  those  of  Martin  of  a  and  0  do  not  in¬ 
troduce  differences  more  than  the  present  experimental 
uncertainties  in  the  local  bond  length  (-0.01  A)  and 
the  mixing  enthalpies  (  >0.5  kcal/mole). 

Using  the  above  procedure  and  adding  all  contributions, 
the  excess  energy  per  bond  from  Eqs.  (10)  and  (12)  is  the 
full  perturbation  theory  (FPT)  result 
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TABLE  1.  The  bond  lengths  d  (in  A),  valence  force  constants  a  and  f)  (N/m).  shear  coefficients  C  of 
the  continuum  (in  10"  ergs/cnr1),  melting  temperatures  Tm  (K),  and  Liedermann  ratios  Xm  for  the  com¬ 
pounds  used  in  this  paper. 


Compound 

d 

a 

P 

C 

Tm‘ 

Xm 

A1P 

2.367 

44.323 

8.068 

122.396 

1773 

0.070 

AlAs 

2.451 

40.849 

8.717 

112.695 

1873 

0.073 

AlSb 

2.656 

34.073 

6.900 

85.351 

1323 

0.062 

GaP 

2.360 

44.764 

10.737 

145.921 

1510 

0.064 

GaAs 

2.448 

39.235 

9.159 

121.844 

1738 

0.071 

GaSb 

2.640 

31.876 

7.347 

89.372 

985 

0.055 

InP 

2.541 

40.363 

6.543 

91.785 

1343 

0.059 

InAs 

2.622 

33.203 

5.752 

78.816 

1215 

0.061 

InSb 

2.805 

28.557 

4.891 

60.721 

798 

0.049 

ZnS 

2.342 

40.429 

5.273 

89.272 

2123 

0.081 

ZnSe 

2.454 

32.200 

4.562 

82.687 

1788 

0.080 

ZnTe 

2.637 

29.445 

4.659 

62.430 

1511 

0.071 

CdTe 

2.806 

26.569 

2.722 

38.453 

1371 

0.067 

HgTe 

2.798 

26.396 

2.746 

40.363 

943 

0.056 

•Reference  30. 


AF  =  (  3a/(6-6„)2/2+fr(6-60)2/2+-a(6+-2y)2/2 
+  0<  6  +  y  >2/2  +•  062/4  +  p[  38  +-  )2/8 

+  p{ 6  +  XlY  )74 ]d 1  -f  /2Cy!rf5/( 3 Vl )  +  A £ch  , 

(16) 

where  the  chemical  contnbution  is  written  as 

A £ch=Fch  Ad  +  //(Ad)2  +  A£, 

=  }Fch(2fe  +  rW  +  ?W(25  +  y)2d2 

+  [f/(AK3)2-GAFj]  ,  (17) 

where  AFch  =  F-  W  AF3  and  A£,  =  l/<  AK3)2-G  AK3. 
Flh  is  a  chemical  force,  which  when  it  is  positive  tends  to 
push  the  C  atom  away  from  the  impurity  atom  A.  This 
force  arises  from  the  difference  in  the  bond  tensions  in¬ 
duced  between  the  AC  and  BC  bonds  adjacent  to  C  be¬ 
cause  the  neighboring  antibinding  states  are  different 
from  those  of  their  respective  host  states.  A Ep  is  due  to 
the  difference  in  the  polarities  AF3  alone  and  is  indepen¬ 
dent  of  the  displacement  Finally,  H  can  be  regarded  as  a 
chemically  induced  force  constant,  which  when  it  is  posi¬ 
tive  tends  to  restrain  the  lattice  from  distortion  and  in¬ 
creases  the  elastic  energy. 

The  equilibrium  requirements  0(A£)/36  =  0  and 
3(AF)/0y  =  O  then  lead  to  the  solution  y  =  g6,  and  6  is 
given  by 

6  i  &o+  )/|  1  +  f «( 1  -  2Q)  +  0(17/4  -  kQ) 

+  167/(1  2£))/9]/(3a,+0,l|  ,  (18* 

where  the  constant  X  is  1  +  3A.t  /4  +  k2/2,  and 


%  =  -4Fch(l-G/2)/l3d(3a/  +-/?,)]  ,  (19) 

with  Q  —  2J /K,  J=a  +  W2  +  8W/9,  and 
K  =4a+-2t/2Cd/(3v/3)+-(l  +-A.^/4+-X.2/2)^+-8i//9  . 

IV.  VALENCE-FORCE-FIELD  MODELS 

In  this  section  we  consider  several  models  based  on  the 
valence  force  field.  These  models  have  been  used  fre¬ 
quently  to  explain  the  impurity  bond  relaxation.5,4,8,9  We 
shall  first  derive  the  explicit  expressions  for  these  models 
and  then  connect  them  with  the  existing  results. 

A.  Model  A:  Third-shell  atoms  and  beyond  are  fixed 
at  their  pure  crystal  positions 

Let  the  bond  lengths  surrounding  the  impurity  again  be 
d\—di\—b)  and  let  the  second-shell  atoms  have  radial 
displacements  of  the  forms  (d/v' 3),  (y,y,0),  etc.  Beyond 
and  including  the  third  shell,  all  the  other  atoms  are  held 
at  their  pure-crystal  positions.  There  are  nine  different 
contributions  to  the  strain  energy  in  VFF  (see  Appendix 
B):  the  a  terms  from  the  four  bonds  surrounding  the  im¬ 
purity,  6a/(6  — 8o)2d2;  the  p  terms  among  the  six  pairs  of 
these  bonds,  ppb  —  b(,)2d2;  the  p  terms  between  the  first- 
and  second-shell  bonds,  2Pib  +  y)2d2\  the  a  terms  from 
the  second-shell  bonds,  2a(S  +  2y  )2d2;  the  P  terms  among 
the  second-shell  bonds,  2 pb2d2;  the  P  terms  between  the 
second-  and  third-shell  bonds,  jpb2d2  +  PIS  +  2y)2d2;  the 
a  terms  from  the  third-shell  bonds,  8 ay2d2;  the  P  terms 
among  the  third-shell  bonds,  4 Py2d2;  and  the  P  terms  be¬ 
tween  the  third-  and  fourth-shell  bonds,  6 Py2d2.  Thus, 
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the  excess  energy  (in  this  case  7  times  the  strain  energy) 
becomes 

A£  =  [  7O/(6-8o)2+70/<6^6o)24-70(6-|-y)2 

+  70(64- 2y  )24-  /382+  {frb+ly)1 

+2ay2+\PY2y2 .  (20) 

The  minimization  of  A£  with  respect  to  8  and  y  leads  Eq. 
(20)  to  y  =  —6/4,  and 

6=5o/[  1  +(a+  \lf}/2)/(6a,  -hPi>]  ■  (21) 

We  note  that  there  is  some  ambiguity  in  the  third  con¬ 
tribution  listed  above  for  the  P  terms  between  the  impuri¬ 
ty  and  host  bonds.  The  value  of  P  could  be  chosen  as  one 
of  these  combinations  p,  Pi,  j(p+Pi),  V PPi<  or  other 
proper  combinations.  Because  the  values  of  p  and  Pi  are 
comparable  and  p  values  are  much  smaller  than  a  (see 
Table  l),  the  results  for  6  and  A£  are  not  too  sensitive  to 
the  choice.  There  is  also  some  ambiguity  in  the  values  for 
d(0,  dy0’  for  the  “undistorted”  crystal.  The  —d2/ 3  used 
is  the  simplest  choice.  A  different  choice  will  not  affect 
the  results  for  5  at  all,  but  will  make  AE  slightly  different. 
In  fact,  model  A  was  first  used  by  Martins  and  Zunger.3 
However,  their  expression  for  6  is  different  from  Eq.  (21) 
because  they  made  different  choices  of  the  two  quantities 
just  mentioned.  Nevertheless,  Sec.  IV  will  show  that  these 
two  expressions  yield  very  similar  results.  These  ambigui¬ 
ties  do  not  occur  in  the  full  theory  in  Sec.  Ill,  where  the 
impurity-host  interactions  are  taken  into  account  natural¬ 
ly  by  the  replacement  energy  A*  [see  Eq.  (10)]. 

B.  Model  B:  Second-shell  atoms  connect 
to  a  fixed  boundary 

This  model  corresponds  to  y=0  in  model  A.  So  we 
nave 


JUJTf 

A£ =(  ja,(6  -  6o)J  +  7^/(6- 60)2  +  |a62  + -f/S62]d2 

(22) 

and 

6=6o/[l-Ma+19/?/4)/(3a/  +  70/)]  .  (23) 

This  expression  will  be  used  to  study  the  effect  of  trunca¬ 
tion. 


C.  Model  C:  Simple  spring  model 

If  all  the  £Ts  in  Eqs.  (22)  and  (23)  are  set  equal  to  zero, 
we  have  the  simple  spring  model  with 


A£  =  [ja1(8-80)2+7a62]d2 

(24) 

8=80/(14-70/0/)  . 

(25) 

The  spring  model  recently  discussed  by  Shih  et  al .*  corre¬ 
sponds  to  Eq.  (25)  with  a=at,  506/60=7. 


D.  Model  D:  VFF  with  the  continuum  connected 
to  the  second-shell  atoms 

/.  Model  D1 

In  this  case,  A E  only  contains  the  first  five  contribu¬ 
tions  listed  for  case  A  plus  the  elastic  energy  in  the  con¬ 
tinuum.  However,  the  P  terms  between  the  second-  and 
third-shell  bonds  are  modified  because  atoms  outside  R  in 
the  continuum  now  have  radial  displacements  proportion¬ 
al  to  the  inverse  of  the  square  of  the  radius.  The  result  is 


A  E 


{a,(6-60)2-i  7/3,(6-6o)2+Ta(5-f-2y)2+|£(6-fy)2-t-7^82+T/^38-l-A.Iy)2+|^(6+A.2r)J+-^75Crfr2‘ 


d2. 


(26) 


where  A.,  and  A,  arc  the  same  as  the  constants  that  appear  in  Eq.  (16).  The  corresponding  equilibrium  condition  can  be 
shown  to  be 

b  6p/!  1  4  fcr(l-2£>)+  190/4  — 0(1+ 3A.,/4  +  k2/4)Q]A3a, +P, /2)\  ,  (27) 


where 

(2a-r-/?+3A,ff/4  +  A20/2)  _ 

^  (4a4  2v led /3V' 3  404  7 A.?0+ yA.^) 

2  Model  D2 

A  comparison  between  F.q  (26>  and  the  full  perturba¬ 
tion  theorv.  Eq  » I #»>.  shows  two  maior  differences.  First, 
all  the  chemiul  term'  are  absent  in  Eq  1261  Secondly, 
the  terms 


j/3/(6-60)2d24-i/382d2 

in  Eq.  (16)  become 

70,(5-6o)2d24  7062d2 

in  Eq.  (26).  This  difference  in  the  strain  energy  will  mask 
the  true  effects  of  chemical  forces  if  6  from  Eq.  (261  is 
compared  with  FPT.  A  better  way  to  study  the  chemical 
effects  is  to  use  the  following  equation: 
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A£=[7a/(6-5)2  +  yA(8-6o)J+Ta<s+2r)2  +  T^5+y)2+7/58J  +  T^36+A.Ir)2+7/3<S+A.jy)2]d2+-^-Cr2<<J ,  £ 

(29)  _ 

ft 

which  is  Eq.  (16)  with  all  the  chemical  terms  neglected.  The  corresponding  6  becomes 

6=80/|l+[a(l-2C)  +  Jf^-^l  +  7^  +  7A.2)C]/(3ai+ft)|  ,  (30) 

' -  s 

where  the  mixing  energy  is  given  by 

A£m  =F^f—ixEAC-lryEBC)  |7 

=&N(zAAr AA  +tABrAB+EBBrBB)—6N(xeAA  +yzBB) 


with  Q  still  given  by  Eq.  (28). 

E.  Model  E:  Continuum  connected 
to  the  fint-ihell  atoms 


In  this  case,  y  =  —  3 8 v^2 )  and  A £  only  includes 
the  first  three  contributions  listed  in  model  A  plus  the 
strain  energy  of  the  continuum: 


A£  = 


ja/(8-5o)2+70f(8-8o)2 
|2 


+  10 


1- 


3vl 

8t/2 


r2+7-ctf82 


(31) 


The  relaxation  parameter  is  given  by 

8-8o/ 


31/3 

2 

1+ 

1 

\cd  + 

1  8t/2 

0 

/(3 a,  +  i0,) 

(32) 


We  note  that  the  continuum  model  used  to  estimate  the 
bond-length  relaxation  by  Baldereschi  and  Hopfield9  cor¬ 
responds  to  Eq.  (32)  without  the  0  terms,  which  yields 
6/6qcs0.4  to  0.3,  rather  than  the  proper  values  around  0.7 
to  0.8. 


=6NrAe,  (35) 

where 

Ae=e>t*  —  t(e,m+e*s)  ■  (36) 

The  mixing  entropy  AS  can  also  be  written  from  a  simple 
generalization  of  the  random  distribution.12  For  modest 
pressure,  AE  is  the  same  as  the  mixing  enthalpy  A Hm . 

Now  the  pair  interaction  energies  can  be  approximately 
related  to  the  impurity-substitution  energies  by 


A,(4  in BObsIIUju — tBB ) 

(37) 

and 

Ag(,4  in  £C)cs12(e^| — eaa  )  . 

(38) 

Thus,  Ae  of  Eq.  (36)  becomes 

AE=yr[Af(i4  in  EC)+A,(B  in  ,40] 

=  -j-[A£(,4  in  BC)+ A£(£  in  ,40]  .  (39) 


V.  ALLOY  MIXING  ENTHALPY 

The  impurity-substitution  excess  energies  A £  provide  a 
first  estimate  of  the  mixing  enthalpies  of  pseudobinary  al¬ 
loys.  Most  current  thermodynamics  theories  of  semicon¬ 
ductor  alloys  are  based  on  an  extension  of  the  binary  solu¬ 
tion  model.’2  In  this  model,  the  mixing  Helmholtz  energy 
of  an  _XC  alloy  is  defined  as 

A Fm=Ftaoy-  (xFAC + yFgc )  ,  (33) 

where  y  =  1  -x,  and  FAC  and  Fsc  are  the  respective  free 
energies  of  the  pure  AC  and  BC  compounds  at  the  same 
temperature.  Because  the  C  atoms  occupy  a  sublattice, 
the  nearest  neighbors  of  A  and  B  atoms  in  the  alloy  are 
the  C  atoms.  Thus,  the  pair  potentials  that  enter  the 
binary  solution  theory  are  now  the  second-neighbor  in¬ 
teractions  Let  NAA,  Na8,  and  SBB  be,  respectively,  the 
numbers  of  the  second-neighbor  A  A,  AB,  and  BB  pairs, 
with  corresponding  pair  interaction  energies  eAA,  e^s,  and 
fgg.  For  tetrahedral  semiconductors,  there  are  a  total  of 
6 N  second- neighbor  pairs  for  a  crystal  containing  N  unit 
cells.  Denote  the  ratios  N AA,  NAB,  and  NBB  to  6 N  as 
rAA,  rAP  =r,  and  rBB,  respectively.  Then  those  ratios  are 
related  to  the  alloy  composition  by  rAA  =x  -r/2  and 
rBB  -y  -  r  /2.  The  mixing  free  energy  has  two  terms, 

6  =  A£m  —  T  AS  ,  (34) 


Usually,  the  experimental  A/fm  is  written  as 

A//m  =x(  1  —x  )tl ,  (40) 

which  is  equivalent  to  assuming  a  random  distribution, 
i.e.,  r=2x(l—  x).  Using  this  expression  for  r  and  com¬ 
paring  Eqs.  (40)  and  (33),  we  see  that  the  mixing  enthalpy 
parameter  fl  is  given  by 

n  =  2[A£(,4  in£C)  +  A£(£  in  ,40]  .  (41) 

This  connection  provides  a  further  check  of  the  theory. 

VI.  NUMERICAL  RESULTS  AND  DISCUSSION 
A.  Chemical  terms 

Table  II  lists  8q=  l—d/Zd,  8=1  —dx/d,  the  excess  en¬ 
ergy  (per  bond)  A £  for  the  full  theory  and  its  correspond¬ 
ing  VFF  model  D2,  and  the  terms  derived  from  the 
metallization  energies,  [Eq.  (19)],  £<*,  H,  A Ep,  and 
A£ch  [Eq.  (17)].  The  appropriate  derivatives  /,  g.  A,. . . 
[see  Eq.  (9>J  are  computed  using  the  atomic-term  values 
that  we  have  generated  from  impurity-level'5  and 
structural  studies/ 

For  substitutions  involving  the  cation  pair  (Ga,Al),  Fcb 
has  the  same  sign  as  Sq,  which  means  that  £ch  prevents 
relaxation  and  thus  tends  to  increase  the  strain  energy. 
The  chemical  forces  H  are  also  significant.  As  a  result, 
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TABLE  II.  Comparison  between  the  full  theory  and  the  corresponding  VFF  model  D2  to  study  the  effects  of  chemical  terms.  All 
AiTs  are  in  units  of  kcal/moleband 
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all  six  cases  involving  this  pair  have  nearly  equal  dt  and 
d2 ,  i.c.,  the  small  bond-length  differences  are  made  even 
smaller  in  the  alloy.  The  excess  energies  all  become  nega¬ 
tive,  mainly  because  A Ep  is  negative.  For  the  systems  in¬ 
volving  the  (Ga.ln)  and  (ln,Al)  pairs,  has  the  opposite 
sign  from  So.  so  8^  and  8q  have  the  same  sign.  The  chemi¬ 
cal  force  favors  bond  distortion.  However,  because  H  is 
positive  and  introduces  an  increase  in  the  denominator  of 
Eq.  (18),  most  of  the  effect  of  So  is  cancelled.  For  cases 
involving  (Ga,ln),  the  polarity  contributions  h£p  are  all 
negative.  The  Fch  Ad  term  is  negative,  but  //(Ad)2  is  pos¬ 
itive,  so  they  cancel  to  a  certain  degree  and  leave  A E 
lowered  primarily  because  of  A Ep.  While  A Ep  is  still  neg¬ 
ative  for  the  (In,Al)  substitutions,  its  magnitude  is  reduced 
considerably.  The  other  chemical  energies 
*ch  Ad  +//(Ad)2  can  be  as  large  as  A Ep,  but  the  overall 
reductions  of  A E  are  only  a  fraction  of  those  for  the 
(Ga,In)  cases.  For  the  several  II- VI  systems  studied,  both 
^ch  and  H  are  small  and  the  net  changes  in  6  have  the 
same  sign  as  80.  However,  because  8q  is  small  in  the 
(Cd,Hg)  substitutions,  ^cb  actually  causes  a  reversal  in 
which  the  short  bond  length  becomes  shorter  and  the 
longer  one  becomes  longer.  This  is  the  only  exceptional 
case  of  this  type  found  for  all  the  systems  studied.  The 
change  of  A E  due  to  chemical  terms  in  the  (Hg,Zn)  substi¬ 
tution  is  also  peculiar — it  increases  mainly  because  A Ep  is 
positive. 

Next,  we  examine  the  anion  substitutions.  For  the 
groups  involving  the  (P,As)  pair,  the  chemical  shifts  are 
all  small,  but  the  trend  is  less  toward  relaxation  and  larger 
A E.  The  groups  involving  (As,Sb)  and  (P,Sb)  pairs  behave 
very  similarly:  Fch  are  significant  and  are  opposing  relax¬ 
ations,  i.c.,  6q  and  bu  have  opposite  signs.  At  the  same 
time,  the  H  values  are  several  times  smaller  than  those  for 
the  corresponding  III  V  cation  substitution  case.  Thus, 
most  of  translates  into  a  real  reduction  of  the  ratio 
6/6,,  and  consequently  introduces  extra  strain  energy.  Al¬ 
though  the  A Ep  energies  are  significant  and  negative, 
F  h  Ad  are  positive  and  the  net  AEch  can  be  either  positive 
or  negative  However,  (he  induced-strain  energy  due  to 
reduction  ot  the  6/b,,  makes  all  A E  positive  for  these  two 
groups  of  systems  For  ll-Vl  systems,  all  the  chemical  ef¬ 
fect  again  are  small,  buf  (he  net  chemical  changes  on  A E 
are  slightly  positive 

The  above  discussion  can  be  summarized  in  Fig.  2, 
where  the  excess  energies  A F  calculated  from  the  full  per¬ 
turbation  theory  and  Model  D2  are  plotted  normalized  to 
the  results  of  the  simple  spring  model  of  Shih,  Spicer, 
Harrison  and  Sher  (SSHS)  (Ref.  8),  i.c.,  Eq.  (24)  with 
ctj  -a,  so  A E/\  ,  nd '  I  b p.  The  calculated  A E  rises  fas¬ 
ter  for  than  for  mainly  because  a/a,=  1.  In 

fact,  if  the  relation  lu  a /a,  —  d//d'A  with  S  of  order  of 
3  to  5  is  used  in  F.q  (24>,  we  obtain  a  percentage  correc¬ 
tion  of  S’fv'4  to  the  SSHS  results,  which  explains  the 
skewed  behavior  of  the  curve  It  is  also  clear  that  A E 
rises  faster  than  b  for  larger  However,  the  zeroth- 
order  iheoiv  of  SSHS  is  clearly  an  excellent  representation 
of  the  global  features  of  \F.  The  results  from  model  D2 
are  dost  r  to  (he  parabolic  form  than  those  from  FPT. 
The  figur--  clearly  shows  the  general  (rends;  the  chemical 
terms  ,  ause  negative  shifts  tn  Af.'  for  cation  substitutions 
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FIG.  2.  The  excess  energies  A E  over  jad2  calculated  from 
the  full  perturbation  theory  (FPT)  and  its  corresponding 
valence-force-field  model  D2.  The  solid  curve  corresponds  to 
Eq.  (24)  with  a— a,. 


and  positive  shifts  for  anion  impurities.  It  is  also  clear 
that  the  chemical  shifts  can  be  very  large.  These  effects 
will  have  important  consequences  on  the  alloy  mixing 
enthalpies  to  be  discussed  later. 


B.  Impurity  bond  length 

Table  III  lists  the  impurity  bond  lengths  dt  calculated 
from  different  models,  while  a  comparison  of  theory  and 
the  available  experimental  data1,14’”  is  presented  in  Table 
IV.  The  actual  size  of  changes  in  dt  induced  by  the 
chemical  terms  can  be  seen  by  comparing  model  D2  with 
the  full  theory.  Except  for  the  systems  involving  the  sub¬ 
stitution  of  (Ga,As)sand  (P,Sb)  pairs  (where  changes  range 
from  0.01  to  0.03  A),  all  the  chemically  induced  changes 
are  less  than  0.01  A.  Comparison  among  models  A,  B, 
and  C  shows  that,  while  extending  the  boundary  helps  the 
relaxation  (compare  model  B  to  A),  i.e.,  5/6o  is  closer  to 
1 ,  the  inclusion  of  the  bond-bending  forces  (the  f)  terms) 
(compare  model  B  with  C)  prevents  it.  The  simple  spring 
model  (model  C),  which  contains  neither  of  these  terms, 
evidently  represents  a  delicate  cancellation  of  these  effects 
and  predicts  results  close  to  those  of  the  full  perturbation 
theory  and  experiment.  Although  the  d\  values  of  model 
C  are  often  very  close  to  those  of  model  A,  there  are  cases 
{e.g.,  Ga(P,Sb)]  in  which  model  C  can  differ  from  model 
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A  by  0.02S  A.  Model  A  produces  about  the  same  d, 
values  as  model.Dl,  where  the  maximum  difference  in  d , 
is  only  0.015  A.  Martins  and  Zunger*  used  the  same 
model  as  model  A;  however,  their  analytic  expression  for 
5  is  different  from  that  given  by  Eq.  (21).  Nevertheless, 
numerical  results  indicate  that  these  two  calculations 
agree  to  0.01  A.  The  slightly  different  forms  of  strain  en¬ 
ergies  used  in  models  D2  and  D1  only  introduce  a  small 
change  in  dlf  with  the  largest  difference  being  less  than 
0.01  A.  The  first-shell  continuum  model  (model  E)  allows 
too  little  relaxation,  so  while  the  other  models  produce  a 
ratio  5/Sq  ranging  from  0.6  to  0.8,  model  E  only  ranges 
from  0.4  to  0.6.  The  reason  that  the  fixed  boundary  in 
model  A  works  is  that  the  effective  shear  coefficient  C 
(see  Table  I)  characterizing  the  strain  energy  in  the  elastic 
continuum  is  large.  However,  model  B  is  too  rigid  and 
does  not  provide  enough  buffer  between  the  impurity 
bond  and  the  fixed  boundary. 

The  comparison  of  the  theoretical  results  with  the 
available  experimental  data  in  Table  IV  indicates  that 
models  B  and  E  are  the  least  accurate.  Models  A,  01, 
that  of  Martins  and  Zunger  (MZ),  and  the  full  theory  are 
comparable  in  that  all  have  an  average  absolute  deviation 
of  0.012  A,  which  is  about  the  experimental  uncertainty  in 
extended  x-ray-absorption  fine  structure  (EXAFS).  The 
agreement  between  theory  and  experiment,  however,  is  not 
uniform.  The  most  surprising  result  in  Table  IV  is  that 
the  simple  spring  model  (model  C)  and  its  cruder  version 
used  by  Shih  et  al .'  (a-a,  so  6/6o=0.75,  labeled  as 
SSHS)  have  the  smallest  variance  in  du  about  0.006  A. 
We  know  this  docs  not  imply  that  the  simple  spring 
model  represents  the  real  picture  of  bond-length  relaxa¬ 
tion.  For  example,  if  we  let  all  the  shear  coefficients  be  0, 
i.e.,  /9=C =0  in  our  model,  then  as  the  range  of  the  boun¬ 
dary  is  gradually  extended,  the  local  bond  length  will 
eventually  relax  to  the  impurity  bond  length  d\=dj,  or 
8  =  6o.  This  can  be  seen  in  model  A  from  Eq.  (21),  where 
8  reduces  to  8o/(l-fa/6a/),  rather  than  8o/(l4-a/3a,) 
as  predicted  by  model  C,  and  in  model  D  from  Eqs.  (27) 
and  (30),  8  becomes  80,  if  the  continuum  is  taken  to  be 
shearless.  Considering  that  various  effects  are  included 
that  may  mask  the  absolute  accuracy  of  dt  predictions 
(e.g.,  while  low-temperature  bond  lengths  are  used  in  the 
calculation,  the  room-temperature  values  of  elastic  con¬ 
stants  are  adopted),  the  agreement  of  various  models  with 
experiments  in  Table  IV  should  be  regarded  as  excellent 
There  are,  however,  many  other  impurity  systems  in 
which  the  simple-spring-model  predictions  differ  consid¬ 
erably  from  other  models,  as  is  shown  in  Fig.  3,  where  8  is 
plotted  against  8q  for  the  full  theory.  Those  points  that 
deviate  significantly  from  the  0.75-slope  line  are  the  sys¬ 
tems  with  (As,Sb)  or  (P,Sb)  substitutions.  Additional  EX¬ 
AFS  measurements  on  these  systems  are  needed  to  test 
these  predictions. 

C.  Mixing  enthalpies 

Table  V  lists  the  mixing  enthalpy  parameters  fi  (in 
kcal/mole)  for  a  number  of  alloys  estimated  from  Eq.  (41) 
for  all  the  models  considered,  along  with  other  theoreti¬ 
cal3, 16_ll,  and  experimental  values.16  ”  As  already  dis- 
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TABLE  V.  Mixing  enthalpy  parameter  0  (in  kcal/mole)  estimated  from  the  full  perturbation  theory  and  several  valence  force 
models  discussed  in  Sec.  Ill,  and  comparison  with  experiments  and  other  theories. 


A 

B 

C 

Dl 

D2 

E 

FPT 

MZ* 

DLb 

FMC 

W" 

Expt.' 

(Ga.Al)P 

0.00 

0.01 

0.00 

0.01 

0.0 1 

0.01 

-0.05 

(Oa.ADAs 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

-0.07 

0.02 

0.02 

0.03 

0.11 

0.0 

(Ga^U)Sb 

0.02 

0.03 

0.02 

0.02 

0.02 

0.03 

-0.15 

0.02 

0.02 

0.03 

0.0 

(Ga,In)P 

3.76 

4.79 

3.0 

3.29 

3.39 

5.24 

2.54 

4.56 

3.63 

2.94 

3.25,3.5 

(Ga.InlAs 

2.97 

3.76 

2.36 

2.61 

2.69 

4.14 

1.60 

2.49 

2.81 

2.42 

1.25 

1.65,  2.0,  3.0 

(Ga.In)Sb 

2.22 

2.83 

1.77 

1.95 

2.00 

3.09 

0.81 

2.53 

1.85 

1.83 

1.47,  1.9 

(In.Al)P 

3.24 

4  22 

2.77 

2.78 

2.87 

4.60 

2.55 

(In.ADAs 

2.86 

3.65 

2.32 

2.49 

2.56 

3.93 

2.17 

3.60 

2.81 

2.37 

2.5 

(In.AllSb 

1.81 

2.33 

1.49 

1.57 

1.61 

2.50 

1.36 

2.06 

1.46 

1.45 

0.6 

<Cd,Zn)Te 

1.80 

2.43 

1.73 

1.43 

1.49 

2.45 

1.24 

2.12 

1.97 

1.63 

1.34 

(Hg,Cd)Te 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

-0.07 

0.7,  1.4 

(Hg,Zn)Te 

1.63 

2.20 

1.56 

1.30 

1.36 

2.23 

1.50 

1.91 

1.81 

1.48 

3.0 

AKPhAs) 

0.81 

1.03 

0.65 

0.71 

0.73 

1.14 

0.76 

Ga(P,As) 

0.95 

1.18 

0.70 

0.86 

0.87 

1.32 

0.94 

1.15 

0.98 

0.66 

0.12 

0.4,  1.0 

In(P,As) 

0.60 

0.78 

0.52 

0.51 

0.53 

0.84 

0.57 

0.72 

0.58 

0.52 

0.4 

AKAs.Sb) 

4.31 

5.45 

3.38 

3.80 

3.88 

5.92 

4.09 

Ga(As.Sb) 

3.77 

4.69 

2.81 

3.40 

3.46 

5.22 

3.67 

4.58 

3.35 

2.76 

4.0,  4.5 

IntAs.Sb) 

2.61 

3.39 

2.23 

2.24 

2.31 

3.67 

2.52 

2.89 

2.29 

2.17 

6.65 

2.25,  2.9 

AKP.Sb) 

8.60 

10.99 

6.99 

7.54 

7.73 

12.00 

8.32 

Ga(P^b) 

8  54 

10.61 

6.36 

7.72 

7.88 

11.66 

8.66 

In(P.Sb) 

5.87 

7.64 

5.08 

4.99 

5.15 

8.04 

5.76 

Zn(S.Se) 

1.04 

1.39 

0.98 

0.85 

0.90 

1.49 

0.90 

Zn(S,Sc> 

1.04 

1.39 

0.98 

0.85 

0.90 

1.49 

0.90 

ZotSe.Te) 

2.47 

3.27 

2.23 

2.09 

2.16 

3.63 

2.26 

2.91 

3.11 

2.12 

3.12 

1.55 

Zn(S.Te) 

7  02 

9.34 

6.45 

5.80 

6.02 

9.72 

6.20 

•Reference  3,  column  A  of  Table  III. 
‘Reference  16. 

•Reference  18. 

‘Reference  17. 

•References  16  and  19. 


cussed,  the  chemical  terms  reduce  the  excess  energies  in 
the  cation  impurities  and  increase  them  for  anion  impuri¬ 
ties.  The  corresponding  changes  in  ft  are  the  differences 
between  the  FPT  and  model  D2.  We  note  that  the  reduc¬ 
tions  of  ft  for  the  (Ga,In)  alloys  are  very  large  ( >  1 
kcal/mole)  and  also  significant  for  (In,Al)  alloys.  Howev¬ 
er,  the  increases  in  fl  for  the  anion  s  ibstitutional  alloys 
are  not  as  large.  Also,  the  f 1  differences  between  models 
D1  and  D2  are  less  than  10%.  Model  A  produces  ft 
values  about  20%  larger  than  model  Dl,  model  B  in  turn 
is  20%  higher  than  model  A,  and  model  E  is  10%  higher 
than  model  B  The  fl  values  in  the  simple  spring  model 
(model  C)  are  seen  to  be  about  the  same  as  model  Dl,  al¬ 
though  the  differences  among  systems  can  be  positive  or 
negative.  Although  MZ  used  the  same  strain  model  as 
model  A,  their  fl  values  do  not  agree  with  our  model  A 
values  because  their  way  of  estimating  fl  is  different.  In 


fact,  the  values  of  MZ  are  closer  to  model  B  than  to  A. 

To  distinguish  the  quantitative  nature  of  different 
theoretical  models,  we  note  that  there  are  also  important 
factors  that  may  mask  the  comparison  between  theory 
and  experiment  for  ft.  One  important  factor  is  that  the 
mixing  enthalpies  extracted  from  phase-diagram  analysis 
are  sensitive  to  sample  and  experimental  conditions. 
These  AH  contain  contributions  from  various  nonideai 

m 

structures  such  as  vacancies,  impurities,  dislocations, 
grain  boundaries,  and  surface  conditions,  in  addition  to 
the  ideal  AHm  considered  here  for  solid  solutions.  Thus, 
our  theoretical  AHm  should  represent  a  lower  bound. 
Another  complication  comes  from  the  version  of  the 
theory  of  solid  solution  adopted.  The  theory  used  for 
analysis  so  far  assumes  a  regular  solid  solution  with 
second-neighbor  pair  interactions  as  was  outlined  in  Sec 
V.  Recent  experimental*'0  22  and  theoretical"’1  studies 
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FIG.  3.  Calculated  bond- relaxation  parameter  8  from  FPT  and  D2  as  a  function  of  So-  The  6=  780  curve  corresponds  to  the 
theory  of  Shih,  Spicer,  Harrison,  and  Sher  (Ref.  8). 


have  suggested  the  possibility  of  compositional  clustering 
or  long-range  correlations  in  alloys.  Extending  the  theory 
to  include  such  effects  will  alter  the  simple  results  for 
£Jim  from  Eq.  (41).  Moreover,  there  is  evidence  from  the 


i«i  theory 


composition  variation  of  the  alloy  hardness24  and  from 
the  optical-phonon  frequencies25  that  the  shear  coeffi¬ 
cients  of  alloys  increase  near  the  center  of  the  composition 
range.  This  will  cause  the  effective  continuum  shear  coef- 


Ibi  ixnawHNt 


HO  4  Plot  of  1  /  .  as  a  function  ft..,  /bm  for  J,  obtained  from  (a)  FPT  and  (b)  the  experimental  II  in  Table  V.  In  la),  the  cir 

.  It-',  ire  foi  anion  alloys,  the  squares  for  cation  substitution.  The  solid  tines  in  both  parts  correspond  to  the  simple  theory  discussed  in 
Appendix  <  I  he  dashed  tines  at  l,  / 1  -  0  separate  the  miscible  from  immiscible  groups. 
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TABLE  VI.  Comparison  of  the  critical  temperature  T,  of  mining  and  melting  temperatures  of  the 
constituents  7  ,  and  T2,  in  the  order  of  their  appearance  in  the  parentheses.  Also  tabulated  are  the 
averaged  absolute  values  of  I  bo  i  and  the  ratio  |  bo  |  /b„ . 

i  bo  I  T, 


System 

{%) 

(K) 

Tc/Tt 

tjt2 

|6ol/6 

(Al.Ga)Sb 

0.6 

0 

0 

0 

0.067 

(AI.Ga)As 

0.1 

0 

0 

0 

0.009 

(Al,Ga)P 

0.3 

0 

0 

0 

0.029 

(Al.In)Sb 

5.5 

342 

0.25 

0.42 

0.679 

(Ga.ln)Sb 

5.8 

204 

0.21 

0.25 

0.716 

(AI,In)As 

6.8 

547 

0.29 

0.45 

0.687 

(Ga.In)As 

6.9 

403 

0.23 

0.34 

0.697 

(Al,In)P 

7.1 

642 

0.36 

0.48 

0.732 

(Ga.In)P 

7.4 

639 

0.43 

0.48 

0.763 

In(P,As) 

3.2 

144 

0.11 

0.12 

0.330 

Ga(PAs) 

3.7 

236 

0.14 

0.15 

0.352 

Ga(P,As) 

3.5 

191 

0.10 

0.11 

0.307 

ln(As,Sb) 

6.8 

635 

0.52 

0.79 

0.840 

Ga(As,Sb) 

7.6 

924 

0.53 

0.94 

0.844 

AKAs.Sb) 

8.1 

1030 

0.56 

0.78 

0.810 

In(P,Sb) 

9.9 

1450 

1.08 

1.82 

1.222 

Ga(P.Sb) 

11.3 

2180 

1.25 

2.21 

1.256 

Al(P,Sb) 

11.6 

2095 

1.19 

1.58 

1.116 

(Cd.Hg)Te 

0.3 

0 

0 

0 

0.033 

(Zn.Hg)Te 

6.0 

377 

0.25 

0.40 

0.659 

(Zn,Cd)Te 

6.2 

312 

0.21 

0.23 

0.564 

Zn(S.Se) 

4.7 

226 

0.11 

0.13 

0.362 

Zn(Se.Te) 

7.2 

569 

0.32 

0.38 

0.615 

Zn(S.Te) 

11.9 

1560 

0.74 

1.03 

1.017 

ficient  C  in  Eq  (14)  to  be  composition  dependent,  which 
will  cause  II  to  increase.  Despite  these  uncertainties,  use¬ 
ful  comparisons  across  the  board  in  Table  V  can  still  be 
made. 

Based  on  the  above  considerations,  we  can  conclude 
that  models  B,  E,  and  MZ  predict  11  values  that  are  too 
high.  We  should  emphasize  that  all  the  LI  numbers  for 
models  from  A  through  MZ  are  directly  calculated 
without  any  adjustable  parameters.  The  fact  that  models 
A,  C,  and  D1.D2,  and  the  FPT  agree  with  the  experiment 
as  well  as  or  even  slightly  better  than  the  one-parameter 
theories,  the  delta  lattice  (DL)  model1'1  and  the  model  of 
Fedder  and  Muller18  (FM),  is  already  quite  an  accomplish¬ 
ment.  The  few  numbers  taken  from  Van  Vechten’s  calcu¬ 
lations17  (VV)  indicate  that  the  dielectric  model  predicts 
results  at  larger  variance  with  experiments.  There  are  two 
important  implications  about  the  FPT  that  can  be  drawn 
from  Table  III.  First,  the  theory  predicts  a  small  but  neg¬ 
ative  11  value  for  several  alloys.  This  not  only  means  that 
there  is  no  miscibility  gap  in  these  alloys  but  also  implies 
a  tendency  toward  ordering,  in  which  the  substitutional 
atoms  tend  to  be  surrounded  by  different  second-neighbor 
species  For  stoichiometric  compositions,  this  implies  a 


tendency  toward  compound  formation.  Secondly,  the 
FPT  tends  to  predict  smaller  1 I  values  than  observed  ex¬ 
perimentally,  which  should  be  expected  according  to  our 
discussion.  To  the  extent  that  the  FPT  predicts  the 
correct  A Hm  values  for  an  ideal  solution,  the  difference 
A7/”p'  —  A//m  may  be  attributed  to  imperfect  conditions 
and  deviations  from  the  ideal  solution  theory. 

Finally,  the  calculated  11  values  in  Table  V  provide 
some  guidance  in  separating  the  completely  miscible  al¬ 
loys  from  immiscible  ones.36  27  In  a  true  random  alloy, 
the  criteria38  for  alloy  mixing  at  a  temperature  T  is  that 
T^TC,  where  the  critical  temperature  Tc  is  given  by 
il/(2Rg),  with  Rg  being  the  universal  gas  constant.30  For 
an  AxBt_xC  alloy  to  be  miscible  throughout  the  whole 
concentration  range,  the  requirement  is  that  both  the 
melting  temperatures  T\  and  T2  of  the  pure  AC  and  BC 
compounds  be  greater  than  7’c.  Table  VI  lists  the  values 
of  Fc  associated  with  (he  il  values  in  the  FPT,  the  ratios 
7r/7|  and  7r/7\,  and  the  average  absolute  values  of  6,, 
for  the  alloys  considered  in  Table  V  In  Table  VI,  T,  is 
set  equal  to  zero  if  fl  is  negative  and  7\  is  chosen  to  be 
the  lower  value  of  the  two  melting  temperatures,  so  the 
criterion  for  not  having  a  miscibility  gap  is  Tr/T:<\. 
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~*i  t-.npi:ua!  iule‘f  stating  that  tins  will  be  satis- 
•  ia ttn.e  mismatch  i>,,  '  between  the  two  alloy 
mis  is  less  than  7 .5%  However,  we  find  that  (see 
'  C  -i  more  precise  rule  is  6n  >  fim ,  where 
,,  'd  V  „  is  the  ratio  of  the  rms  bond-length 
ie  flues  nation  :o  the  bond  length  at  the  melting 
t:  >■  !  The  values  of  7  „  for  the  compounds  in- 
tr  '  the  associate  tm  sallies  estimated  from  Eqs. 
iC  !  are  tabulated  in  Table  I.  The  model  used  in 
•:  C  telris  V,  ;  • '  %'bm  )2.  This  suggests  that 

1  m  plot  7  •'  T,  as  a  function  of  iS0|  /bm, 

i  tn  Fr:  4  lor  the  Tc  calculated  from  FPT. 

*  *  ■; m i It- r  the  AT  versus  8„  curve  in  Fig.  2  be- 

fT.  ft  i  .  proportional  to  the  sum  of  the  A E 

o'  i  •;  s-.tueiits  [see  Eq.  (41)].  However,  if 
■i !  ;  iigau.si  6,. .  alone,  the  FI’T  points  are 
.<»> .red,  u.n  l  those  of  SSHS  would  not  even 
-•ms-il.  simple  quadiatic  form  because  the  lower 
(e  npvnmire  7  ;  is  not  a  smooth  function  of  ]  S0  I  . 
\.  vuggests  that  iy  bm  <  !  is  a  better  criterion 
tup, in  t  also  clearly  shows  the  chem- 
j  .it, on  substitution  alloy  points  lie 

,  ■  i i . •  ..  u-\  .uui  all  the  anion-substitution  alloys 
I  •  :  it: or  n»  c. iwsve  the  curve,  corresponding 
•i  ;  ’  .e  ••".iff.  in  IE  due  to  the  chemical 

,.  i  .  j.  ..  bused,  on  the  SSHS  model  is  an 
.  1  -  presentation.  From  the  figure,  we 

i.M  ;  hr"  at-  v,  should  have  miscibility  gaps  and 
h  n'!o\>  nit  pi  filleted  to  be  miscible,  although  on 
(crime,  because  the  actual  mixing  enthalpies  are 
i at.  these  i-itvi  calculated  values.  The  figure  also 
ir  S  I  ■'  has  a  miscibility  gap  but  a  smaller 
•'  ■  '  ,  -.hr  it’.Sb1  alloys,  despite  the  faci 

•  niut  i*  uigei.  All  these  predictions  arc 
;■  : he  r.auaKe  experimental  evidence. 

i i  St'MMARN  AN1>  CONClft  SION 

a  -i  ■  ole  to.-uv  of  defect  substitution  en- 
1  i"  •  ,i  iu,  sulrstnution  eneigy  is  compactly 
■e  u  energy  and  a  distortion 
*  •  ,  •>  '  ■•••.  i  T  | see  Fq  ‘ 2*|.  However,  a 

i  s  u-q  litres  an  improvt- 
ts  of  Harrison’s  bonding 

■  :  ■  the  ii  ■  .  'iistants.  f  he  most  m- 

■  i ' ' ,  ■  •  i  <  ‘  fh.c  '  cry  p.usented  ;u  this  paper 

u  ,  ;  f  v.  i., h  enables  us  to  absorb  the 

'  •  -  •  ‘ff  int<  the  calculation  and, 

-'■■■■  the  chi  rr’-al  effects.  The  ori- 
:  ■  o'  n  rnpunty  f-ond  relaxation 

•  ! nr  nanisms  [see  F.q  (17)];  a 

ii  .oner  helps  or  htndcis  lattice  re- 

». '  ,  ,i  h.i.s  the  same  oi  oppo- 

"fcth  difference  d-d,  between 
.  hen. i  a'  energy  that  depends  on 
tviwcrn  the  iinpuritv  and 
■•i  r * ! ,.  ,ic  i  !a> t i,  force  constant  // 
t,-.  estiiii'.  the  lattice  from 
ht  lo  .  'n'iuid.incs  between  the 
,  m  ir-pir  "  and  'he  rest  of  the  elastic 
"i-  .  \  hi  h,-  %  alone"  force  field 


are  derived  and  their  results  are  compared  with  the  full 
perturbation  theory  and  available  experimental  data.  Wc 
found  at  least  five  models,  including  the  FPT,  that  pro¬ 
duce  the  correct  impurity  bond  lengths  with  variances  for 
the  compounds  studied  about  equal  to  the  ex  per  mental 
uncertainties  in  EXAFS1  4  14  (-0.01  A).  However,  some 
models  are  oversimplified  and  will  certainly  not  predict 
other  properties  equally  well.  However,  more  experimen¬ 
tal  lattice  constant  measurements  to  further  test  the 
theory,  particularly  on  (As.Sb)  and  (P,Sb)  substitution  sys¬ 
tems  for  which  there  are  larger  differences  between  dif¬ 
ferent  models,  are  needed.  It  would  also  be  instructive  to 
see  if  the  predicted  reversal  for  (Hg.CdlTe  is  found. 

The  excess  energies  of  impurity  substitution  are  also 
shown  to  provide  good  estimates  of  the  mixing  enthalpies 
fl  of  pseudobinary’  alloys.  The  chemical  shifts  are  found 
to  have  a  negative  net  contribution  to  f l  for  most  cation 
substitutions,  but  positive  contributions  for  anion  substi¬ 
tutions.  The  chemical  reduction  of  ft  in  (Ga,ln)  alloys  is 
larger  than  1  kcal/mole  (30—100%).  Several  VFF 
models  and  the  full  perturbation  theory  produce  results 
for  ft  that  are  as  good  as  the  best  theories  with  one  adjust¬ 
able  parameter.  However,  the  full  theory  tends  to  yield 
answers  on  the  low  side  of  the  experimental  values,  which 
we  argue  is  as  it  should  be  because  there  are  nonideal 
structures  that  also  contribute  to  ft.  The  calculated  ft 
values  and  the  melting  temperatures  are  used  to  predict 
the  existence  of  alloy  miscibility  gaps,  and  the  results 
correlate  well  with  experiments. 

Finally,  we  wish  to  comment  on  the  accuracy  of  the 
theories  that  are  connected  to  the  present  model.  The  per¬ 
turbation  theory  has  already  been  stretched  beyond  its  ex¬ 
pected  region  of  validity  and  predicts  d ,  to  within  experi¬ 
mental  uncertainties  (  —  0.01  A)  even  for  cases  with  large 
bond-leng'h  differences  (8o  — 0. 1 ).  The  accuracy  can  only 
be  improved  if  the  full  nonperturbation  theory  outlined  in 
Sec  II  is  used.  This  calculation  is  needed  for  the  strong 
substitution  cases  that  were  not  considered  in  this  paper; 
examples  are  (B,Ga),  (B.In),  iN.P),  (N.As),  (N.Sb),  (O.S), 
(O.Se),  and  (O.Te)  substitutions.  Although  we  believe  that 
for  the  properties  treated,  the  model  with  a  continuum  at¬ 
tached  to  the  second  shell  is  as  accurate  as  the  perturba¬ 
tion  theon  u'-ed.  it  remains  to  be  seen  if  this  is  true  for 
other  properties,  especially  strain  coefficients.  Finally,  the 
present  theory  has  been  extended  to  study  alloys'''’  by 
embedding  clusters  in  an  effective  medium.  This  enables 
us  tc  study  the  bond  length  and  energy  variations 
throughout  the  whole  concentration  range.  However,  a 
quantitative  calculation  still  awaits  an  improvement  of  the 
accuracy  of  Hamson's  theory  A  similar  procedure  is 
also  being  extended  to  a  study  of  the  alloy  electronic 
structure  for  which  a  cluster  CPA  (coherent-potential  ap¬ 
proximation)  involving  both  potential  and  structural  disor¬ 
der'1  will  be  used 
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APPENDIX  A:  ELASTIC  ENERGY  IN  CONTINUUM 

In  Sec.  Ill  the  elastic  energy  outside  a  sphere  of  radius 
R  centered  at  the  impurity  is  assumed  to  be  a  continuum 
with  a  radial  displacement  u  cct/r2.  If  the  displacement 
at  R  is  «0.  the*1  u(r)  =  u0(R2/r2J?.  The  energy  density  in 
the  continuum  is  given  by 

6e(r)  =  jCi^e^  +  +  e^ ) 

+  C| ~^~ejueyy  ~^^xxea  3 

+  iC44(e^r  +  )  , 


where 

3u.  i  < 

e„  =  — — =  R2u0(r2-3x2)/r5  , 
dx 

dux  d u.  . 

'-=17+af=-6/jWr . 

Thus,  the  total  elastic  energy  in  the  continuum  is 
Ajf'-  f~6e(r)d3r 

=4irRul(j;Cu  —  t^12  +  TC44) 

=  CRu20  , 

where  the  effective  shear  coefficient  is  given  by 

C  =  jrU.6C,,  — I.6C12+4.8C44)  • 


APPENDIX  B:  DISTORTION  ENERGY 

In  this  Appendix  we  count  the  detailed  contributions  of 
the  bond-stretching  terms  A(r,  r, )  and  "bond-bending” 
terms  A(r(-r t)  for  iy=j  in  VFF  [Eq.  (13)J  that  enter  Eq. 
(16)  in  FPT  and  in  the  VFF  models  in  Sec.  III. 

t.  a  and  0  terms  from  the  first-shell  bonds 

The  four  bond  vectors  pointing  away  from  the  central 
impurity  according  to  Fig.  1  are 

r,  =(  1—6,  1-6.  l-eW/v^  , 

r2  =  (  1  — 6,  -1-6,  -\-b)dVl . 

Thus,  Afr(T|)=  — 26dJ  and  A( ri *r2)  =  —  y6d2.  The  a 
terms  contribute 

4x  3a;(  —  2bd2)2  /%d2  =  barb2d2  , 

and  the  0  terms  contribute 

6x30,(  -  \bd)2/Sd2=P,b2d2  . 

If  an  A  atom  is  replaced  by  a  B  atom,  as  was  done  in 
FPT.  the  a,  and  0,  are  replaced  by  a  and  0,  respectively. 


2.  a  terms  from  the  secood-thell  bonds,  0  terms 
between  the  first-  and  second-shell  bonds 
and  among  the  second-shell  bonds 

For  these  terms  we  need  to  consider  the  four  bond  vec¬ 
tors  pointing  away  from  C.  They  are 

rj  =  (  —  1  +6,  -  1+6,  —  l+6W/v/3  , 

r2  =  (-l+6, 1+6+y,  1  +6+yk//v/3  , 

r3  =  (l  +  5+y, -1+6,  l  +  y+bkf/vl, . . .  . 

Then  A(r2T2)  =  y(6+2yW2,  A(rtT2)  =  -  y(6+ yW2,  and 
A(r2Tj)=  |8d2.  Thus,  the  a  terms  from  the  second-shell 
bonds  become 

4x  3  X3a[  }(6+2y  )d2]2/802  =  2a(6+26)2d2 , 
the  0  term  between  the  first-  and  second-shell  bonds  are 
4x3x30tT(6  +  yW2]2/8d2  =  2(y  +  6)20d2 . 
and  the  0  terms  among  the  second-shell  bonds  are 
4x3x30(y602)2/8d2  =  2620tf2  . 


3.  a  terms  for  the  third-shell  bonds,  0  terms 
between  the  second-  and  third-shell  bonds 
and  among  the  third-shell  bonds  adjacent 
to  the  second-shell  atoms 

For  these  terms  we  need  to  consider  the  bond  vectors 
pointing  away  from  B  in  Fig.  1.  They  are 

r2  =  (  1  —6,  —  1—  y  —  6,  -I— y— 8 )d/Vr3  , 

ij=(l  +  y',  1+3 y’-y,  1  +3y'-yW/t/3  , 

r3  =  (-l-y",  l+3y"-y,  -l+y"-yW/vl  , 

and 

r4=(  —  1  —  y",  — 1+y"— y,  l  +  3y"— ykft/3  . 

Thus,  we  have 

A(r,-r2)=  -  ~(3b  +  5y')  , 

A(r2Tj)=  ~(6  +  2y  —  5y")  , 

A(r,-r«)=  -  \d2y"  , 

A(r,Tj)=  —  2y  —  y'  +  3y")  , 

A(r3T3)=y(6y"), 


and 

A(r,-r,)=~(14y'-4y)  . 

For  model  A,  y'  =  r"=0,  so  the  a  terms  of  the  third-shell 
bonds  become 

4x 3 X 3a[2( Ar3-r3)2  +  (Ar|-ri )2]/8d2  =  icty2d2  , 


A 
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the  0  terms  between  the  second-  and  third-shell  bonds  are 

4x3x3/31(Ar,T2)2  +  2(Ar2Tj)2]/8d2 

=  j0d2[52  +  {(5+2y)2]=  y^82d2+^5-f-2y)Jd2  , 

and  the  &  terms  among  the  third-shell  bonds  adjacent  to 
the  second-shell  atoms  are 

4x3x3^[2tAr1T3)2/(8d2)  +  (ArjT4)2]=4/3y2d2  . 

For  continuum  and  the  only  contribution  from  this  group 
are  the  0  terms  between  the  second-  and  third-shell  bonds. 
Since  the  displacements  in  the  continuum  areproportional 
to  l/R2,y‘  —  8 t/2y /( 1 9 vT9 )  and  y"  =  sV2y/(l\vTl). 
Thus,  these  ft  terms  become 

^0  ~-(38+  5y')2  +  ~-l8-b2y  —  5y")2 

-Hh^r 


4  0  terms  for  bonds  adjacent 
to  the  third-shell  atoms 

These  terms  only  enter  model  A,  so  r'  =  r"  —  0.  There 
are  two  different  groups,  one  like  those  adjacent  to  C’  and 
another  like  those  meeting  at  C".  The  four  bond  vectors 
pointing  away  from  C'  are 


r,=(-  1,  —  1  +y,  —  1  4  y)d/^3  , 
!;=-(  -  1,  1,  1W/V3  , 
r,-  ( I,  —  !,  1) d/\r3  , 


r4  —  ( 1 ,  1,  ~Ud/V3  . 

Thus,  the  only  contribution  from  this  group  is 
4  1  3/j[Air  r.ljVSd2  -  20y2d2  . 

The  four  bond  vectors  around  C"  are 


r,=(  —  1+r,— 1— y.-iw/vl, 
r2=<  —  1, 1,  ltf/v^, 
rj=(l,  -1+y,  l+ykf/v^  , 


r4=(l,  1, -lkf/t/3, 

which  only  results  in  the  first-order  term  A(r2-r})  =  \yd2. 
Thus,  the  group  contributes  to 

4  X  3  X  30[2(  Ar2T3)2]/8d2=4/?y2cf 2  , 

and  the  combined  contribution  from  these  two  groups  is 
60y2d\ 

APPENDIX  C:  CRITERION  OF 
MISCIBILITY 

Starting  with  Eq.  (24)  and  using  the  SSHS  model 
a  —  ai,  one  finds  the  mixing  enthalpy  parameter  ft  to  be 

Cl—\d{dAC—dBC)2N0  ,  (Cl) 

where  N0 .  is  Avagodro’s  number  and  a—  \(aAc+<*Bc)- 
Then  relate  the  mean-square  bond-length  fluctuation  ( f 2 ) 
at  the  melting  temperature  Tm  to  Tm  for  a  compound  by 
equating  the  average  potential  energy  per  unit  cell  to  half 
of  the  thermal  energy: 

<Fpol>~4x|a(f2>  =  |(2x3*,7'm),  (C2) 

where  kB  is  the  Boltzmann  constant.  Defining  a  Lieder- 
mann  ratio  of  melting  Xm  by 

((?))xn~Xmd  (C3) 

and  choosing  the  mixing  criterion  to  be  Tc/Tm  <  1,  where 
Tm  now  is  the  smaller  value  of  the  two  melting  tempera¬ 
tures  of  the  constituent  compounds,  we  require  that 

Te  (Ik  3  ^AC-dBC)2 


Tm  4  RtaX2md2  8  X2md2 


|  6o  i  /8m  <  1  , 


where  6m  =  1 . 63Ym  and  |8q|  is  the  percentage  bond- 
length  difference. 
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By  starting  from  realistic  band  structures  of  the  constituent  materials,  the  electronic  structure  of 
Si.Ge,  ,  alloys  are  obtained  in  the  coherent-potential  approximation  (CPA).  Various  quantities,  in¬ 
cluding  'he  bowing  parameter  of  the  fundamental  band  gap  and  the  energies  of  several  optical  gaps, 
the  masses,  and  the  linewidths  of  the  E0  and  £,  transitions,  are  calculated  on  the  basis  of  both  diag¬ 
onal  and  off-diagonal  CPA.  All  of  the  band-energy  and  linewidth  predictions  are  in  good  agreement 
with  experiments.  Furthermore,  the  theory  yields  an  alloy-scattering-limited  electron-drift  mobility 
in  qualitative  agreement  with  experimental  results. 
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I.  INTRODUCTION 

Semiconductor  alloys  offer  the  freedom  to  design  ma¬ 
terial  properties  by  choosing  appropriate  alloy  constitu¬ 
ents.  In  some  cases,  the  physical  properties  of  the  alloys 
can  be  quite  different  from  those  of  the  constituents.1-4 
In  recent  years,  there  has  been  a  renewed  interest  in 
Si,GC|_x  alloys’  and  superlattices.5- 10  Because  silicon  is 
the  most  technologically  advanced  semiconductor,  the  re¬ 
sults  on  Si-Ge  systems  have  many  potential  applications. 

The  lattice  constants  of  silicon  and  germanium  differ 
by  ~4%.  Hence,  the  strain  introduced  in  the  formation 
of  SixGe,_*  alloys  can  affect  the  band  structure10  and  the 
transport  properties.8  Prior  authors  used  virtual-crystal 
approximation  (VGA)  (Refs.  II  and  12)  and  coherent- 
potential  approximation  (CPA)  (Ref.  13)  to  study  the 
band  structure  and  related  properties.  Either  because  of 
less  accurate  band  structures  of  the  constituent  materials, 
or  because  of  the  approximations  involved  in  the  alloy 
formalism,  these  calculations  predicted  only  trends  of 
specific  quantities,  not  quantitatively  accurate  results.  Be¬ 
cause  the  r-state  site  potentials  ie,)  for  silicon  and  ger¬ 
manium  differ  by  approximately  1.5  eV,  VCA  cannot  ac¬ 
curately  describe  effective  masses  and  other  finer  details 
of  the  band  structure  Because  of  the  use  of  poor  basis 
functions,  earlier  CPA  work15  predicted  alloy  broadening 
of  conduction-band  states  substantially  differing  from  ex¬ 
periment  The  purpose  of  this  paper  is  to  correct  these 
flaws  and  treat  transport  phenomena. 

Because  of  a  substantial  difference  between  the  site  po- 
lentials  and  lattice  constants  of  silicon  and  germanium, 
we  incorporated  both  chemical  and  structural  disorder  in 
the  calculation  of  the  electronic  structure  of  SixGe|_x  al¬ 
loys.  Thus,  both  diagonal  and  off-diagonal  CPA  are  in¬ 
cluded  in  the  predicted  band  structure  and  relatied  quanti¬ 
ties  Parts  of  the  band  structure  have  been  used  to  study 
the  Si  2p  core  exciton14  and  the  alloy  mobilities.15  A 
comprehensive  report  of  the  calculations  and  results  is 
presented  here. 

The  rest  of  the  paper  is  arranged  as  follows.  The  de¬ 


tailed  procedure  of  fitting  silicon  and  germanium  band 
structures  is  given  in  Sec.  II.  The  VCA,  CPA,  and  off- 
diagonal  CPA  calculations  are  described  in  Sec.  III.  The 
results  and  interpretation  of  the  alloy  band  structures  and 
mobility  are  given  in  Sec.  IV. 

II.  BAND-STRUCTURE  BASIS 

In  order  to  derive  an  accurate  alloy  band  structure,  one 
must  start  from  a  realistic  band  structure  of  the  constitu¬ 
ent  materials.  Chen  and  Sher  have  developed  a  method,16 
following  a  prescription  of  Kane17  and  Chadi,18  which  in¬ 
cludes  all  long-range  interactions,  and  then  they  have 
fine-tuned  the  band  structure  with  an  adjustable  local 
Hamiltonian.  Because  the  details  have  already  been  pub¬ 
lished,1619  the  underlying  method  will  be  presented  here 
in  brief. 

Gaussian  orbitals  of  the  type  a  (a  can  be  s,  px,  py,  or 
p,)  for  each  sublattice  in  a  cell  are  used  to  construct  the 
corresponding  Bloch  basis.  In  this  basis  set,  the  overlap 
matrix  and  the  Hamiltonian  derived  from  empirical  pseu¬ 
dopotentials  can  be  calculated.17, 18  It  is  possible  to  cast 
the  problem  in  a  basis  set  of  Gaussian  orbitals  in  which, 
in  crystal  units  (c.u.),  the  same  exponential  factors  apply 
for  all  III-V  compounds.19  In  this  universal  basis,  the 
overlap  matrix  and  the  kinetic  energy  matrix  are  the  same 
for  all  III-V  compounds.  Then,  by  a  unitary  transforma¬ 
tion,  the  basis  set  is  orthonormalized.20  The  Hamiltonian 
in  this  new  basis  set  denoted  H0(k).  The  band  structure 
resulting  from  this  method  reproduces  the  results  of  ela¬ 
borate  band-structure  calculations  within  a  few  percent 
throughout  the  Brillouin  zone  (BZ).  To  establish  accu¬ 
rately  certain  important  band-structure  features  adjacent 
to  the  gap,  an  extra  small  8x8  Hamiltonian  matrix 
J/|( k)  is  added  to  //0(k).  This  f/|(k)  has  the  form  of  a 
tight-binding  (TB)  Hamiltonian,  in  which  only  the 
nearest-neighbor  interactions  are  included,  and  simulates 
the  effect  of  nonlocal  pseudopotentials  and  an  expanded 
orbital  set.  The  total  Hamiltonian  H(k)  in  this  orthonor¬ 
malized  basis  set  is  diagonalized  to  obtain  the  band  cner- 
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gies  and  the  corresponding  wave  functions. 

Following  this  procedure  with  the  same  exponential 
factor  0=0.26  in  the  Gaussian  orbitals  for  both  silicon 
and  germanium,  the  matrix  H0( k)  is  obtained.  For  sil¬ 
icon  and  germanium,  H  t  contains  six  adjustable  parame¬ 
ters:  namely,  the  corrections  to  the  term  values  A,  and  A, 
and  to  the  nearest-neighbor  interactions  Vn,  V9,  and 
Vv.  The  values  of  A,,  Ap,  Va,  and  V„  are  determined 
from  fitting  the  three  experimental  energy  gaps21-2*-39  at 
nk  =  0):  Tj-r |,  T |5-  r'25,  and  Tj-  r'25,  and  the  photoelec¬ 
tric  threshold  values  —  5.07  and  — 4.80  eV  for  silicon  and 
germanium,  respectively.29  The  remaining  parameters  Vv 
and  Vxy  are  obtained  from  the  experimental  values21  ~ 28  of 
the  gaps  X\c-Xio  and  L{c-L'io.  Some  adjustments  in 
these  input  quantities  are  made  to  obtain  an  overall  good 
band  structure  with  more  accurate  effective  masses. 
Table  I  lists  the  empirical  pseudopotential  form  factors 
and  the  parameters  used  to  obtain  the  band  structure. 
The  calculated  band  structures  and  experimental  values 
are  given  in  Table  II.  From  Table  II,  one  can  see  that  an 
excellent  fit  to  the  silicon  and  germanium  band  structure 
is  obtained:  All  the  calculated  values  lie  within  the  exper¬ 
imental  uncertainties.  The  optical  difference  between  L 
and  Lj„,  r'25  and  f  15  are  in  excellent  agreement  with  the 
known  optical  transition  values. 

Although  the  calculated  effective  transverse  masses 
agree  very  well  with  experiment,  the  effective  longitudinal 
mass  for  germanium  is  less  than  the  experimental  value. 
This  is  due  mainly  to  our  attempt  to  have  a  common  0 
and  the  choice  of  local  pseudopotentials,  causing  //0(k)  to 
be  the  same  in  crystal  units  (c.u.)  for  both  germanium  and 
silicon.  Because  of  the  common  //0(k),  the  alloy  disorder 
is  contained  in  these  adjusted  parameters.  This  H0  would 
also  be  useful  for  the  interface  and  superlattice31  prob¬ 
lems.  If  we  grant  ourselves  the  freedom  to  adjust  Vv, 
longitudinal  effective  mass  in  germanium  can  be  fitted  to 


TABLE  I.  Pseudopotential  form  factors  and  the  band  param¬ 
eters  (in  eV). 


Parameter 

Silicon 

Germanium 

vtVi) 

-2.872 

-2.872 

V(^4) 

0.124 

0.124 

v  v/i) 

0.638 

0.638 

v(VT\) 

0.109 

0.109 

A, 

-16.175 

- 16.922 

A, 

— 16.109 

-  14.971 

v* 

-0.111 

0.131 

v* 

0.040 

0.150 

0.025 

0.030 

v*, 

0.050 

0.100 

the  experimental  value.  When  is  changed,  the 
L  | e,L  j„  will  also  change.  We  have  chosen  not  to  do  this 
because  little  is  gained  for  the  extra  complexity.  For  an 
indirect-gap  semiconductor,  the  important  effective  mass 
used  in  transport  studies  is  the  conductivity  mass 

Because  m*t  »m*,  in  germanium,  m*  will  not  be  much 
different  if  a  less  accurate  value  of  mf*  is  used.  More¬ 
over,  the  SixGe,__,  alloys  which  have  potential  device  ap¬ 
plications  are  in  the  silicon-rich  region,  where  the  effec¬ 
tive  mass  at  the  L  edge  is  not  expected  to  affect  the  fur¬ 
ther  studies. 

It  is  important  to  note  that  an  excellent  fit  to  the  exper¬ 
imental  values  can  be  obtained  with  only  seven  adjustable 
parameters  VB,  Vv,  Vv),  with  0  being 

universal  in  c.u.  The  calculated  band  structure  of  silicon 
and  germanium  are  shown  in  Fig.  1(a)  and  Kb),  respec¬ 
tively.  The  characteristic  indirect  gaps  are  clearly  seen. 


TABLE  II.  Calculated  symmetry-point  band  energies  of  silicon  and  germanium  compared  with  ex¬ 
periments  and  empirical-pseudopotential-method  (EPM)  results,  (all  energies  in  eV). 


Silicon  Germanium 


Bands 

Calculated 

Expt.*-EPMh 

Calculated 

Expt.'-EPM' 

r, 

-12.60 

12.4±0.6 

-  12.56 

-12.610.3 

l!u 

-10.26 

—  9.310.4 

-10.74 

-10.610.5 

L  i„ 

-6.99 

—  6. 810.2 

-7.65 

-7.410.3 

-8.29 

-7.69 

-9.20 

-8.65 

*4, 

-2.55 

-2.86 

-2.55 

-3.29 

L 

-  1.11 

-1.210.2 

-1.13 

-1.110.2 

r2s. 

0.0 

0.0 

0.0 

0.0 

Lw 

2.24 

2.23 

0.76 

0.76 

r;. 

4.10 

4.0010.05 

0.99 

0.99 

r  iv 

3.43 

3.40 

3.24 

3.22 

Xu 

1.34 

1.17 

0.95 

1.16 

Lv 

4.34 

4.34 

4.16 

4.25 

E, 

1.11 

1.11 

0.76 

0.76 

(0.8,0,01 

(0.8,0,01 

10.5,0.5,0.5) 

(0.5,05,0.5) 

0  89 

0.91 

1  09 

1  59 

m* 

0.16 

0.19 

0.077 

0.082 

m* 

0.35 

0.50 

0.28 

0.34 

'References  21  28. 
Reference  30 
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ALLOY  CONCENTRATION.  * 
FIG.  1.  (Continued). 


These  band  structures  compare  favorably  with  the  best  re¬ 
sults  available,  and,  in  contrast  with  those  obtained  in  the 
usual  empirical  TB  approaches,  produce  good  conduction 
bands.  However,  if  they  were  included,  the  spin-orbit  in¬ 
teractions  (neglected  in  this  work)  would  add  some  fine 
details,  especially  the  splitting  between  the  heavy-  and 
light-hole  bands  near  T. 


scattering  potential  parameter,  the  difference  between  the 
j-electron  site  potentials  in  Si  and  Ge,  was  1.49  eV,  close 
to  our  value  1.46  eV,  the  calculation  predicted  too-large 
linewidths  in  the  E0  spectrum  and  essentially  no  effect  on 
the  electron  mobility.  With  the  availability  of  a  set  of 
goM  basis  functions  and  constituent  band  structures, 
more  realistic  band  structures  of  the  alloy  can  be  obtained. 

In  the  current  model,  we  have  a  TB  Hamiltonian, 
which  contains  matrix  elements  to  all  ranges.  The  sim¬ 
plest  alloy  model  is  to  assume  that  the  important  disorder 
resides  only  in  the  diagonal  matrix  elements  e,  and  (p.  In 
our  model,  the  ff1  and  ep*  differ  by  1.46  eV,  whereas  e£' 
and  ef*  differ  by  0.2 1  eV.  For  the  present,  we  neglect  the 
disorder  in  the  off-diagonal  element.  Mathematically,  we 
have 


tfnioy  =H+2vi>  (1) 

i 

where  /  is  a  fee  lattice  vector  identifying  a  site,  and  V (  is 
the  8x8  diagonal  matrix  with  elements  U,— 

Up—ep  —  ep  in  the  orthonormal  local  orbital  s| Ija);  j 
denotes  the  two  atoms  in  the  unit  cell  labeled  /,  a 
represents  s  or  p  symmetry,  and  ej  and  (p  are  the 
concentration-weighted  average  values  of  s  and  p  silicon 
and  germanium  term  value  energies. 

The  one-particle  alloy  Green’s  function  is  defined  as 


1  alloy 


(Z)» 


i 


Z-H , 


alloy 


(2) 


We  are  after  the  configuration  average  of  this  Green’s 
function,  which,  in  effective  medium  theory,  is  replaced 
by  an  effective  Green’s  function  G, 


G(Z)= - — -  , 

Z  -H-E(Z) 


(3) 


where  2  is  the  self-energy.  In  CPA,  we  can  now  write 
I  =  with  2/  being  an  8x8  matrix  in  the  basis 

!  Ija )  having  the  form 


III.  ALLOY  CALCULATION 
A.  VCA 

Because  we  have  the  same  H0{  k)  matrix  for  both  sil¬ 
icon  and  germanium,  it  is  only  the  //,  matrices  of  the 
constituents  in  scaled  VCA  which  distinguishes  them.  In 
this  approximation,  the  diagonal  elements  of  the  alloy 
Hamiltonian  Hi  It)  are  simply  the  concentration-weighted 
average  of  the  corresponding  elements  of  the  pure  silicon 
and  germanium  Hamiltonians,  whereas  the  off-diagonal 
elementsof  Hi  k)  are  obtained  by  assuming  a  1/d2  depen¬ 
dence  Hi  k)  can  be  diagonalized  to  obtain  the  VCA  band 
structure  for  various  concentrations  x.  The  VCA  band 
structures  for  x  =0. 1  and  0.5  are  shown  in  Fig.  1(c)  and 
1(d),  respectively. 


B.  Diagonal  CPA 

An  earlier  work  on  CPA  band  structure  of  SiGe  al¬ 
loys1'  is  based  on  a  local  but  energy -dependent  pseudopo- 
tential  approximation.  While  the  value  of  the  principal 


2  = 
where 


d  = 


d  Q 
0  A 

2,  0 
0  2, 
0  0 
0  0 


0  0 

0  0 

2,  0 

0  2p 


(4) 


Here  2,  and  2p  are  the  s  and  p  parts  of  the  self-energy. 
The  2,  and  2p  are  determined  from  the  conditions  that 
the  average  atomic  /  matrix  with  respect  to  the  CPA 
Green’s  function  G  is  zero.  With  our  ansatz  for  2,  the 
matrix  equation  {t )  =0  reduces  to  two  coupled  equations 
<t,  >=0  and  (tp)=0,  where  the  average  is  the 
concentration-weighted  average  (Q )  =xQs'  +yQa',  and 
the  r  is  defined  as 


^  =  (Gf-2a)[l-Fa((/f-20))-'  , 


( a  —  s  or  p  ,  0  =  Si  or  Ge  >  .  1 5 ) 
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In  the  above  expression  Fa  is  the  diagonal  matrix  element 
of  G  in  the  local  basis 


Fa{Z)s(lja\G(Z)\  ija  >  . 


I,  and  2p  are  coupled  because  F,  and  Fp  each  contain 
both  2,  and  2p. 

An  iterative  average-  f-matrix  (IATA)  procedure19  is 
employed  to  solve  the  CPA  equation.  This  procedure  im¬ 
proves  2a  upon  a  guessed  solution  2°  through  the  follow¬ 
ing  equations: 

2,=2?+(/10>(l+/f(/J°>)-1  , 

(6) 

2p  =  2“  +  <rp°)(l+f^(f°))-1  , 

where  ( t°  )  and  F°a  are  similar  to  those  in  Eq.  (5)  except 
that  2“  now  replaces  2a-  The  most  time-consuming  cal¬ 
culation  is  then  the  computation  of  the  local  Green’s 
functions  ff  and  F%,  given  by  the  BZ  summation,  e.g., 


F°(Z) 


_  1 

Z-H(  k)-2?  „ 


where  the  inverse  of  an  8  X  8  matrix  is  involved  for  every 
k  This  can  be  simplified  by  observing  that  2„  has  the 
same  form  as  2  in  Eqs.  (4)  and  (5)  and  that  the  4x4  ^ 
matrix  can  be  written  as  A  =  2pl  + (2,—  2p)£  where  1  is 
the  identity  matrix  and 


1  = 


10  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


Defining  the  matrix 

\l  Q 

Q  l 


e  =  (2?-  2°) 


F?(Z)  and  F°(Z)  can  now  be  calculated  from 


(7) 


F°<z>=  2{g22(k.Z)+gn(k,Z)+gu(k,Z)]  , 


where 

with 

u;a(k)u„^k) 

?  Z~e,(k)-2j 

In  Eq  (9),  r,(k)  is  the  band  energy  in  VCA  and  |  f/^lklj 
satisfy  the  following  Eigen  equation: 

Xf/^lk'f/^tkt ..--eJk)U„<k)  ■ 
e 


>-2 


(8) 


(9) 


Because  the  o  matnx  has  only  two  nonzero  elements,  the 
matrix  inversion  in  Eq.  !8'  is  obtained  analytically. 

A  substantial  reduction  in  computer  time  is  made  possi¬ 


ble  by  using  an  analytical  continuation  method.32  In  this 
method,  2,  and  2p  are  calculated  as  a  function  of  com¬ 
plex  Z,  and  then,  using  the  analytical  properties  of  the 
self-energy  and  Green’s  functions,  they  are  interpolated 
for  real  Z.  Because  the  functions  2f,  2p,  and  G  are 
smooth  for  complex  Z,  the  CPA  iterations  and  BZ  in¬ 
tegrations  can  be  carried  out  with  substantially  less  com¬ 
puter  time. 

For  the  concentrations  x  —0. 10  and  0.50,  the  L  and  X 
(A)  gaps,  respectively,  are  preferred.  The  L  to  X  (A) 
crossover  takes  place  near  x  =0. 15.  The  CPA  correction 
to  L  and  X  edges  at  x  =0. 10,  0.15,  and  0.50  should  be 
good  enough  to  study  the  quantitative  variation  of  band 
gap  in  Si1Gel_x  alloys.  Hence,  the  calculations  are  car¬ 
ried  out  for  these  three  cases.  In  addition,  because  the  ex¬ 
perimental  results  are  available  for  x  =0. 109,  CPA  calcu¬ 
lations  are  also  done  here  for  comparison.  As  expected, 
2,  is  much  larger  than  2p  for  all  the  cases.  The  self¬ 
energy  as  a  function  of  energy  is  plotted  in  Fig.  2  for  an 
x  =0.50  alloy. 


C.  Off-diagonal  CPA 

As  mentioned  earlier,  silicon  and  germanium  differ  in 
their  lattice  constant  by  —4%.  In  order  to  include  the  ef¬ 
fect  of  the  structural  disorder,  the  CPA  calculation  is  re¬ 
peated  next  with  off-diagonal  (OD)  disorder  included.  By 
an  application  of  the  molecular  coherent-potential  approx¬ 
imation  (MCPA),33  Hass  el  al.  included  OD  disorder  in 
the  CPA  calculation  of  A'x  A  \'_XB  semiconductor  alloys.34 
Assuming  that  B  atoms  occupy  the  sites  of  an  ordered 
zinc-blende  virtual  lattice,  they  modeled  the  dominant 
structural  effect  as  the  difference  in  A'-B  and  A"-B  hop¬ 
ping  matrix  elements.  Hence,  the  chemical  and  structural 
disorder  effects  are  treated  as  random  variations  of 

e*,  V\,Vi  ,  where  the  symbols  have  their  usual  mean- 

■  55 

ing. 

The  extension  of  the  method  to  Si,Gei_x  alloys  is  not 
straightforward,  mainly  because  silicon  and  germanium 
can  occupy  both  sublattices;  hence,  there  can  be  no  or¬ 
dered  virtual  lattice  in  this  case.  If  we  choose  the 
tetrahedral  unit  cell  as  the  molecular  unit  for  MCPA,  we 
see  that  the  disorder  is  not  cell  diagonal.  However,  by 
choosing  an  appropriate  basis  set,  we  can  make  the  inter- 
cell  interaction  be  the  highest-order  effect.  We  start  with 
a  hybrid  basis  |  Ih  ),■  obtained  from  the  spy  hybrid  orbi¬ 
tals.35  The  hybrids  1  through  4  (i  —  1— 4)  are  obtained 
from  orbitals  centered  on  a  sublattice  I  site,  and  the  states 
5  through  8  (i  =5—8)  are  those  from  the  orbitals  located 
at  the  four  nearest-neighbor  sites  on  sublattice  II.  The 
Bloch  basis  states,  corresponding  to  AX,T2  symmetries, 
j  k)/  located  on  a  I  site  (i  =  1—4)  and  II  site  (/  =5—8),  are 
obtained  from  the  corresponding  hybrid  states  given  by 
the  relation 

|/>i  =  2Cij\lh)i  ,  (10) 

where 


c= 


£■ 

Q 


0 

£t 
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FIG.  2.  Variation  of  the  imaginary  part  of  the  self-energy  2,(CPA)  as  a  function  of  energy  for  x  =0.50. 


and 

1111 
j  1  1  -1-1 

£l=  2  1-1  1  -1 
1-1-1  1 

An  explicit  definition  of  these  orbitals  can  be  found  in 
Ref.  20.  In  this  new  basis,  the  self-energy  2  at  the  given 
site  takes  the  form 

Co  Cl 

^  ei  ?o  ’ 
where 

2,  0  0  0 

0  2,  0  0 

2o=  0  0  2,  0  ’ 

0  0  0  2, 

2;  0  0  0 

0  2,  0  0 

e2=  o  o  2;  o 

0002; 

2;  0  0  0 

0  2;  0  0 

2(1  “  0  0  2;  0 

0002; 


and 

2;  =  1132,-1-2,),  (13) 

where  (2,, 2, ),(2j,2;)  are,  respectively,  the  self-energies 
associated  with  diagonal  and  off-diagonal  disorder  corre¬ 
sponding  s  and  p  symmetries. 

The  self-energies  can  be  obtained  again  from  the  I  AT  A 
iteration  procedure: 

2  =  20+«7'»(l+F«r»r'  ,  (14) 


with 

lS(k)  +  2]0=£:(Z,k){?  , 

{{T))=x(TA)+y(TB)  , 

and 

<  TA)=x*tft  t-Axiyt^a+6x2y1tAlBj 

+4 pVZb, -hvV,  .  >4=Si  (16) 

with  a  similar  expression  for  ( TB ).  Physically,  for  a 
given  A  atom  at  the  center,  the  other  four  atoms  in  the 
molecular  unit  cell  can  be  all  A  atoms,  three  A  atoms  and 
one  B  atom,  two  of  each,  one  A  atom  and  three  B  atoms, 
or  all  four  can  be  B  atoms.  (  TA  )  or  <  TB )  represents  the 
configuration-averaged  t  matrices,  and  (( T )>  is  the 
concentration-weighted  average  of  the  configuration.  By 
exploiting  the  symmetry,  as  seen  in  Eq.  (11),  one  can 
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reduce  this  problem  to  solving  two  2x2  coupled  matrix 
equations.  Equation  (14)  can  be  iterated  to  obtain  2,,  1f, 
i;,  and  2P  After  every  iteration,  we  get  a  new  set  of  2,, 
I  2;,  2'p,  and  2*:  The  new  set  has  not  been  tested  to 
see  if  Tk  is  still  given  by  Eq.  (13).  In  our  calculation,  we 
did  not  iterate  to  obtain  a  new  2'* ;  the  error  introduced  by 
this  approximation  is  expected  to  be  very  small.  As  in  the 
case  of  diagonal  CPA,  the ‘computation  can  be  substan¬ 
tially  reduced  by  the  method  of  analytical  continuation.32 

IV.  DISCUSSION 
A  E0  and  E{  optical  transitions 

The  VCA  values  of  E0  (r^-Tjj,,)  and  E'0  (r13<.-r2j[,) 
and  their  measured  values  are  plotted  as  a  function  of  x  in 
Fig.  1(e).  Because  the  measurements36  are  made  at  room 
temperature,  the  experimental  values  are  smaller  than  the 
values  calculated  from  the  zero-temperature  band  struc¬ 
ture.  Inclusion  of  the  relativistic  effects,  which  are  not 
present  in  our  calculations,  is  expected  to  form  a  more  ac¬ 
curate  basis  for  comparison  with  the  experiments.  As 
seen  from  Fig.  1(e),  the  theoretical  and  the  experimental 
values  both  have  a  linear  variation  with  x.  Similar  calcu¬ 
lations  of  E |  —  L'iv)  also  have  a  linear  variation  on  x 

and  are  in  qualitative  agreement  with  experiments.36 

From  the  CPA  self-energies  2,  and  2p,  it  is  straightfor¬ 
ward  to  calculate  the  correction  to  the  VCA  bands.  The 
calculated  complex  band  structure  is  plotted  for  x  =0.50 
in  Fig.  3.  The  CPA  corrections  are  shown  only  in  the  vi¬ 
cinity  of  the  band  gap.  'the  shaded  portion  represents  the 
half-width  of  that  energy  state.  Because  s  scattering  is 
dominant  in  these  alloys,  we  see  that  the  major  disorder 
lies  in  the  conduction  band.  The  topmost  valence  band, 
with  its  rich  p  content,  is  least  affected.  The  CPA  band 
structure  is  used  to  calculate  the  E0  and  Ex  peak  posi¬ 
tions  for  x  =0.10,  0.109,  0.15,  and  0.50  concentrations. 
The  calculations  and  the  data  from  Ref.  36  show  a  small 
bowing  that  is  not  seen  on  the  scale  of  Fig.  1(e). 


FIG  3.  Calculated  CPA  complex  band  structure  of  the 
St„  .(Je,. .  alloy.  Only  the  bands  in  the  vicinity  of  the  energy  gap 
are  shown  The  shaded  portion  represents  the  alloy  broadening 


-4.3  -4.0  -3.5  -3.0 
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FIG.  4.  Imaginary  part  of  2,(CPA)  and  2,(MCPA)  as  a 
function  of  energy  for  x  =0. 50  alloy. 

The  self-energies  2,,  2p,  2'p,  and  2j  are  calculated  in 
MCPA  for  the  x  =0.50  alloy.  As  in  the  case  of  CPA,  the 
self-energies  associated  with  s  symmetry  are  much  larger 
than  the  ones  associated  with  the  p  symmetry.  While  2', 
is  found  to  be  very  small,  2p  is  at  least  an  order  of  magni¬ 
tude  smaller — almost  zero.  However,  the  Im2,  obtained 
by  CPA  and  MCPA  differ  considerably.  As  seen  from 
Fig.  4,  the  difference  increases  as  one  goes  away  from  the 
band  edge.  Therefore,  the  lifetime  associated  with  the  al¬ 
loy  disorder  is  decreased  by  the  inclusion  of  OD  structural 
disorder.  In  addition,  the  OD  disorder  lowers  the  conduc¬ 
tion  band,  introducing  an  extra  bowing.  The  E0  and  E{ 


TABLE  III.  Calculated  values  of  £0,£i  and  their  respective 
half-widths  A(£0)  and  A(£|)  (all  energies  are  in  eV). 


X 

Quantity 

VCA 

CPA 

MCPA 

0.10 

Eo 

1.290 

1.248 

A<£0) 

0.011 

E , 

2.016 

1.995 

A(£0) 

0.001 

0.109 

Eo 

2.028 

2.009 

A  <£0> 

0.013 

E, 

1.318 

1.275 

A(  £,) 

0.002 

0.15 

E0 

1.442 

1.382 

A  <£0) 

0.032 

Et 

2.083 

2.051 

A  (£,) 

0.002 

0.50 

Eu 

2  517 

2.418 

2.391 

A<£„) 

0.186 

0.206 

£. 

2  578 

2  510 

2.498 

Ai  £.) 

0  0308 

0.0319 
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values  are  reduced  by  27  and  12  raeV,  respectively.  The 
VCA,  CPA.  and  MCPA  values  of  E0  and  £(  are  listed  in 

Table  III 

The  half-width  of  the  alloy  states  is  calculated  from  the 
imaginary  part  of  the  CPA  self-energies.  The  half-width 
of  the  lowest-lying  conduction  band  of  Si0  5Geo  5  alloy  is 
plotted  in  Fig.  5  as  a  function  of  Kx  in  the  [100]  direc¬ 
tion.  The  calculated  half- width  is  186  meV  for  the 
state  and  decreases  to  zero  at  the  band  edge.  Because  of 
the  negligible  alloy  broadening  of  the  topmost  valence- 
band  state,  the  half-width  corresponding  to  the  E0  transi¬ 
tion.  A(£ 0),  is  186  meV,  which  is  approximately  one-half 
of  the  previously  published  CPA  results.13  The  CPA 
value  of  the  half-width  corresponding  to  the  £,  transi¬ 
tion,  A(£|),  is  31  meV.  Because  of  the  increase  in  the 
imaginary  part  of  the  self-energies,  the  MCPA  values  of 
the  half-widths  of  the  E0  and  £(  transitions  are  206  and 
32  meV,  respectively.  Because  the  complete  £0  peak  is 
not  shown  in  the  published  electroreflectance  spectrum,34 
it  is  difficult  to  estimate  the  corresponding  half-width. 
However,  one  can  conclude  from  the  spectrum  of  the 
x  =0.458  alloy  that  the  half-width  of  the  Ex  transition  is 
considerably  smaller  (~50  meV)  than  that  of  the  £0  tran¬ 
sition.  The  agreement  between  the  experimental  and  the 
theoretical  values  can  be  regarded  as  good  because  there 
are  errors  in  estimating  the  width  from  the  published 
spectra,  and  we  have  neglected  the  extrinsic  broadening 
due  to  the  apparatus  used  in  the  experiments. 

In  order  to  make  a  more  accurate  comparison  with  the 
experiments,  the  CPA  values  of  A(£0)  and  &(E\ )  are  cal¬ 
culated  for  the  x  =0. 109  alloy.  The  calculated  half¬ 
widths  of  the  E0  and  £,  transitions  are  13  and  2  meV, 
respectively.  From  the  spectrum,  we  estimate  the  corre¬ 
sponding  values  to  be  8—15  and  3—6  meV.  We  see  that 
CPA  values  are  in  excellent  agreement  with  these  experi¬ 
ments.  Because  x  is  small,  the  inclusion  of  off-diagonal 
disorder  is  not  expected  to  change  the  calculated  values 
significantly. 


200  f~  ■ — t "  . . r~  ■  ~ t  - i 


K, 


FIG  5  Variation  in  the  width  of  the  lowest-lying  conduction 
band  as  a  function  of  K,  in  the  (KX)]  direction  for  the  x  =0  50 
alloy 


B.  Energy  gap 

The  fundamental  gaps  of  these  alloys  are  calculated  as  a 
function  of  concentration.  The  VCA  gap  is  an  increasing 
function  of  x  with  a  slope  discontinuity  at  x~0. 11.  The 
conduction-band  minimum  changes  from  L  point  to  X 
(A)  point  at  this  crossover.  In  addition  to  the  band  gap, 
the  effective  electron  masses  and  the  band  edge  K0  are 
also  calculated.  When  X  (A)  gap  is  preferred,  the  band 
edge  moves  linearly  from  k  at  (0.9,0,01,  «o  is  to 
(0.8,0,01,  =  1.  The  effective  masses  at  a  given  minimum  in¬ 
crease  linearly  from  their  pure  germanium  values  to  the 
corresponding  pure  silicon  values. 

Using  CPA  self-energies,  the  band  gap,  band  masses, 
and  the  band  edge  are  also  calculated.  The  position  of  the 
band  minimum  did  not  change  by  virtue  of  the  inclusion 
of  off-diagonal  disorder.  While  the  effective  transverse 
mass  remains  almost  the  same  as  the  VCA  value,  the 
longitudinal  mass  has  a  maximum  of  12%  enhancement. 
Because  the  real  part  of  CPA  self-energies  is  negative  in 
the  forbidden  gap  region,  an  extra  bowing  is  introduced  to 
the  VCA  energy  gap.  Because  of  this  bowing,  the  L-X  (A) 
crossover  takes  place  near  x~0. 13.  The  VCA,  CPA,  and 
experimental36  bowing  parameters  are  0.06,  0.18,  and  0.24, 
respectively.  The  calculated  energy  gap  is  plotted  as  a 
function  of  x  in  Fig.  6. 

Because  of  the  negligible  change  in  the  effective  masses, 
the  corresponding  values  in  the  pure  materials  are  used  in 
the  calculation  of  the  alloy-scattering-limited  electron  mo¬ 
bility.  The  CPA  -Y-gap  Eg  and  L-gap  Eg  are  fitted  to  a 
polynomial  form.  The  generalized  Brooks’s  formula  that 
is  applicable  to  the  alloys  with  an  indirect  gap  and  multi¬ 
ple  bands  is  used.15  The  calculated  electron-drift  mobility 
and  the  experimental  Hall  mobility1  are  plotted  in  Fig.  7, 
where  the  theory  explains  the  qualitative  behavior  of  ex¬ 
perimental  results.1  As  observed,5  even  a  few  atomic  per¬ 
cent  alloy  concentration  can  reduce  the  drift  mobility  sub¬ 
stantially.  It  can  be  seen  that  the  rate  of  decrease  near 
x  =0  and  jc  =  1  are  quite  different.  This  is  because  the  L 
edge  has  more  s  content  than  the  X  edge.  Because  the  s 


X 


FIG  6.  Variation  of  the  VCA  energy  gap  (dash-dolled  line) 
and  the  CPA  energy  gap  (solid  line)  as  a  function  of  x. 
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FIG  7  Calculated  drift  mobility  (solid  line)  and  the  experi¬ 
mental  (dashed  line)  Hal!  mobility  (Ref.  1)  as  a  function  of  x. 


scattering  is  dominant  in  these  alloys,  the  L  electrons  are 
scattered  more  than  the  X  electrons.  Precisely  for  this 
reason,  one  observes  a  dip  in  the  mobility  near  the  L  to  X 
(A)  crossover.  For  x  <0.13,  the  minimum  gap  is  the  L 
gap  After  the  crossover,  the  minimum  gap  is  the  X  (A) 
gap,  and  the  reduced  alloy  scattering  increases  the  average 
mobility.  For  still  larger  x,  the  mobility  decreases  because 
of  the  increased  alloy  disorder.  All  these  features  are 
clearly  seen  in  Fig.  7.  While  our  calculations  include  the 
intervalley  scattering  mediated  by  alloy  disorder,  the  ef¬ 
fect  of  other  scattering  mechanisms  is  expected  to  increase 
the  dip  near  the  crossover. 

The  calculated  alloy  scattering  rate  for  the  holes  is 
several  orders  smaller  than  that  for  electrons,  because  (1) 
the  valence  band  edge  has  dominant  p  content,  (2)  the  p- 
scattering  parameter  (Aep—  0.21)  is  only  y  of  Af(,  which 
alone  decreases  the  scattering  rate  for  holes  by  a  factor  of 
50,  and  (3)  finally,  the  imaginary  part  of  the  self-energy  is 
proportional  to  the  density  of  states,  which  approaches 
zero  at  the  band  edge.  Hence,  the  hole  mobility  in  this 
system  is  insensitive  to  alloy  disorder 

In  MCPA,  the  conduction  band  is  pushed  down,  be¬ 
cause  of  an  increase  in  the  imaginary  part  of  the  self- 
energy,  giving  rise  to  an  additional  bowing  in  the  funda¬ 
mental  gap  For  an  *=0.50  alloy,  the  gap  is  reduced  by 
7  meV.  The  bowing  parameter,  including  the  MCPA 
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correction,  is  0.21,  which  is  in  excellent  agreement  with 
experiment.37 

It  is  interesting  to  compare  the  results  of  our  calcula¬ 
tions  with  those  of  Hass  et  al. 34  In  their  calculations  on 
the  Ga)_,In,As  alloy,  CPA  introduced  an  extra  bowing 
in  the  fundamental  gap.  However,  after  the  MCPA 
corrections,  the  total  scattering  was  diminished  and  the 
results  were  similar  to  VCA  results.  These  results  were 
explained  in  terms  of  the  relative  strength  and  sign  of  the 
atomic  term  values  and  V We  extend  their  argument 
to  SijGe,-*  alloys.  The  hybrid  level  of  silicon  is  higher 
than  that  of  germanium.  Because  of  its  shorter  bond 
length,  the  V2  of  silicon  is  larger  than  that  of  germanium. 
Thus,  in  this  case,  both  effects  combine  to  give  more  dis¬ 
order  in  the  conduction  and  valence  bands.  Therefore,  the 
scattering  is  enhanced  in  these  alloys.  This  explains  the 
increase  in  the  imaginary  part  of  the  self-energy  due  to  in¬ 
clusion  of  OD  disorder  in  our  calculation. 

In  conclusion,  we  have  incorporated  both  chemical  and 
structural  disorder  into  the  calculation  of  the  CPA  band 
structure  of  Si^Ge^,  alloys.  The  calculation,  based  on  a 
realistic  band  structure  of  silicon  and  germanium,  sug¬ 
gests  that  the  band  gap  is  an  increasing  function  of  x  with 
a  slope  discontinuity  at  x~K).  13.  The  linewidths  of  the 
£0  and  Ex  transitions  calculated  by  CPA  and  MCPA  are 
in  good  agreement  with  experiments.  The  effects  of 
structural  disorder  on  diagonal  CPA  for  the  Ge-Si  alloy 
system  will  be  difficult  to  test  experimentally.  MCPA  de¬ 
creases  the  band  gap  only  slightly.  However,  unlike  the 
GalnAs  case,  the  imaginary  part  of  the  s  self-energy  in¬ 
creases  with  increasing  energy  and  serves  as  a  warning 
that  all  systems  will  not  exhibit  the  same  behavior. 
Therefore,  in  Ge-Si  alloys,  the  place  to  look  for  experi¬ 
mentally  significant  differences  between  CPA  and  MCPA 
predictions  is  in  the  high-lying  transitions.  However,  even 
in  the  £0  transition,  the  MCPA  linewidth  at  x  =0.5  is 
still  only  10%  larger  than  the  corresponding  CPA  value. 
The  calculated  ailoy-scattering-limited  electron-drift  mo¬ 
bility  is  in  qualitative  agreement  with  the  observed  Hall 
mobilities. 


ACKNOWLEDGMENTS 

This  work  was  supported  in  part  by  U.S.  Air  Force  Of¬ 
fice  of  Scientific  Research  Grant  No.  AFOSR-84-0282 
and  Defense  Advanced  Research  Projects  Agency 
(DARPAI/AFOSR  Project  No.  PR:FQ8671-851 100.  One 
of  us  (A ,-B.C.)  would  like  to  thank  Professor  W.  E.  Spie  r 
for  his  hospitality  at  Stanford  University. 


'H  M  Manasevit,  I  S  Gergis,  and  A  B  Jones,  Appl  Phys. 

Lett.  41.  464  (1982);  J  F.lcctron.  Mater.  12,637(1983). 

"S.  Krishnamurthy  and  J  A  Monarty,  Superlat  Micros! met.  I, 
209  (1985  .  Phys  Rev  B  32.  1027  (1985) 

J  C  Bean.  1  C  Feldman.  A  T  Flory.  S  Nakahara.  and  F  K. 

Robinson.  J  Vac  Sci  Trchnol  A  2.  436  Il984i 
HR  People.  J  C  Bean.  [)  V  Fang.  A  M  Sergent.  II  F  Storm- 


33 


BAND  STRUCTURES  OF  Si.Ge,.,  ALLOYS 


1035 


er,  K  W  Wechl,  R  T.  Lynch,  and  K  Baldwin.  Appl.  Phys. 
Lett  45.  1231  (19841. 

9F  Cerdetra,  A.  Pinczuk,  and  J.  C.  Bean,  Phys.  Rev.  B  31,  1202 
(1985). 

10T.  P.  Pearsall,  F  H.  Poliak,  and  J.  C.  Bean,  Bull.  Am.  Phys. 
Soc.  30.  266  (1985). 

nK.  E.  Newman  and  J.  D.  Dow,  Phys.  Rev.  B  30,  1929  (1984). 
IJM.  Z.  Huang  and  W.  Y.  Ching,  Superlat.  Microstruct.  1,  137 
(1985).  , 

l3D.  Stroud  and  H.  Ehrenreich,  Phys.  Rev.  B  2,  3197  (1970). 

I4S.  Krishnamurthy,  A.  Sher,  and  A.-B.  Chen,  Phys.  Rev.  Lett. 
55,  320(1985). 

>»S  Krishnamurthy,  A.  Sher,  and  A.-B.  Chen,  Appl.  Phys.  Lett. 
47.160(1985). 

I6A  B.  Chen  and  A.  Sher,  Phys.  Rev.  B  22,  3886(1980). 

,7E.  O.  Kane,  Phys.  Rev.  B  13.  3478  (1976). 

I8D.  J  Chadi,  Phys.  Rev.  B  16.  3572  (1977). 

'’A.-B  Chen  and  A.  Sher,  Phys.  Rev.  B  23,  5360  (1981). 

20A.-B.  Chen  and  A.  Sher,  Phys.  Rev.  B  26,  6603  (1982). 

2IR.  R.  L.  Zucca.  J.  P.  Walter,  Y.  R.  Shen,  and  M.  L.  Cohen, 
Solid  State  Commun.  8,  627  (19701. 

22R  R  L.  Zucca  and  Y.  R.  Shen,  Phys.  Rev.  B  1,  2668  (1970). 
23M.  Welkowsky  and  R  Braunstein,  Phys.  Rev.  B  5,  497  (1972). 
24R  A  Poliak,  L.  Ley,  S.  Kowalcyzk,  D.  A.  Shirley,  J.  Joanno- 
polos,  D.  J.  Chadi,  and  M.  L.  Cohen,  Phys.  Rev.  Lett.  29, 
1103  11973). 


2,L.  Ley,  R.  A.  Poliak,  F  R  McFeely,  S.  P  Kowalcyzk,  and  D. 
A.  Shirley,  Phys.  Rev  B  9,  600  (1974). 

26W.  D.  Grobman  and  D.  E.  Eastman,  Phys.  Rev.  Lett.  29, 
1508(1972). 

27D.  E.  Eastman,  W.  D.  Grobman,  3.  L.  Freeouf,  and  M.  Erbu- 
dak,  Phys.  Rev.  B  9,  3473  (1974). 

28W.  E.  Spicer  and  R.  C.  Eden,  Proceedings  of  the  Ninth  Inter¬ 
national  Conference  on  the  Physics  of  Semiconductors,  1968 
(Nauka,  Leningrad,  1969),  Vol.  1,  p.  61. 

NG.  W.  Gobeli  and  F.  G.  Allen,  Phys.  Rev.  245,  A137  (1965). 

R.  Cheliskowsky  and  M.  L.  Cohen,  Phys.  Rev.  B  14,  556 
(1976),  and  references  therein. 

31  J.  A.  Monarty  and  S.  Krishnamurthy,  J.  Appl.  Phys.  54,  1892 
(1983). 

32K.  C.  Hass,  B.  Velicky,  and  H.  Ehrenreich,  Phys.  Rev.  B  29, 
3697(1984). 

33F.  Ducastelle,  J.  Phys.  C  7,  1795  (1974). 

^K.  C.  Hass,  R.  J.  Lampert,  and  H.  Ehrenreich,  Phys.  Rev. 
Lett.  52,  77  (1984). 

35W  A.  Harrison,  Electronic  Structure  and  the  Properties  of 
Solids  (Freeman,  San  Francisco,  1980),  p.  65. 

36J.  S  Kline,  F.  H.  Poliak,  and  M.  Cardona,  Helv.  Phys.  Acta 
41,968(1968). 

37R.  Braunstein,  A.  R.  Moore,  and  F.  Herman,  Phys.  Rev.  109, 
695(1958). 


